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ABSTRACT 


The transient behavior of a generic military capability under wartime 
environment was analyzed and, under certain assumptions, a Wartime Sustainability 
Model (WSM) was developed analytically using various stochastic and inventory 
techniques. A simulation of the WSM was also developed to incorporate variations 
for repair policies such as repair prioritization and limited repair capability. These 
variations are extremely difficult to model analytically. The adequacy of these two 
models was verified using a numerical example. Finally, the feasibility of using 
OPUS-8, a steady-state spares optimization model developed by Systecon AB, 
Stockholm, Sweden, as an approximation to the analytical version of the WSM for 


the case of no repair capability was also investigated. 


lil 


TABLE OF CONTENTS 


I. INTRODUCTION 223 r. . «os pene el iete et cue i ner ene ene 1 
A. MOTIVATION FORTHE THESIS@. 7. os... 7 eae 1 
B. GOATS one ee te a Onn te nee ee oo 2 
C. PREVIEW . 2.2.2... ccs eee eae e eee renee: ==. ian 5 


Sete ee eee 8 8 tO eBEE Sle RUDMMEST ECT ae tare’. eRe s,s. eee 5 
A. SUSTAINABILITY DEFINITION =... ...... 5 
B.. POLICIES... ..c.5:4 SRR he cece eee as et 6 

1, War and Mobilization Plan(WMP). ................2008% 6 
2. Optimization of Investment In Spares ................... 8 
3. The Inadequacy of A Minimum Buy Policy ............... 8 
4. Cannibalization Roheys...... . : Re... 2G ee 9 
>. Repair Pnontization"Policies.. .. <7 ee eee 10 
C. MEASURES OF EFFECTIVENESS (MQEs) ..............- 10 
1. Mission-oriented MOEs (MMOEsS) ............. | ae ate 11 
2. System-oriented MOEs (SMOEs) ..............-.-008- 12 


Dro oS VEVeS TUG AND LOGISTICS SUPPORT ........ 14 


HI. AN ANALYTICAL WARTIME SUSTAINABILITY MODEL (WSM) . 15 


A. ANALYTICAL MODEL ASSUMPTIONS ................. 15 
Be RATE RE wR AACA Se ak ol nee en i ON 16 
CG DEMAND SIBORSSP ARE Serie Oe A ve ea sa eas 16 
1. Nonhomogeneous Poisson Assumption.................. 17 

a. The Distribution of LRU Failure Arrivals............ 17 

b. The Distribution of the Repair Process.............. se 


c. The Distribution of Failed Units of an LRU Still in 


BSc [Oc Teer te Wo veer em Rie EU EUMNGDs ss sss a a ee i 
2. Nonhomogeneous Compound Poisson Assumption ......... 20 
a. The Distribution of LRU Failure Arrivals ............ 20 


b. The Distribution of Failed Units of an LRU Still in 


RSC TAA acd Pe el Cho ho ire ere ene D ESSERE Sy pl 

ir MULTIPERINDENTURE PIPELINES ..........0. 0.6... ge 
1, A Case of Two-Indenture System ...................06. ZS 
aeeecanmibalization PolicyeccGe sw. ee ee 24 

Dee NOLCannIDalizallOn: sameness 425 es Se es ee ee ewe es 26 

EE. MULTIPEE-ECHELONIFRIPELINES seme eeeeeieseme:... - 28 
1. A Case of Two-Echelon Logistics Support ............... 28 

dee MNGeal RGD alieICI De UTS ©. x cueiiet en en ea aati 5 oa) 3 29 


(1) Nonhomogeneous Poisson Pipelines ............ 29 


(2) Nonhomogeneous Compound Poisson Pipelines .... 30 


b. “Total Pipeline). ties ee 31 

F. THE SMOES oc). cere regres le ee ea 36 
1. EB,(t), Expected Number of Backorders for LRU, ......... 3H 

a. EB,(t) for Nonhomogeneous Poisson Pipelines ........ a7 


b. EB (t) for Nonhomogeneous Compound Poisson 


PIPEHIMES Lee ig. eee or Sins eae eee ee Chao a tree 38 

2. EB(t), the Total Backorders of the Deployed Systems ...... 38 

3. OperationalAvailability (A(t) emer ene «ae oe 39 

a. Ao(t) without Cannibalization ................00.. 39 

b. AO) With Cannibalizadom wo sss ene 41 

G. THE MM. er se ene nce ane een” ot 44 
1 EINIMIG SP on ne 6 Oe eee soca ee eee 44 

2 PCENMCS <9NMCGS) tt ne eee ee 44 

H. LIMITATIONS OF THE ANALYTICAL APPROACH ........ 45 


IV. A SIMULATION VERSION OF THE WARTIME SUSTAINABILITY 


MOE rec cc ee ee Coe 0 ee Tee ee arama ect 46 
A. REASON FOR THE SIMULATION APPROACH ........... 46 
B. MODDED DESIGN ete a eet mn. MMP etre cs. 46 

1. LRU Failuré*Arrivals . Seer ee) oat 5. re 47 


5 PaitinemiwectMicanonme. ........... newer oey..... 48 

Ae Ne Dalla LOWS ce eee eee eM, ke ee 49 

aye ase We paiimem want nett. we ee 50 

BEI POtAINCIOER ara mreetite UN ic) 92 aidniiierslossieeente swe s 50 

Se Keplications tof statistical Significance” 2. 22.22... ee ce. 50 

6. Random Number Streams and Generation .............. 51 

Ua pe RepaimmleMOnitizatlon POlCleS sa ene eee as oso. ee a 

a.) pEinst @ome First semved (FGFS) Guawye........... 51 

b. Least Available Item First (LAIF) ................. a4 

C. IMPLEMENTATION OF WSM USING MODSIM ........... a2 
ees A LIS MIGAIINEERENCES. 2255 eeeen erie es eee 54 
1. The Resulting Distribution for Failure Arrivals ........... 54 

2. Replications for Statistical Confidence ...............04. a7) 

V. WSM BY AN OPUS-8 APPROXIMATION ................0-008- 58 
A. THE STEADY-STATE ASSUMPTIONS OF OPUS-8 ......... 58 
1. Current Features and Assumptions .................... 58 

2. Limitations Of Steady-state Assumptions................. 61 

B. OPUS-8’S SUSTAINABILITY OPTION .................... 61 


VI. VERIFICATION OF THE WARTIME SUSTAINABILITY MODELS . 63 


A. THE EXAMPLE ..........2. See 64 
1. War and Mobilization Plan (WPM) \23 eee ee oe 64 
2. System Structure and Logistic Support .................. 65 
B. RESULTS FROM ANALYTICAL CALCULATIONS ......... 67 
1. Computing the Demand Rates and Repair Flows ......... 68 
2. Calculating the Base Kepaimmipenine ee 69 
3. Caleulatme the Total Depot Piscine ser ee 70 
4. Calculating thediotalePipehneseny see... = ee 1s 
5. Calculating Expected Number of Backorders ............. ie 
6. Calculating Operational Availability.................... 73 
7. Calculating ye NMGS (0) ag i en We 
8.  Calcilanme*P(ENMGS*< NIMES)... oe eee os ee 74 
9. Computations for Analytical Case Two ................. 74 
C. RESULTS FROM SIMULATION RUNS ................... 1D 
1. ~eSunulation Variation'@ne™.., [a eee 76 
2. Simulation Variation#fWwoe.. 2.2] )-.cee eee. ae Ti 
3... Simulation. Vanation Tiree" 0a a eee 1 
4. Simulation Vanation Four 2° 2-- > a es ees. 77 
5. Simulation Varlationtcivés. .. ee eee ee. 78 
D. RESULTS FROM OPUS-8 APPROXIMATION ............. 78 
E. ANADLYSESOFWRESWIPIS oie. 2. ei eee eee 80 


vill 


[Pe NITAINSe VOM MROtAIMICINEe sercciys . <0). ss eee ee we es 85 


DEAT AISeSPOl BACKORUCK Sar. sen. se (Sa Sesesils hyo Wow iptayerd a Sele 0 9] 

3. Analyses on Operational Availability, Ao(t), ............ 100 

EE SAI SCS Ol EIN (1) cee 5 owas sans Sot sailsns saa erie ekeeeree ie eo 104 

VII. SUMMARY, CONCLUSIONS AND RECOMMENDATIONS ..... 106 
SO AEM eu scenes OPS Sa ckiae , Wake Yardage 106 
ECON CIEE SI@IN Sa. a. fy ck hs oe «Rn cise ets Wes De 107 
Nee MVE NIDAMIONS: 5c. 2. Wiens ee Sa oles eee Soe 108 


APPENDIX A. MODSIM-II PROGRAM FOR WSM SIMULATION 
A) ede cian 2 Pe AR ee re ruins eel eenigr ata, eter Bip eit See 110 

APPENDIX B. MATHLAB PROGRAM FOR ANALYTICAL WSM.... 146 

APPENDIX C. NUMERICAL RESULTS FROM ANALYTICAL VERSION 


Rema eNO Vigan ke ee rate ana iO) ues Hee wyemna dda abi 154 


Aes DNS RY NOUV Me rage es og eh S Goha os Mictapie pn acnaer sity Gel Grn 158 

APPENDIX E. OPUS-8 INPUT REQUIREMENTS AND RESULTS .... 164 
A. OPUS-8 INPUT FORMAT AND REQUIREMENTS ......... 164 
B. RESULTS USING OPUS8’S ENDURANCE OPTION TO 


Ae tN OAM Ee VV-OIVA oa) 2 a oo 9iiy @ SoS aatanaei's wera ee Boers ene < 171 


ACKNOWLEDGEMENTS 
I would like to express my gratitude to all the persons and institutions that have 
given me this unique opportunity to widening my horizon at the Naval Postgraduate 
School. I am particularly grateful to the Defence Material Organization of the 


Singapore Defence Ministry for this sponsorship. 


My sincere thanks to my advisor, Professor Emeritus M. B. Kline who has given 
me great insights in the execution of this thesis. I am equally indebted to the 
guidance of my co-advisor Professor A. W. McMasters and second reader, Professor 


M. P. Bailey. 


Last but not least, my heart-felt appreciation to my wife, Violet, who always 
supported and encouraged me especially during those times when the going was 


tough. 





I. INTRODUCTION 


A. MOTIVATION FOR THE THESIS 

The Ministry of Defence (MINDEF) of the Republic of Singapore recognized 
the need to optimize its investment in spares to maintain the desired level of 
readiness of its weapon systems and, therefore, has purchased under license a 
software package called "OPUS" [Ref. 1]. Installed in 1985 and updated to the 
present version, OPUS-8 has since been used to perform steady-state optimal 
allocation of spares for many types of weapon systems. This has led to a better 
understanding of the kind of effects a given allocation of spares have on Operational 
Availability (Ao), a widely accepted system-oriented steady-state Measure of 
Effectiveness (MOE) used for peacetime deployment. 

However, during a wartime period, utilization of a military capability becomes 
highly dynamic and this can result in inevitable interruptions and delays in the 
logistics support for the military capability. OPUS-8 does not address this problem 
directly because of its steady-state assumptions. Moreover, the operational planners 
and users are more comfortable with the use of mission-oriented MOEs to define the 
system’s readiness during such dynamic periods. 

For these reasons, both the Joint Logistics and Operations Analysis 


departments of MINDEF have specifically expressed the need to establish sound 


policies and approaches to relate spares and related logistics resources to a mission- 
oriented MOE under a dynamic wartime environment. 

Existing sustainability models such as Dyna-METRIC [Ref. 2] and the Aircraft 
Sustainability Model [Ref. 3] have been developed with the specific purpose of 
studying the effects of spares on the mission-readiness of an aircraft squadron as a 
detachment unit in a remote location where logistics support is limited. Although 
the basics of these models are readily accessible, the use of these models are 
restricted to US and NATO military agencies. Hence, the author was motivated to 
investigate the concepts and algorithms of these models with the aim of developing 


a specific Wartime Sustainability Model (WSM) for use in MINDEF. 


B. GOALS 

A goal of the thesis was to formulate relevant policies, analytical assumptions 
and rationales for the development of an analytical model for the WSM. 

Another goal was to develop a simulation version for the WSM to study the 
effects of policies such as cannibalization, repair prioritization and limited repair 
resources on the sustainability of a military capability during an anticipated wartime 
period. Such policies are extremely difficult to model analytically. 

Both models are verified with a numerical scenario and the results contribute 
towards a better understanding of the transient behavior of a military capability 


under a dynamic combat environment. 


Another objective of this thesis was to use the existing OPUS-8, a steady-state 
model, as an approximation to the proposed WSM under the special case of no 


repair capability. 


C. PREVIEW 

The basic layout of this thesis and the relationships between the chapters and 
the appendices are illustrated by Figure 1-1. 

Chapter I introduces the reader to the motivation and goals of this thesis. 
Chapter II establishes a scope or framework for the development of the WSM in a 
systematic and effective manner. Chapter III shows the development of the 
analytical model for the WSM under specific stochastic assumptions. Chapter IV 
describes the simulation version for the WSM under different situations using a 
fourth generation simulation language, MODSIM II [Ref. 4], to implement the 
simulation. The MODSIM Program listing of the simulation model is provided in 
Appendix A. Chapter V explores the possibility of using the current features of 
OPUS-8 to approximate the WSM. Chapter VI contains many numerical results and 
graphs based on an example scenario, and analyses of these results are provided. For 
efficient computation, the exact analytical expressions derived in Chapter III were 
coded into a computer program using PC-MATLAB ([Ref. 5] syntax, providing 
numerical solutions to the required numerical examples. This PC-MATLAB 
program listing is given in Appendix B. Examples of computer outputs from the 


analytical model, the simulation model and the OPUS-8 Approximation are also 


provided in Appendices C, D and E, respectively. Chapter VII summaries the thesis 


effort and provides conclusions and recommendations for further model development. 


Introduction 
Chapter I 










Framework For Wartime 
Sustainability Model (WSM) 
Chapter II 


Analytical Wartime 
Simulation Model 


(WSM) Chapter III 





Il. FRAMEWORK FOR A WARTIME SUSTAINABILITY MODEL (WSM) 


A. SUSTAINABILITY DEFINITION 

Sustaining a military capability in a wartime environment is a complex 
combination of many military aspects such as manpower, tactics and logistics support. 

For the purpose of this thesis, the term "sustainability" applies only to the 
prediction of the effects of logistics support policies on the ability of the military to 
Sustain its capability during an anticipated wartime period. For example, the military 
capability can be a squadron of aircraft deployed at a particular base supported by 
many logistics resources. 

Because of the short period of the time allocated for this thesis, the design of 
the WSM is restricted to being an assessment tool for analyzing the effects of a 
prespecified allocations of spares and repair resources on the mission readiness of 
the military capability. In other words, the WSM to be developed in this thesis is not 
an optimization tool capable of recommending further requisition of an optimal 
number of spares and repair resources when the initial allocation is determined to 
be inadequate. Further development of the WSM Is envisioned when the author 
returns to Singapore. 

Policies and Measures of Effectiveness (MOEs) which have influences on the 


design and development of the WSM are discussed in the following sections. 


B. POLICIES 

There are many problems involved in trying to establish sound policies to guide 
the determination of wartime logistics support. One of them is the lack of a common 
terminology for this area which leads to much confusion and misunderstanding 
between the operational planners and the logisticians who are responsible for 
implementing these policies. Another is that even when it is accepted that 
operational requirements must be defined first before logistics resources can be 
determined, these requirements are often not specific enough to be translated into 
reasonable logistics objectives. 

Research carried out by the Logistics Management Institute of U.S.A. 
emphasized that proper determination of war reserve spares requires clear and sound 
policies which must be endorsed by both the users and the logisticians 
[Ref. 6]. The same reasoning should apply in the context of MINDEF and 
the author feels that the following policy issues should be addressed as part of 


developing an acceptable WSM. 


1. War and Mobilization Plan (WMP). 
Each major military capability in the Singapore Armed Forces (SAF) has 
a War and Mobilization Plan formulated by the MINDEF strategic planners. In the 
context of the WSM, the WMP should address the following issues: 


a. the total wartime period which includes the pre-tension period, the surge 
period and post-tension period; 


b. the utilization of this capability during these periods; 


Cc. war attrition; and 


d. the minimum critical numbers for the various facets of the capability to 
sustain the wartime period. 


An example of what might be included in the WMP is the utilization 


profile shown in Figure 2-1 for a squadron of aircraft during an anticipated wartime 


period. 
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FIGURE 2-1. TYPICAL WARTIME UTILIZATION PROFILE 


It is assumed in this thesis that in the earlier stages of planning to acquire 
the assets needed to achieve the military capability, war attrition analysis has already 
been performed to ascertain the minimum critical mass of the assets required to 
Sustain an anticipated war. As a consequence, the effects of attrition will not 


included in the design of the WSM. 


2. Optimization of Investment In Spares 
Given the limited allocated funds for defense, the need to maximize a 
military capability for a given dollar amount cannot be overemphasized. Although 
not specifically addressed in this thesis, there is a need to buy the correct assortment 
of spares and repair resources for a particular asset investment to achieve the 
maximum MOEs with a limited amount of funds. The proposed WSM can be used 
in a limited way to answer this policy. But, when expanded into an optimization tool, 


the WSM will be capable of fully addressing this issue. 


3. The Inadequacy of A Minimum Buy Policy 
In the past, it has been common for the SAF to buy war-reserve spares 
based on a minimum-buy policy. However, such a minimum quantity was often 
bought based on the suppliers’ recommendations. Another method was to buy 
spares for each item to meet a specified service level based on the assumption of 
Poisson distribution for the demands for spares over the specified period of conflict. 
Both approaches disregard the need for optimization with respect to any MOEs. In 


fact, the minimum-buy policy prevents resources from being available for a desired 


MOE. This has been shown to be true by a study carried out under the supervision 
of the Joint Logistics Department of MINDEF which highlighted the urgency to 
minimize the excessive amount of unused spares in the SAF inventory system. 

Also, the minimum-buy policy fails to take into account intense fluctuations of 
demands during a wartime period and therefore is not able to accommodate the 
demands for spares and repair resources during a surge period. The proposed WSM 
Is a more effective and systematic approach for surge protection during a wartime 


period. 


4. Cannibalization Policy. 

During peacetime deployment it makes sense to discourage cannibalization 
of systems for spare parts since this practice can create havoc within any logistics 
management accounting system. However, in time of war where the chief objective 
is tO maximize the utilization of all available systems to accomplish a mission, 
studies have shown that cannibalization does improve operational availability (Ao). 
However one study recommended a policy to control cannibalization by setting an 
upper limit on the number of systems to be cannibalized [Ref. 6]. This study also 
showed that a correct choice of this limit on cannibalization can maximize certain 
MOEs. The effects of cannibalization are considered in the development of the 


WSM. 


5. Repair Prioritization Policies 

An appropriate choice of a repair prioritization policy is important when 
repair resources are limited, since a decision to repair one item also means a 
decision not to repair another when both are competing for the same repair 
resources. Repair priority policies such as First Come First Served (FCFS) and Least 
Availability Item First (LA/F) are commonly adopted. Of course, other policies may 
be more suitable under certain circumstances. 

With the above in mind, the properly defined WMP enables the 
operational planners to develop relevant Measures of Effectiveness (MOEs) which 
can be used as design criteria for the development of the WSM. Details of such 


MOEsS are discussed in section II-C. 


C. MEASURES OF EFFECTIVENESS (MOEs) 

The WMP provides a basis for determining the relevant Measures of 
Effectiveness (MOEs) which can serve as design goals for the development of the 
WSM. MOEs can be further classified into System-oriented MOEs and Mission- 
oriented MOEs and their relationships are depicted in Figure 2-2. The dependence 
of these MOEs on the performance characteristics of the military asset is also shown. 
However, the effects of performance attributes on a military asset will not be 


considered in the WSM since this is beyond the scope of this thesis. 


10 





Mission-oriented 
Measures Of Effectiveness = Examples 
ENMCS, P[LENMCS<NMCS] 
ft fh Kill Probability, 
Attrition rate. 


MISSION EFFECTIVENESS ANALYSIS 


=> 
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System-oriented 
Measures Of Effectiveness = Examples 
Ao(t) (Wartime), 
f f Ao (Peacetime), 
MTTR, MLDT, MTBF. 


SYSTEM EFFECTIVENESS ANALYSIS 


= 
= 


Performance-oriented 
Measures > Examples 
Altitude, Range, 
f f Speed, Radar Cross 
Section. 


PERFORMANCE ANALYSIS 





FIGURE 2-2 : RELATIONSHIPS BETWEEN THE MOEs FOR 
A MILITARY CAPABILITY 


1. Mission-oriented MOEs (MMOEs) 
The proposed MMOE for the WSM is NMCS where 
NMCS = the maximum allowable number of Not-Mission-Capable Systems that 


can be tolerated without reducing the military capability during an 
anticipated wartime period. 


11 


NMCS a mission objective which should be specified by the operational planners and 
its value should set to at least a 95% upper confidence limit. It is important to 
emphasize that the proposed WSM presumes that the user had undertaken 
comprehensive studies (analyses or wargaming) to obtain the NMCS. Therefore, this 
objective must be in the WMP for the desired military capability. The WSM 
computes the Expected number of Not-Mission-Capable Systems (ENMCS) based on 
a given allocation of spares and repair resources. 

The military capability is not severely downgraded as long as ENMCS is less 
than NMCS. Reference 6 preferred this MMOE to avoid possible statistical 
confusion associated with confidence-level oriented objectives. 

Another MMOE candidate is the Confidence Level of having not more 
than a specified number of failed systems throughout the anticipated wartime period. 


It can be formulated as Probability [ENMCS <= NMCS]. 


1. System-oriented MOEs (SMOESs) 

Many steady-state spares models including OPUS-8 use Operational 
Availability (Ao) asa key SMOE. However, in a dynamic environment, Ao becomes 
Ao(t) which is time-dependent. Ao(t) is computed based on the ratio of the expected 
number of Mission-Capable Systems (MCS) and the total number of deployed 
systems. It is shown in Chapter III that ENMCS and Ao(t) have a direct 


relationship. 
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The main focus of the thesis is on the effects of spares on both Ao(t) and 
ENMCS._ In Chapter III, other two SMOEs, the expected number of demands for 
spares and the expected number of backorders for each Line Replaceable Unit 
(LRU), are shown to be basic building blocks of the WSM for achieving desired 


Ao(t) and ENMCS levels. 


D. SYSTEM STRUCTURE AND LOGISTICS SUPPORT 

In order to be operationally ready, the deployed systems have to be supported 
by many logistics resources. A typical organization for logistics support of a major 
system having multiple indenture-levels is depicted in Figure 2-3. The base repair 
and supply facilities are shaded to emphasize the proximity of these locations to the 
deployment sites. 

When a Line Replaceable Unit (LRU) of the system breaks down at the base, 
the system will enter a base repair network which has support capabilities (men, 
machines, spares and supplies) to remove and replace any system’s faulty LRU with 
a good one. Depending on the nature of the failure, the failed LRU is sent either 
to the base repair facility, the intermediate repair facility or the depot repair facility 
for the appropriate repairs. 

Demands for spares of an LRU are intense during a wartime period and 
therefore the system requires a sufficient number of spares of each LRU so that the 


system will never be unavailable due to lack of spares. 
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DEPOT LRU 
AND SRU 


REPAIR AND SUPPLY 


INTERMEDIATE LRU 
AND SRU 


REPAIR AND SUPPLY 


WARTIME MOBILIZATION 
PLAN AND POLICY 





FIGURE 2-3 : A TYPICAL LOGISTICS SUPPORT ORGANIZATION 
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II. AN ANALYTICAL WARTIME SUSTAINABILITY MODEL (WSM) 

This chapter first describes the assumptions required by the WSM analytical 
model to obtain the relevant exact analytical expressions for the SMOEs and 
MMOEs. In the development of the formulae, proofs for well-known theorems will 


only be referenced to the appropriate sources. 


A. ANALYTICAL MODEL ASSUMPTIONS 
The following assumptions are made to enable the development of an analytical 
model for the WSM: 
a. Failures of the Line Replaceable Units (LRUs) are generated at the base 
level. These failures are repairable at the base with a certain probability, 
otherwise they will be shipped to the depot for repair. Also, it is assumed 


that there will be no preventive maintenance in a wartime environment. 


b. All faulty LRUs shipped to the depot can be repaired (i.e., there are no 
condemnations or attritions at the depot). 


c. The distribution of the interarrival times of any LRU failure follows either a 
nonhomogeneous Poisson distribution or a nonhomogeneous compound 
Poisson distribution. 


d. The failure of one type of LRU is statistically independent of those which 
occur for any other type of LRUs. 


e. The repair times and transportation times are statistically independent of one 
another. These times can be stationary or non-stationary. 


f. Repair resources are assumed to be unlimited. 
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B. FAILURE ARRIVAL RATES 
Let's look at a particular LRU of the system. Since the system has a 
varying utilization rate during an anticipated wartime period, the arrival rate of the 


LRU failure, D(t), at an instant in time has the following formula: 


D@) = ff. Gg do AsO); (3.1) 

where 

f = the basic item failure rate of the LRU. Here, interarrival times are also 
assumed to follow an exponential distribution (i.e., it 1s constant). 

q = quantity of the LRU found in each system. 

u(t) = the utilization rate per day of all the deployed systems at time ¢. Its value 


varies with time based on the WPM as specified by the operational 
planners. For example, if the systems are required to operate for 12 hours 
a day at time ¢, then u(t) has a value of 0.5. 


ns (t) the number of available systems at time ¢ This also changes with time 


since some systems may become unavailable due to the lack of a spare. 


C. DEMANDS FOR SPARES 

In this section we present the derivations of the exact analytical expressions for 
the distribution of the number of failed units of a particular LRU being held up in 
repair at a given time ¢. This number of failed units generated demands for spares 
for the same LRU which were hopefully satisfied by the inventory in stock. Systems 


become Not-Mission-Capable when there is inadequate inventory. 
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1. Nonhomogeneous Poisson Assumption 


a. The Distribution of LRU Failure Amrtvals 
Let {N(t), t 2 0} be the counting process of a particular LRU’s 
failure arrivals which follows a nonhomogeneous Poisson process. Also assume that 
there are no arrivals before time t=0 (1.e., the arrival process starts empty; that is, 


N(0)=0). The expected value, denoted by M(t), for the process N(t) is 


M(t) = E[N(t)] = l ' D(s) ds . C2) 


b. The Distribution of the Repair Process 
Corresponding to an LAU failure arrival that occurs at time s, let the 
repair time random variable denoted by Y, have a distribution function G(s, s+y) 
which is dependent on s. The actual repair times y are also assumed to be 
independent. 
c. The Distribution of Failed Units of an LRU Sull in Repair 
We wish next to determine the number of arrivals still in repair, 
denoted by Z(t), for the case of unlimited repair resources. 
Suppose the Poisson arrivals have a constant arrival rate D and 
undergo repair with repair times independent and identically distributed (iid) from 
a Stationary distribution F. Ross [Ref. 7], Example 4b, p.237, showed that Z(t) also 


has a Poisson distribution with expected value of 
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E{Z(t)] = A(t) = {, D[{1 - F(s)] ds . (3.3) 


Takacs [Ref. 8], p.160, further showed that ast => ~, the limiting distribution of 


the number of arrivals still in repair follows a Poisson distribution with mean 


lim A(t) = DEY). (3.4) 


pee 


so long as the limit is finite (i.e., 1/D = E[Y)). 
This is widely known as Palm’s Theorem. 

Palm’s Theorem was extended further by Hillestad and Carillo 
[Ref. 9] to handle a particular LRU having a nonhomogeneous Poisson failure and 
undergoing a repair process. The assumption that the failure arrivals and repair 
times are independent of one another is needed. Under these conditions, Z(t) also 


has a nonhomogeneous Poisson distribution with the expected value, denoted by A(t), 


as 
At) = [ [1 - Gs)IDG) ds (3.5) 
where 

5 = time when the repair was initiated, and 


time of interest. 


f 
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Finally, the distribution of Z(t) is given by 


mia) = ki = e 0 or (3.6) 

Instead of P{/Z(t) =k], the notation P/k ; A(t)] will be used from now 
on where A(t) represents the expected number of failed units of @ particular LRU 
still in repair. A(t) can also be interpreted as the expected number of outstanding 
demands for spares for that same LRU. If there are enough spares in stock, then 
these demands do not cause a system to be Not-Mission-Capable. 

Since Z(t) is a Poisson random variable at time ¢ its variance, 
Var[Z(t)], is equal to E/Z(t)](i.e., A(t)). The ease of computation resulted from this 


closed-form expression will be demonstrated in Chapter VI. 


2. Nonhomogeneous Compound Poisson Assumption 
a The Distribution of LRU Failure Amvals 
Intense failure arrival rates of a particular LRU during a wartime 
Scenario can be expected. Also it is reasonable to expect an arrival to consist of a 
batch of units of the LRU instead of just one. 
W,, is defined to be the batch size at the nth arrival of the Poisson 
process N(t). It is assumed that the W, ,n = 0, J,... are ud random variables having 


a common compounding distribution { C; = P[W=j],j = 90, 1,...}, with its expected 


19 


value defined as E/W]. Assume further that the family { W, } is independent of the 
arrival process. 

With the above conditions, we define {X(t), t = 0} as the resulting 
counting process of the number of failure arrivals for a particular LRU, where 


N(t) 


Xt) = YW, (3.7) 


Ross [Ref. 10], p.49, has shown that X(t) is a compound Poisson process having the 


following expressions for its mean and variance, 


E[X(t)] = M(t)E[W] (3.8) 
and 
Var[X(t)] = M(t)E[W?) . (3.9) 


b. The Distribution of Failed Units of an LRU Still in Repair 
Again, we are interested in Z(t), the number of failed units of a 
particular LRU still in repair. If the failed units of the LRU in each batch have 
independent repair times from a distribution G (this can be treated as a special case 
of an M/G/o queue with batch arrivals and general service), then the resulting Z(t) 


also follows a compound Poisson process with 


- CY) e-X [A()p 
PIZ() =k) = SS - = (form Oe: | om) 
[rene J: 
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where {C,%} is the j-fold convolution of the compounding distribution {C,} with 
itself. The proof for Equation 3.10 is given by Feeney and Sherbrooke 
[Ref. 11]. Therefore Z(t), the number of failed units of LRU still in repair, 
is nonhomogeneous Poisson with compounding distribution { C,, j = 0, 1, 2, ... }. 
Also Z(t) may not have a closed form expression due to the convoluted term, {C,”}, 
in Equation 3.10. However it is shown in Feller [Ref. 12], p.291, that if W 
is assumed to follow a logarithmic probability function with probability density 
function (PDF) of 

P{[W=j] = -j[1 - a] / [jlog(a)], for 0 < a < 1; j=1,2.... (S211) 
and its the expected number given as 

E[W] =- [1-a] /alog(a) (Az) 
where 


I/a = E[W*]/E(W] = VMR, 


and VMR is known as the Variance-to-Mean Ratio of the distribution with a constant 
value always greater than 1, then the distribution of X(t) follows a negative binomial 
distribution, having parameters { a, R(t) } where 


R(t) = - M(t) / In(a) . (Sal) 
Thus, according to Feller [Ref. 12], Z(t) will also have negative binomial distribution 
with parameters 


{ a, R(t) } where 


ZA 


R(t) = - A(t) / In(a) . (3.14) 
Also, its PDF is 


RUG) ice 
k 


P[Z(t)=k] = ( 1) RO (1 - a)* fork=0,1,... (3.15) 


D. MULTIPLE-INDENTURE PIPELINES 

In the previous section, exact analytical expressions for the number of demands 
for spares for each LRU have been derived. Therefore, the number of units of the 
LRU in repair at the repair base, repair depot as well as in transit between the base 
and the depot can be computed using Equations 3.6 and 3.15. From here onwards, 
the term Pipeline is used to replace the phrase "expected number of failed units in 
repair and/or in transit of a particular LRU". For example, the expected units in 
repair at the repair base will be referred to as the Base Pipeline. Pipeline can also 
be interpreted as the expected number of demands for spares for that LRU. 

This section further develops the expressions for the pipelines to accommodate 
a system with more than one level of item breakdowns (i.e., multiple-indenture 


system). 


ee 


1. A Case of Two-Indenture System 

A system can be broken down into many indenture levels depending on 
the complexity of the system design. For the purposes of this thesis, two levels of 
item breakdowns are analyzed (i1.e., a two-indenture system). 

The terminology of Line Replaceable Units (LRUs) is used for the major 
items whose failures cause the system to be down if no spare LRU is available. Shop 
Replaceable Units (SRUs) are second-level items which make up the LRUs. The 
failure of an SRU creates a "hole" in an LRU which may cause a failure of the LRU. 
Therefore, a failed SRU may indirectly lead to a system being Not-Mission-Capable. 

Muckstadt [Ref. 13] has shown that for steady-state cases, Little’s 
Formula, a classical queuing result, can be adapted to analyze each repair pipeline 
given that an LRU has to wait for SRUs for its repair. The approach used in the 
WSM 1s similar except that steady-state assumptions can now be relaxed. 

Let Q.(t) represents the quantity of LRU. waiting for SRUs at time ¢ and 
let EQ,(t) be its expected value. This quantity is also equivalent to the number of 
LRU. being held up at the repair facility due to the lack of SRUs. Because of this, 
additional demands for spares are generated. For ease of explanation, the following 
analysis considers only a particular base but the results are applicable to any base or 
depot. Based on the nonhomogeneous Poisson assumptions made in subsection III- 


Cl, let A°.’(t) be the Base Repair Pipeline for LRU;. Then according to Muckstadt 


[Ref. 13], the Base Repair Pipeline of LRU;, when taking into account the pipeline 


due to lack of SRUs, is revised to become 


No) = 8) + FO). Cm 


The analysis carried out in this subsection is only valid for LRUs which 
undergo repair. When there 1s no repair allowed for the LRU, then there will not 
be any shortages due to waiting for SRUs. The following analyses also assume the 
independence between the LRU and its SRUs’ demand distributions even though it 


is more realistic that the failure of an LRU 1s due to its SRUs’ failures. 


a. Canzubalration Policy 
We now look at how cannibalization can affect the computation of 
EQ,(t). First, shortages of SRUs can be consolidated into the smallest possible 
number of LRU;. This is accomplished by using a serviceable SRU of a failed LRU 
to repair another LRU which requires the serviceable SRU. Let P!(n,t) be the 
probability that the shortages of the SRU; are less than or equal to a quantity n at 


time t. Then Pi(n,t) is 


St)en 
Pint) = > [ P( number of SRU, failures=m at time ¢ )] (3.17) 
m=0 


where S/(t) represents the inventory level of SRU, at the repair site at time ¢. 
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The exact form of the probability P(.) in Equation 3.17 depends on whether the 
demands for SRU; are Poisson or Compound Poisson. Now let P;(n,t) be the 
probability that the number of LRU, waiting for its SRUs to be repaired is less than 
or equal to m Failure independence among the SRUs and among the LRUs are 
assumed. Assume that the LRU, is made up of different types of SRUs connected 
in series, then P,(n,t) can be obtained by the multiplication of all its SRUs’ P}(n,t). 


This is an application of reliability theory for a serial model. Therefore P.(n,t) is 


Git) = IT Pi(n,t) , (3.18) 
jel 


assuming that there are a total of J different SRUs in LRU. . 
Then the expected number of unavailable LRU, due to shortages of SRUs can be 
computed from the conditional expectation of the above distribution 
K 
EQ) =X [i - Paw), (3.19) 
neQ 
where K is the total number of SRUs available which is equivalent to 


K ns(t) J + S\(t), and 


S(t) = the inventory level of SRU, at time ¢. 


From its distribution function given in Equation 3.18, the probability density function 
(PDF) of the number of LRU, that have to wait for SRUs at time t can be computed 


as 


PDF (nt) = P.(n,t) - P.(n-1,t). (3.20) 
In other words, PDF,(n,t) is the probability of having exactly n units of LRU,, each 


waiting for at least one of its SRUs to be repaired. 


b. No Canmnibalization 
Reference 13 also showed that the computation of EQ,(t) when no 

cannibalization of SRUs is allowed can be developed as follows. The probability that 
LRU, is short of its SRU,, given that there are B(t) Backorders of SRU; at time ¢ can 
be written as 
B(t)/[ns(t) + S'(t)] if we assume that there is only one of each type of SRU; in LRU, 
The formulation is true only if: 

(a) it is finite, 

(b) failure times are exponential , and 

(c) there is no information about (0, f}. 
The probability that LRU; has a shortage of SRU; is computed by summing over all 
possible B(t) = k - ns(t)+S'(t) multiplied by the probability that they occur, which is 


SK = OS Ol py , Ay) = FBO (3.21) 
kens(t)+S(t)+1 ns(t) +5 1(t) ns(t) +S 1(t) 
where 
A(t) = expected number of demands for SRU; at time ¢, 


and because the numerator 
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lk - ns(t) SHIP s NO) 
k=ns(t)+S(t)+1 
is equivalent to EB'(t), the expected number of Backorders of SRU, at time t. This 
equivalence will become evident later in section III-Fla. Therefore, the probability 
that LRU; does not have any shortage of SRU; is the complement of the distribution 
shown in Equation 3.21. 

Assuming that all failures are independent and all the SRUs are 
connected in series to form LRU,, the probability that LRU, does not have any 
shortage of all its SRUs is just 
TI 1 - _ EBX!) = P{ LRU, is Mission-Capable ] . (3.22) 
jel ns(t)+S\(t) 

Again, the complement of the distribution in Equation 3.22 is equivalent to the 
probability that LRU, has shortage of all its SRUs. If we assume that there is only 
one unit of each SRU; found in LRU; then the expected number of LRU; having 
shortages of SRUs is given by 


, : EBi(t) 
ie * = +§ J ® = = —SEEEq SSE SSS e (3.23) 
20) = lenses - TG - 


Zi 


E. MULTIPLE-ECHELON PIPELINES 
In this section we examine the logistics support for the entire deployment of 
systems. A Total Pipeline is obtained by combining the Pipelines at the repair 


facilities and in-transit between these facilities. 


1. A Case of Two-Echelon Logistics Support 

In this analysis, it is assumed that there is only one base repair facility and 

one depot repair facility. Let a failed LRU be repairable at the depot level with 

fixed probability commonly known as Not-Repairable-This-Station (NRTS) and at the 

base level with probability (J - NR7TS), independent of where other failed LRUs are 

repaired. Then, the single nonhomogeneous Poisson arrival stream with mean value 

function M(t) becomes two independent nonhomogeneous Poisson streams [Ref. 9] 

with mean value functions expressed as 
NRTS M(t) for the depot, , 


(3.24) 
(1-NRTS) M(t) for the base. J 


When the Repair Pipeline at each repair facility is assumed to follow a 
nonhomogeneous Poisson process, it will be denoted by 
A°(t) = Base Repair Pipeline 


= expected quantity of LRUs being repaired at the base (including the LRUs 
waiting for SRUs) at time f, 
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A‘(t) = Depot Repair Pipeline 
= expected quantity of LRUs in repair at the depot (including the LRUs 
waiting for SRUs) at time ¢ 


and the formulae for A(t) can be obtained from Equations 3.16 and 3.24. 
a. Total Repair Pipeline 
The following analysis reveals that there are two cases for which the 


Total Repair Pipeline can take on exact analytical results. 


(1) Nonhomogeneous Poisson Pipelines 

When the Base Repair Pipeline and the Depot Repair Pipeline are 
nonhomogeneous Poisson processes and are shown to be independent (allowing 
Equation 3.24), then the Total Repair Pipeline is also a nonhomogeneous Poisson 
properties with its mean value equal to the sum of the values of the Base Repair 
Pipeline and the Depot Repair Pipeline. The proof for this has been provided in Ross 
[Ref. 7]. Therefore, the mean of the Total Repair Pipeline, A,(t), for a particular LRU 
at any time ¢ is just the summation 


pay A(t) + AX). (3.25) 


Figure 3-1 illustrates Equation 3.25 for the data from the numerical example in 


Chapter VI. 
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FIGURE 3-1: NONHOMOGENEOUS POISSON PIPELINES. 


(2) Nonhomogeneous Compound Poisson Pipelines 
In the case of nonhomogeneous Compound Poisson pipelines 
where the Base Repair Pipeline and the Depot Repair Pipeline have negative binomial 
distributions with the same VMR, then the resulting Total Repair Pipeline also has a 


negative binomial distribution with parameters { a, R(t) } where 
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a = 1/VMR } 
(3.26) 
R(t) = - [AX(t) + AX(t)] / log(a) 3} 


Sherbrooke [Ref. 14] also obtained similar results for demands with 
geometrical compounding distributions. However, if the Base Repair Pipeline has a 
different VMR from the Depot Repair Pipeline even if both have logarithmic Poisson 
distributions, then the total quantity would no longer be negative binomial (i.e., the 
expression given in Equation 3.26 no longer applies and there is no closed form 


expression). 


b. Total Pipeline 

It is evident that more systems will be become unavailable if there are 
more LRUs being repaired at the repair facilities or in transit between these 
facilities. Therefore, the total expected number of failed units for an LRU still in 
repair (Total Pipeline) has to be ascertained in order to compute the MMOESs of the 
WSM. 

In a real situation, the Total Pipeline is a complex convolution of both 
the Base Repair Pipeline, Depot Repair Pipeline and other logistics Pipelines (usually 
transportation). Because of this convolution of many pipelines, the outcome for the 
Total Pipeline does not have an exact closed-form expression. This subsection derives 
the expression for the Total Pipeline when all the Pipelines are assumed to have 


nonhomogeneous Poisson properties. However, the same approach is equally 


ol 


applicable to Pipelines with nonhomogeneous compound Poisson properties under 
specific conditions explained earlier. 

The following time-dependent notations will be used for an arbitrary 
LRU, and, to simplify these notations, the subscript / used for LRU type will not be 
included. 
A” (t) = Base Repair Pipeline generated from base k; 
A“\(t) = Depot Repair Pipeline generated from base k; 


A(t) = Forward Pipeline generated from base k, 
In transit from base k to depot; 


A'\(t) = Return Pipeline generated from base k, 
In transit from depot & to base; 


S“(t) = Supply level at the depot; 
S,(t) = Supply level at the base k; 
T™, = Transportation time from base k to depot (Forward Time); 


Tl, = Transportation time from depot to base k (Return Time). 


Only one repair base and one repair depot will be considered. Figure 

3-2 illustrates the case. The transportation Pipelines are included in the analysis. 
When no inventory of the LRU is held at the depot, a defective LRU sent 
to the depot must be repaired at the depot before it can be returned to the base for 
stocking. For convenience, let A(t) = A‘,(t) + A%,(t) and refer to this term as the 


Depot Forward Pipeline. 


Sue 


In this case, the expected Total Pipeline is given by 


A(t) = A(t) + AN (t) + AY (t) . (3.27) 
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FIGURE 3-2 : A SIMPLE TWO-ECHELON MAINTENANCE SYSTEM. 
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On the other hand, when inventory is available at the depot, the base 
can now depend on the depot inventory to reduce the waiting time for LRUs that 
were sent to the depot. Instead of A'%,(t), the base will only feel the shortages due 
to the unfilled orders placed at the depot. Let this expected number of Backorders 
at time t be EB S(t), A‘%,(t) ]. Then the depot portion of the Total Pipeline at base 
k, denoted by A",(t), is the sum of these Backorders and the Return Pipeline. Thus, 


Mito) =~ EBS Se ae) |p oe (3.28) 


Therefore, the expected Total Pipeline at base k is 


Agj= NMO+Awdw . (3729) 


However, time dependency does present difficulty in determining the 
A’ (t) since this quantity is conditioned on the stock available at the depot, 
completion of the depot repair and 7",. A better way of looking at this quantity is 
to study the quantity of demands during one 7", since all these will be unfilled (.e., 


the shortages that exist started at time [t - 7",]. Thus 


Nn) = { D0) ds + EBS%¢ - 71), Nee - 7) (3.30) 


I 
t -T, 


where 


D°\(s) = The depot failure arrival rate of a particular LRU generated from base k 
Starting at time s. 
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For multiple bases served by only one depot, the approach is the same except that 


the depot supply level, S°(t) now is used by all the bases. Therefore 


K 
Nr) = DNA) (3.31) 


and the expected number of shortages in the Depot Forward Pipeline is given by 
EB (S(t), A‘°(t)). 

The next step is to allocate these shortages to all the bases depending 
on the choice of the repair policy. One alternative suggested by Sherbrooke 
[Ref. 15] for steady-state cases is to allocate based on the relative quantity 


of demands placed at time 4. That ts, 


d 
D,(t) : 
EB) (t) = —~—— EB4S%s), A%0) . (3.32) 


y Doo) 


k=1 


However, sudden changes in demand at the bases may render this technique rather 
unstable. A more stable method suggested by Hillestad [Ref. 16] is to employ time- 
averaged demands at the bases. To do this, let DA‘%,(t,6) be the average demands 


in interval /t-5,t] where é is a small interval of time. Therefore 


DAY(t,6) = f Day at (3.33) 


> 


Hence the new allocation rule is: 


d 
DA p(t,6) 


EBA) - EBS (1), A%(t)) . (3.34) 


K 
d 
2 DA, (t,6) 


The influence of past demands can be fine-tuned by changing the value of 6. 
Reference 16 also observed that a larger 6 value would increase the influence of past 
demands and a good value would be near the average time the depot takes to 
reallocate the shortages. Therefore using 6 = ¢ will keep the allocated number of 
Backorders below the expected Total Pipeline from time {0 to t}. The approach is 


also widely used by Sherbrooke [Ref. 15]. 


F. THE SMOEs 

It has been shown in section E of this chapter that the expected Total Pipeline 
experienced by a base is dependent on the various pipelines of the logistics flow 
depicted in Figure 3-2 and the available inventory levels at the base and the depot. 
Combining these elements allows us to assess various measures describing the 
availability of LRUs required to keep the systems operationally ready. 

Although these SMOEs do not directly describe the system mission readiness, 
they do form a basis for computing MMOEs and also are useful in identifying LRUs 
that cause the system to be unavailable for missions. This subsection describes the 


analytical approach to derive these SMOEs. 
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1. EB,(t), Expected Number of Backorders for LRU; 
The computation of the expected number of backorders (also referred to 
as "Backorders" in this thesis) for a particular LRU,, denoted by EB.(t), depends on 
Whether the pipelines follow the nonhomogeneous Poisson process or the 


nonhomogeneous Compound Poisson process. 


a. EB;(t) for Nonhomogeneous Poisson Pipelines 
Suppose the initial inventory level for a particular LRU, is S(t). The 
purpose of S.(t) is to protect the systems against the Total Pipeline demands for 
spares of the LRU. When that demands exceeds S(t), the LRU, is said to be in a 
"packordered” state. Since systems require this LRU to be Mission-Capable, any 
Backorder situation affects the MMOEs. 
Suppose that the expected Total Pipeline A,(t) for LRU, follows the 


nonhomogeneous Poisson distribution, the EB.(t) at time ¢ is given by: 


EB(t)= 0 [k - S(t) Pik; AQ) 


k=Si(t)+1 
Si(t) 
= (A) -S()] + 30 1S) - Al Pik 5 AC) . (3.35) 


The latter expression is easier to program on a computer. 


o/ 


To complete the analysis, the variance of the Backorders is given by 


VB(t)= S[k - S()P PIk ; A@) - [EB()) 


k=S;(t)+1 
S;(t) 

= A(t) + [A(t) - S(OP - (EB)? - S[k - SCOP Plk s ACI. 26) 
k=O 


Again, the latter expression 1s easier to program on a computer. 


b. EB;(t) for Nonhomogeneous Compound Poisson Pipelines 
The exact expression of EB(t) for nonhomogeneous Compound 
Poisson pipelines can be carried out in a similar manner as shown in the previous 


paragraph and has been done by Sherbrooke [Ref. 15]. 


2. EB(t), the Total Backorders of the Deployed Systems 
The Total Backorders, denoted EB(t), is an SMOE which measures the 
expected total number of all Backorders experienced by all the deployed systems of 
the military capability. If there are J types of LRUs then EB(t) at any time t can be 
computed as 
I 
EB(t) = > EB(t) . (3239) 
iel 
With a policy of no cannibalization, EB(t) is a parameter used in the 
computation of the most conservative value for Ao(t). The exact expression for Ao(t) 


is Shown below in Equation 3.41. 
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3. Operational Availability (Ao(t)) 

Suppose that there are a total of / types of LRUs for each system and the 
total number of deployed systems are assumed to be identical both in hardware and 
functions. In the computation of Ao(t), all of the LRUs are considered essential for 
the operation of the system and therefore their failures cause the system to fail. It 
iS important to note here that the computation to be described later can be extended 
to handle different types of systems. Similar to the OPUS-8 algorithm, the overall 
Ao(t) for the whole deployment can be obtained by taking the weighted average of 
all of the individual systems’ Ao(t). 

a. Ao(t) without Cannibalization 

Here the policy of not consolidating shortages of LRUs among the 
systems is adopted. The probability that a single military capability consisting of a 
total number of ms(t) systems at time 4 is short of LRU. , given that there are b(t) 
Backorders of LRU, at an arbitrary base at the same time , can be written as 
B(t)/ns(t) under the same conditions used to derive Equation 3.21. Assume for now 
that there is only one of each type of LRU in each system, then the probability that 
these ns(t) systems has a shortage of LRU, is computed by summing over all possible 
b(t) = k - S(t) values multiplied by the probability that they occur, which is 
— [k ~ S(t)] 


KS(ye1 S(t) 


Pik ; A(t} = ae | (3.38) 
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Assuming that the LAU taihires are independent and noting that the probability of 
no shortage of each is the complement of the probability given by Equation 3.38, 
then using the reliability theory of items-in-series, the probability that the whole 
deployment is Mission-Capable is given by 


E'B(t) 


rus(f) 





l 
Ao(f) = pa 1 - 
t~| 


As mentioned in the previous subsection, there is a relationship between Ao(rt) and 


EB (rt) which can be shown by expanding Equation 3.39 as follows. 











I e( l-1 1 . 
@ige= fb - (9 ee eS ee (3.40) 
~~ vet yeret MES(f) — ras(T) 


where 


EB(t) = Backorders of LRU, at time +. 


When individual EB (r) are very small as compared to ms(t), then the second and 
higher terms of Equation 3.40 become negligible. Therefore, Equation 3.40 can be 


reduced to the approximate expression 


Aoi) = [1 - ane) (3.41) 
rg(t) 





lf there are several units of LRU, , denoted by g, . in each system, Hillestad [Ret.16] 


showed that Equation 3.39 can to be moditied to 


40) 


ns(t)q, - l 


1 ns(t)q; ( 


qd. 
A = ae oe , rae 
0 TE ag P80 3.42) 
q, 


where PB,(y) is the probability that LRU, has y shortages at time ¢ with the following 
conditions 


Sj(t) 


2, Pik; A] ify =0 


PBi(y) = 


Py + Sit); AO) ify > 0. 


b. Ao(t) with Cannibalization 

Again, suppose that there are J types of LRUs in each system and 
ns(t) is the number of identical systems. Under the cannibalization policy, shortages 
of LRUs can be consolidated into the smallest possible number of systems. 
Cannibalization is often practiced in military establishments to minimize the number 
of Not-Mission-Capable Systems (NMCS) in time of war. The rationale for doing 
this is that when a system becomes unavailable, it does not matter if there is only one 
missing LRU or more. Therefore, the good LRUs in a failed system can be used as 
spares to keep other systems Mission-Capable. 

Suppose that there is only one of each type of LRU in the system and 


let P'(j) be the probability that the number of shortages of the LRU, are less than or 
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equal to j. Also, assuming that the A(t), the Total Pipeline of the LRU., is from a 
nonhomogeneous Poisson process, then 
5 (9 +7) 


PEG) = ge A) (3:48) 


h =) 


Now let P(y) be the probability that the number of NMCS is less than 
or equal to j. Here j is the maximum number of NMCS which can be tolerated by 
the users. With the cannibalization policy, we can concentrate all the shortages of 
LRUs to the minimum number of NMCS (e., there are j NMCS if all its LRUs also 
have jy shortages). Assuming that each system is made up of / types of LRUs 
connected in series, then P(y) can be obtained by the multiplication of all its LRUs’ 


P'(j). When there is only one of each LRU in each system, P(j) 1s 


PU) = I Pj). (3.44) 


In the case where there are q; units of LRU; in each system, these 


umits can be used as spares for cannibalization if the system becomes NMC, then 
. 

EG) =e FP Gy) - (3.45) 
iol 

From the theory of stochastic processes, the expected value of a nonnegative discrete 


distribution can be computed from the sum of the values of the complementary CDF 


of that distribution. Since the upper limit of the summation is the number of 
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systems, s(t), being deployed at time ¢, then the expected number of NMCS at time 
t, ENMCS.(t), is 


ENMCS,(t) = > [1 - Py] . (640) 
jx 


Then Ao(t) with Cannibalization, denoted by Ao,(t), 1s 


[ns(t) - ENMCS.(t)] 
ns(t) 


Ao{t) = (3.47) 


G. The MMOEs 

The SMOE, Ao(t), computed in the previous subsection can also be used to 
derive the MMOEs which are measures more readily understood by the operators of 
the system. Although ENMCS is of main interest, the formulation for P/ENMCS < 


NMCS] is also given here. ‘The definitions of both MMOEs were given in section 


II-C1. 
1. ENMCS 
Without Cannibalization, the expected number of NMCS at time f¢ 1s 
ENMCS(t) = ns(t) - [ns(t) Ao(t)], (3.48) 


Depending on whether there is one or more of each LRU, the expression of Ao(t) 


can be obtained either from Equations 3.39 or 3.42. 
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In the case of Cannibalization, the expression for ENMCS(t) is exactly 


the same as Equation 3.46. 


2. P(ENMCS < NMCS) 

The analyses described in subsection III-F3 are useful in the formulation 
of P(/ENMCS < NMCS) with or without cannibalization. There j represents the 
maximum number of NMCS. Therefore, Equations 3.44 or 3.45 can also be used as 
P(ENMCS < NMCS). Since this MMOE is also time dependent, there is a value of 


P(ENMCS < NMCS) for each time interval of the wartime period. 


H. LIMITATIONS OF THE ANALYTICAL APPROACH 

The most obvious limitation of the analytical model for the WSM 1s the need 
to make certain simplifying assumptions so that closed-form expressions can be 
obtained. 

However, in a realistic wartime period the demands for spares are intense and 
highly unpredictable, and there will always be queuing at repair facilities due to the 
limited availability of spares and repair resources. Also, in the course of deployment 
management will try to take actions which it thinks may optimize the MOEs. 

Therefore, a Monte Carlo simulation model is needed to analyze the effects of 
more realistic scenarios and to verify whether the actions taken intuitively by 


management actually improve the ability of the logistics support system to meet the 


operational demands. The simulation approach for the WSM is developed In 


Chapter IV. 
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IV. A SIMULATION VERSION OF THE WSM 


A. REASON FOR THE SIMULATION APPROACH 

As stated in section H of Chapter III, the limitations of the analytical approach 
for the WSM warrants the need to perform Monte Carlo sampling to simulate more 
realistic scenarios so that the effects on the MOEs of the WSM can be studied. This 
chapter describes the design and the implementation of a simulation version of the 


WSM. 


B. MODEL DESIGN 

The simulation model design for the WSM has a structure similar to the 
analytical approach as shown in Figure 3-2. The design takes the form of a multi- 
server, multi-job-class queuing simulation with time-varying demand (arrival) rates 
for each LRU of the system. Unlike most queuing simulations which are designed 
to study the steady-state behavior of a system with constant parameters, the WSM 
simulation studies the behavior of a group of systems deployed as one military 
capability and its logistics support over a finite time period during which the demand 
rates vary dynamically. 

The flow of the various processes of the model are discussed in the following 
Sections. Issues concerning the need to replicate runs and their statistical effects on 


the SMOEs and MMOEs, sampling interarrival times, handling of random number 
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Streams, repair selection priority rules, and base and depot repair policies are also 
addressed. The main processes to be modelled for the WSM Simulation are 


provided in Table 4.1. 


TABLE 4.1 : DESIGN OF FLOW PROCESSES 
FOR THE SIMULATION MODEL 


[step || pRocess PROCESS PROCESS si PROCESS CHARACTERISTICS 
| LRU Failure Arrivals | Computation based on input parameters. | 


Ze Failure Removal A process with unlimited resources and 
random or fixed parameters. 


4, Failure Rectification | A process constrained by number of 
available spares and cannibalization polices. 
4, Base Repair A process with limited or unlimited resources 
and random or fixed parameters. 
». Base-to-Depot A process with unlimited resources and fixed 
Transportation parameters. 
Depot Repair A process with unlimited resources and 
random or fixed parameters. 
4, Number of Demands |} Simulation output (combination of Base and 
Depot demands). 
Number of Simulation output. Depends mainly on 
Backorders - EB.(t) processes 3 and 7. 
Number of Mission- | Simulation output. Depends mainly on 
Capable Systems processes 3 and 8. 
hoe Ao(t) - (SMOE) Computation based on process 9, 
ti 4 ENMCS(t) -(MMOE) | Computation based on process 9 


1. LRU Failure Arrivals 











Each LRU failure in a system causes the system to be Not-Mission- 


Capable (NMC). A time-dependent variable, to be defined as the Mean Time 
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Between Demands (MTBD) for each LRU, is computed based on the inverse of the 
demand rate, D(t), which has a formula given by Equation 3.1. MTBD is not a 
constant parameter since the value of D(t) changes according to the status of the 
utilization rate and the number of Mission-Capable Systems (MCS) at a particular 
time interval of the wartime period. To generate a given interarrival time for each 
LRU of the entire deployment of systems, we sample from an exponential random 


number stream using M7BD as the mean parameter. 


2. Failure Removals 
When the system becomes Not-Mission-Capable (NMC), time is spent on 
isolating the cause of the failure and the removal of the faulty LRU. An exponential 


time or a constant time is used to generate this isolation and removal time. 


3. Failure Rectification 

An NMC system can be recovered by replacing the faulty LRU with a good 
spare. However, the NMC system remains in the same status if there is a shortage 
of spares for that LRU. 

When the policy of no cannibalization is adopted, the model treats the 
failure of an LRU as a system failure (i.e, one system becomes NMC). If there is 
an available spare for that failed LRU in the inventory of stock, the NMC system will 
be restored to a Mission-Capable (MC) status. When no more spares are left, the 
failed system must then wait for a repaired LRU of the same type to be returned to 


the inventory storage facility. Repair prioritization polices that were discussed in 
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subsection II-B5 were used to determine which NMC system is entitled to the newly 
arrived repaired LRU. Additional provisioning of spares to satisfy unfilled 
Backorders is not considered in this simulation in view of the fact that WSM is 
presently only an assessment tool. 

The cannibalization policy is more complex to model. As long as there are 
spares in the inventory, no cannibalization of other NMC systems will be carried out 
to restore a system that just had a failure. In this case, that system is considered MC 
after a fixed time needed for replacing the faulty LRU with a spare. When there is 
no spare left in the inventory, cannibalization will start with the first NMC system. 
A search for the required LRU will be made from the NMC systems until a good 
LRU can be found. Then the system that became NMC last can be repaired using 
a good LRU. If no LRU is found, then the system will remain NMC until a repaired 
LRU is returned to the stock. A limit can be placed on the maximum number of 


systems that can be cannibalized if such a requirement is stated in the WMP. 


4. Repair Flows 
Similar to the analytical approach, a failed LRU can either be repaired at 
the base or sent to the depot for repair depending on its Not-Repairable-This-Station 
(NRTS) value (see Equation 3.24). A uniform(0,1) random stream is used to 


generate the random variable. If the random variable is less than the NATS value, 
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then the failed LRU is sent to the depot for repair. Otherwise, it goes to the base 


for repair. 


a. Base Repair 

If the faulty LRU is to be repaired at the base, the model allows the 
specification of a number of base repair stations in order to study the effect of 
limited base repair on the SMOEs and MMOESs. In the example of a one-base, one- 
depot maintenance system to be studied later in Chapter VI, these repair stations are 
all located at one base and all failed LRUs will form a single queue to be served by 
these stations. Different policies to prioritize repair are also allowed. The base 
repair facility will service a failed LRU, either with an exponential time or a 


constant time, using an input parameter which is specified for that LRU. 


b. Depot Repair 
Repair at the depot requires three time components - fixed 
transportation time from a base to a depot, a fixed or an exponential depot repair 


time and a fixed transportation time from a depot to a base. 


5. Replications for Statistical Significance 
A single replication of the model will cover the entire wartime period. 
Many replications of the simulation must be conducted to achieve satisfactory 
Statistical confidence intervals for the outputs. To measure the results’as a function 


of time, each replication is further divided into smaller equal time intervals. 
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6. Random Number Streams and Generation 
All of the random processes indicated in Table 4.1 require Monte Carlo 
sampling. To minimize the number of random number streams, all the LRUs will 
use the same stream for their arrival generations. The other processes will each have 


a dedicated random number stream. 


7. Repair Prioritization Policies 
An appropriate choice of repair prioritization is important when repair 
resources are limited since a decision to repair one LRU also means a decision not 
to repair another LRU when both are competing for the same repair resources. Two 
widely known policies will be considered. 
a. Furst Come Fust Served (FCFS) 

The FCFS priority selects the LRU which has been waiting the longest 
for repair. Therefore, this prioritization policy does not consider those LRUs that 
may have the greatest demands for spares or repairs. 

b. Least Available Item First (LAIF) 

The Least Available Item First (LAJF) repair priority looks for the 

LRU with the highest EB,(t) to service first since filling a backorder will improve the 


Ao(t). 
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C. IMPLEMENTATION OF WSM USING MODSIM 

MODSIM II is a fourth generation programming language recently acquired by 
the Naval Postgraduate School for developing complex simulation models. Its 
features, programming syntax, and structure are described in the MODSIM Manuals 
(Ref. 4]. MODSIM was chosen to develop the WSM Simulation in view of its 
powerful features to handle an object-oriented discrete event simulation. The WSM 
Simulation is implemented using a PC-based version of MODSIM II. Hence the 
simulation program was limited by the PC base memory of 640K and this was one 
main reason for not being able to incorporate aspects such as condemnation of LRUs 
and an intermediate repair level into the simulation. 

Figure 4-1 depicts the main design flow and linkage between the various 


modules of the WSM simulation program as used in this thesis. 


The program codes and listings are given in Appendix A. 


52 


TGQOW NOILWINWIS WSM HOT FAUNLONULS Nwudsodd 





“NUT PpeTTey Axeu 
abebua OF UOTISIS SI[PI | AtedsyusfTutd Nsw 

> pouaen 
Cqouot3eIs LOaArAO 








pow’ 7p edeuG 


uotTAeTnwts doqys 


pue swtq dsosy 


CS 


T° AUNNSIA 


Sqy[nsel Teoyzasyqeqas Teusd 





qViodayTeusd ASW 


STCAIEAUT 


"3898 YOETTIOD MONQI00ey NSW 


pextjy 32 


: SPOUIEH 
“UOFADITTOSO 
soyasyaeqs 103 IDelqo 


Pon’ 310deuG Pon’ a1OdeuUl 


a IOGOW 
LYIOdaAa 





fqojicodey Loarao 


SOTATATIORS 





uoTJeTNUTS 47e9S 





POW" WSMN 


a TOGOW NIWN 


ATNGON AXOOLS 


*woqsAs UuAMOpP 


g9IojysgelI oq sieds gen AIOQUSAUT IAL ASW 


“ezTTeqyUuURD OY 


swaj3SAS UMOP IOJ YoOoT uoFAIeZTTeqyuuey II3L 


: spouyen 
‘-er1eds 103 BbuTatTemM nw CqouesqyrTaAul LOarEoO 


*suloqsAg 
UMOP 98IOJSEI OQ oIedsS 
IojJ AIOQUPAUT YOIYD 


AIOQUSAUTHYOIYD NASW 
>poureoNW 
EeNINGAIOAVUSAUL LoOarGo 


‘a9ieds yores Ao0jJ sasnand 


POW’ HOOYSI 


POW’ YOOISAd 





Sisqowered jyanduy [je pesy 


BNUYT TTe JO AaGLW eI jnduos 
93eI UOFIBZTTTIN synduod 


AINGON WHO'D! pow sTeqoIOI | POW" BTBqOTOM | 


sqandut pesy —————— 


— LdOad 








*ztedez eyeTdwos 
O23 owt Afyedaez eTdwvg 





BVOFAIBISAIVIS TIAL 


"MmOTZ IPFedeI epyoeq ryedeyesooyud T1gEL 









‘we qeAs peTyej 9107Ese I 
07 FqoueqyAUT 47eIS 





weqsAsyoeudD TIAL 


“guotqoe AyedezZ IoOJ ADeITAd 


“BNUT PETFRI 
JO STBVAFIIV STduVE 





we qeAg 1e3uq TIAL 
: SpouyaNn 








qojuewdtnba Lograo 


pon’ 3dbaa 


sny1 T1e z0J 35efqGo 


Pon’ ydbal 


INdNIGWSuY 
pueweg 


NOLLWZITILN [F8TOOF UOFILTNUTS 
suoyqounj yAeN 
edInosey enend 

6TOOQ soy Ae7yAeIAS 
*eq0efO0/seTqeyTIeA j10OReTEeUeH WOpueY - 

eqo[TH {Te sare[seqg 


POnwTSs 
POWYUIeN 
pond 
POWIBIS 
poppuey 
rabeuew O/I - POWO!L 


:saInpesoig 


SLOSUCHO NI-LIINGA WISGOW 


D. STATISTICAL INFERENCES 

This section describes the statistical analyses that were performed on the inputs 
and the outputs of the simulation runs to verify that these two elements were 
Statistically sound. STATGRAPHICS [Ref. 17] was used to enhance the 


analyses. 


1. The Resulting Distribution for Failure Arrivals 

The arrival of failures for each LRU is generated from an exponential 
random stream using MTBD as the parameter. Since MTBD is not a constant, the 
distribution for the interarrival times of failures for each LRU is not expected to 
belong to any one single family of commonly known distributions such as exponential 
and gamma. To verify this, data for the interarrival times of each LRU were 
collected for statistical analyses to assess whether they could be fitted to any of the 
commonly known distributions. 

An exponential distribution was first used to attempt to fit 3001 data points 
obtained for the interarrival times of LRU D of the example in chapter VI and, at 
first glance, Figure 4-2 suggests that the fit seems to be good. However, on further 
examination using the Kolmogorov-Smirov (K-S) goodness of fit test, the exponential 
distribution was rejected based on the test results given in Table 4.2. A Chi-square 
goodness of fit test was also performed with the same conclusion. The results are 
shown in Table 4-3. Further attempts to fit the data with other distributions such as 


Gamma, Normal and Lognormal were also unsuccessful. 
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FIGURE 4-2: EXPLORATORY ANALYSIS USING EXPONENTIAL 
DISTRIBUTION FITTING ON 3001 DATA POINTS OBSERVED FOR 
THE INTER-ARRIVAL TIMES OF LRU D. 
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TABLE 4.2 : K-S TEST ON INTERARRIVAL TIME DISTRIBUTION 


>. 









Estimated KOLMOGOROV statistic DPLUS 0.0434898 






0.0107071 


Estimated KOLMOGOROV statistic DMINUS 







Estimated overall statistic IN = 0.0434898 





Approximate significance level = 2.3493E-5 
(Therefore not a good fit) 





Lower Upper Observed Expected 
Limit Limit Frequency Frequency Chisquare 


at or below 4.286 4.37296 
4.286 21.429 7.31331 
21.429 38.571 6.64243 
38.571 55.714 3.67673 
55.714 72.857 2.17470 
72.857 90.000 1.00432 


90.000 
107.143 
124.286 
141.429 
158.571 
175.714 
192.857 


107.143 
124.286 
141.429 
158.571 
175.714 
192.857 
210.000 


- 21550 
1.50532 

00728 
1.13197 
2.12192 
3.66202 
6.36820 


210.000 24 .63549 


Chisquare = 64.8321 with 12 d.f. Sig. level = 2.92816E-9 
(Therefore not a good fit) 





All the above are indications that one cannot expect to obtain exact 
analytical expressions for the SMOEs and the MMOEs when the distribution of the 
failure arrivals of each LRU does not belong to a single family of known 


distributions. 
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2.  Replications for Statistical Confidence 

In any simulation, many replications are required so that a reasonable 
statistical confidence can be obtained for all the outputs. In the case of the WSM 
simulation, each replication "opens" at time t=O under the same set of initial 
conditions and then "closes" at the end of the same wartime period, which is fixed. 
Under these conditions, the method of independent replications (see Banks and 
Carson [Ref. 18], p.421-422) is used for this type of terminating simulation. 
The whole simulation was repeated many times, with each replication following the 
conditions described above and using a different random number stream. Therefore, 
all the data points from the replications collected at a particular time interval for a 
particular output are statistically independent and identically distributed. Banks and 
Carson concluded that the classical methods of confidence interval estimation using 
the Student t distribution can be applied to the output. In this simulation, a 95% 


confidence interval was used as the stopping criterion for the whole simulation run. 
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V. WSM BY AN OPUS-8 APPROXIMATION 

Due to the unavailability of a proper analytical tool to assess the sustainability 
of a military capability under wartime environment, the MINDEF logistics staff is 
currently using OPUS-8 to roughly estimate the requirements of war-reserve spares. 
Despite the fact that OPUS-8 is a steady-state model, it remains a popular evaluation 
and optimization tool in MINDEF for assessing the adequacy of spares to support 
a military capability. In the past, the author was personally involved in trying to 
adapt OPUS-8 as a rudimentary sustainability tool but the approach had not been 
successful due to the lack of understanding of the dynamic behavior of a military 
asset under a wartime environment. Having gained understanding of the WSM 
through the development of the analytical and simulation models, the author 
attempts in this chapter to formulate a proper methodology to allow OPUS-8 to be 


used as an approximation tool for the WSM. 


A. THE STEADY-STATE ASSUMPTIONS OF OPUS-8 


1. Current Features and Assumptions 
OPUS-8’s main strength lies in its superb ability to handle complex multi- 
echelon logistics support and multi-indenture system structure. Figure 5-1 illustrates 
a multi-echelon logistics support with multi-system deployment and Figure 5-2 the 


multi-indenture breakdown of a system. 
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FIGURE 5-1: MULTI-ECHELON LOGISTICS SUPPORT 
WITH MULTI-SYSTEM DEPLOYMENT. 


LEGENDS: 


O-LEVEL - Organization Level Maintenance 
I-LEVEL - Intermediate Level Maintenance 
D-LEVEL - Depot Level Maintenance 
ILOC_ - Intermediate-Level Location. 
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FIGURE 5-2 : MULTI-INDENTURE SYSTEM STRUCTURE 





Sophisticated techniques such as network and dynamic programming are 
imbedded in the OPUS-8 optimization algorithms [Ref.4]. The user has a choice of 
different SMOEs (Effectiveness) to be used together with spares investments (Cost) 
as the optimization criteria. A single run will produce a Cost-Effectiveness curve, 
with each optimal point on the curve providing an investment amount corresponding 
to an optimized range and depth of spares. In addition to the optimization 
techniques, OPUS-8 uses many well-known steady-state analytical formulae from 
stochastic and inventory theories based on the following assumptions: 


a. the number of demands in a given time interval follows a Poisson distribution 
with a known constant mean value. However, at the operational level, the 
number of demands in a given time interval is allowed to follow a Geometric 
Compound distribution with a given value on the Variance-to-Mean Ratio 
(VMR). With this type of demand distribution batching of demands can be 
modelled. 


b. The mean values of the Turn Around Times are known. (No other 
assumptions are made on the distributions of the Turn Around Times.) 
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c. A failure of one type of item is statistically independent of those that occur 
for any other type of item. However by introducing Removal Rate Factors 
for item multiple errors in a mother item can modelled. 


d. The Turn Around Times for items in the repair cycle are statistically 
independent. 


e. No queues are assumed to be formed at repair facilities, nor is any controlled 
bulk service permitted in any phase of the repair cycle" [Ref. 4]. 


2. Limitations Of Steady-state Assumptions. 

Steady-state assumptions are fine for peace-time prediction of spares 
requirements since the fluctuations in demands when analyzed over a long-time will 
tend to stabilize. But in time of war when the operational profile of a military asset 
changes dynamically with time, these steady-state techniques are no longer valid. 
Also, the assumptions of no queuing and limited batching are not too suitable for 


sustainability analyses. 


B. OPUS-8’S SUSTAINABILITY OPTION 

All of the options of OPUS-8 are based on steady-state assumptions and they 
would require extensive modifications to achieve the capabilities of the WSM. 
However, the sustainability option’ of OPUS-8 has features which are favorable for 
use aS an approximation to the WSM without modifying any of its codes. 


Unfortunately, its main limitation is that no deployed or on-site repair capability is 


| OPUS-8 uses the term "Endurance" instead of sustainability [Ref. 4]. However the 
meanings of both are equivalent in this thesis. 
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allowed for the systems deployed as a military capability during a given the wartime 
period. 

Similar to the WSM, MOEs such as the Total Pipeline and Backorders for each 
LRU and the expected number of NMCS are obtainable from the sustainability 
option of OPUS-8. Therefore, these MOEs are to be interpreted directly as the 
WSM’s MOEs if approximation is possible. 

When there is no repair, there is only one pipeline generated from the location 
where the systems is deployed. The pipeline will increase linearly with time if the 
demand rate for spares is a constant. In the case of the nonhomogeneous Poisson 
assumption, the pipeline increases monotonically with time during the wartime 
period, with its value dependent on the demand rate at the time of analysis. 
Therefore, at any time t the OPUS-8 option can produce an approximate result to 
the WSM at the same time, by using a weighted-average utilization value of all the 
utilization rates that have occurred over the time from zero to t. This means that 
multiple runs must be performed to analyze the entire anticipated wartime period. 
The runs are ordered so that each subsequent time period overlaps the previous time 
period. The rationale behind this approximation becomes evident when a numerical 
example is used to compare the results between the OPUS-8 approximation, the 


analytical model and the simulation model. These results are given in Chapter VI. 
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V. VERIFICATION OF THE WARTIME SUSTAINABILITY MODELS 


In Chapter III, an analytical model was developed for the WSM using certain 
assumptions to obtain exact analytical formulae. However, when these assumptions 
were relaxed, the distributions of the failure arrivals for the LRU were found not to 
belong to any single family of known distributions as explained in subsection V-D1. 
As such, exact analytical expressions for the SMOEs and MMOEsS were not attainable 
and therefore a simulation model was built for the WSM. 

To verify that the results from the analytical model conform with the results 
from the simulation model, a numerical example using the similar set of input 
conditions was run under both models. This example was also performed with the 
simulation model to study the effects of wartime logistics policies such as 
cannibalization, limited repair resources and repair prioritization. 

A special case of no repair capability during the anticipated wartime period was 
also investigated by both models. This case was also examined by a third approach 
using the sustainability option of OPUS-8 as an approximation for the WSM. 

This chapter covers the comparison of the results in detail and draws 


conclusions about significant trends. 
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A. THE EXAMPLE 
All the analyses mentioned in this chapter use the same numerical example. 


The following input parameters are used. 


1. War and Mobilization Plan (WPM) 

A squadron of twenty-four aircraft are deployed at only one base. The 
WPM in this case study anticipate a 720-hour wartime period with an initial 
deployment of twenty-four of the same aircraft type. Only the surge and post-tension 
activities were considered in this example although peacetime and pre-tension 
activities can also be part of the WMP. This was done to minimize the number of 
factors that can affect the output. 

The surge period begins at time t=0 and lasts for 168 hours. Three sorties 
per day are expected from each aircraft, with each sortie having an average flight 
mission time of 1.6 hours. This corresponds to a daily utilization rate of 20% for 
each aircraft. The scenario then anticipates a daily utilization rate of 10% for the 
rest of the wartime period (post-tension period). An abrupt (step) change was used 
for the transition from the surge period to the post-tension period. This was done 
to simplify the analytical expressions for the analytical model. In this WMP, it is 
assumed the operational planners were mainly interested in estimating the ENMCS(t) 


during the anticipated wartime period. 


de 


System Structure and Logistic Support 


Each aircraft has a one-indenture level breakdown of ten LRUs with input 


parameters as shown in Table 6-1 (not all the parameters were used each numerical 


example). 





INVLVL 


BRTIM 


DRTIM 
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fe aee777-7 || O.21') 2 iL 72.0 HG) A) 1.0 
D 86.9 | 0.59 | 1 5 72.6 TG) | tale. O 
E 222.71] 0.04 Ih 1 2 48.0 1.0)i| 1.0 
F Hod.7 | O.Odal) fl 4 48.0 iO | a. 0 
eo 1754.4 | 0.24 )\2 i 720 1, Gy et. 0 
iH 223.2 || 0.09 | 1 2 48.0 1.07) 1.0 | 
I 58.8 0.06 | 1 6 48.0 | se, 16 
| 5 91.6 |} 0.28) 1 2 | 48.0 } 120.0 1.0) 1.0 | 


TABLE 6-1: INPUT PARAMETERS FOR LRUs 






Mean Time Between Failures in hours. This is the reciprocal of the 
failure rate. This characteristic is constant for each LRU. 


Not Repairable This Station. A value to indicate that the proportion 
of the repair of each LRU which will flow to the depot. The 
complement of this value is the proportion being repaired at the 
base. No condemnation is assumed. 


Quantity Per Application. This is the quantity of each LRU found in 
each aircraft. 


Initial Inventory Stock Level for each LRU. 
Mean Repair Time (hours) at the base repair facility. 


Mean Repair Time (hours) at the depot repair facility. 
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REMTIM = Mean Removal Time (hours) at the system level which includes fault 
isolation time and removal time. This time component was used only 
for the simulation model. 


REPTIM = Mean Replacement Time (hours) at the system level to replace a 
faulty LRU with a good spare. This time component was used only 
for the simulation model. 


(120 hours) 


DEPOT REPAIR 
TIME 


(120 hours) (120 hours) 


BASE TO DEPOT (48 - 72 hours) DEPOT TO BASE 
TRANSPORTATION TRANSPORTATION 
TIME BASE REPAIR TIME 


TIME 


i OPERATIONAL 
oe. 





FIGURE 6-1 : MAINTENANCE STRUCTURE FOR THE EXAMPLE. 
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The logistics support for the aircraft consists of one repair base and a 
repair depot having the repair turn-around times given in Figure 6-1. The 
transportation time from base to depot has a value of 120 hours denoted by TBDTIM 
and transportation time from depot to base also has a value of 120 hours denoted 


by TDBTIM. 


B. RESULTS FROM ANALYTICAL CALCULATIONS 

The wartime period was divided into twenty equal time intervals for the 
analytical calculations. PC-MATLAB [Ref. 5] was chosen because of its powerful 
features in matrix manipulation and ease in plotting the results. It was also preferred 
because it takes much less time to program the procedures as compared to 
FORTRAN. The exact arava expressions derived in Chapter III can be 
expressed in MATLAB almost exactly as they are written in conventional 
mathematical forms. Two cases were used; the first (Analytical Case One) had 
unlimited base repair and depot repair capabilities and the second (Analytical Case 
Two) had no repair capability during the 720-hour wartime period. 

Appendix B presents the listing of the program which implement the analytical 
model for Analytical Case One in PC-MATLAB using the computational steps 
described in the following subsections. To illustrate the process, the computational 
details were carried out for one LRU, LRU-D. The results of all other LRUs were 


computed in the same manner using PC-MATLAB. 
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1. Computing the Demand Rates and Repair Flows 
The ns(t) used in the analytical calculation of the demand rates based on 
Equation 3.1 is assumed constant at a value of 24. This is a worst-case calculation 
for the demand rates. The distribution of the failure arrivals for each LRU follows 
a nonhomogeneous Poisson process with only two demand rates based on the 
utilization rates. 
Based on Equation 3.1, the peak demand rate, D,, for LRU-D is 


D, = fxqxu(t) x 24 


= (1/86.9)x1x0.2x24 = 0.055 failures/hour. (6.1) 
The lower demand rate, D, is 
D, = (1/86.9) x 1x 0.1.x 24 = 0.028 failures/hour. (6.2) 


From Equation 3.24, D, has two components as follows 


D’, = D,.(1-NRTS) (flow to base) 

= 0.055 X 0.41 = 0.023 failures/hour, (6.3) 
D*, = D,.NRTS (flow to depot) 

= 0.055 X 0.59 = 0.032 failures/hour. (6.4) 


D, also has two components as follows 


D°’, = D,.(1-NRTS) (flow to base) 
= 0.028 X 0.41 = 0.011 failures/hour, (6.5) 
D°, = D,.NRTS (flow to depot) 


0.028 X 0.59 = 0.016 failures/hour. (6.6) 
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2. Calculating the Base Repair Pipeline 
From Table 6-1, the base repair time for LRU-D has a constant value of 
BRTIM = 72 hours. A time chart depicted in Figure 6-2 is used to indicate that the 
computation of the expected Base Repair Pipeline, A°(t), for LRU-D is time 
dependent. Equation 3.5 is used for all the computations of A°(t) where the term /1 


- G(s,t)] has a value of one for all the calculations. This is due to the fact that 


BRTIME is a constant. 





In time region 1 where {0 < t < 72}, A°(t) increases linearly with time since all the 
failed units of LRU-D are still undergoing repair at the base at a rate of D”,. 
Hence, 


t 


A(t) = [> ds = D?.t = 0.0234. (6.7) 
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In time region 2 where {72 < t < 168}, only failed units of LRU-D that are in repair 


from time {t-72} to time ¢ are still in repair at time ¢ at a rate of D®,. So, 


5 


7 2 eGo (6.8) 


t 
A(t) = [> oD 
i: 
For time region 3 where {168 < t < (168+72)}, the demand rate changes to D,”. 
However, there is still a residual number of units in repair resulting from the 
influence of D,°. Therefore, 
168 +72 


t 
A(t) = | Do ds + [ D3 ds 
1 


= D® (168+72 -t) + D3(t ~ 168) = 3.672 - 0.012r. (6.9) 


In time region 4 where {168+72 < t < 720}, only D,” is causing failures of LRU-D. 
In the same manner as Equation 6.8, only failed units that are in repair from time 


{t-72} to time ¢ are still in repair at time t. Therefore, 
t 
A(t) = [ D3 ds = D3.72 = 0.792. (6.10) 
t- 
3. Calculating the Total Depot Pipeline 
The combined turn-around time at the depot is also a constant with a 


value equal to the sum of 


TBDTIM + TDBTIM + DRTIM = 120 + 120 + 120 = 360 hours. 
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Using the same approach described in the previous subsection, a time chart is also 
used to illustrate the time dependent computation of the Total Depot Pipeline or the 
equivalent term A°(t) which was specified by Equation 28. This is shown in Figure 


6-3. 





FIGURE 6-3 : TIME CHART FOR COMPUTATION OF THE TOTAL DEPOT PIPELINE 
FOR LRE-D. 


In time region A where {0 < t < 168}, A°(t) increases linearly with time since all the 
failed units of LRU-D is still being transported to or in repair at the depot under the 
peak demand rate D*,. Hence, 


t 


MO) = | “ ds = Do 1 . (6.11) 


In time region B where {168 < t < 360}, a portion of A°(t) is increasing linearly with 
time under the influence of D*,. The other portion due to D°, has to wait for repair 


with a time constant of 168 hours since no depot repair is possible. 


ol 


Therefore, 


168 


t 
vee d . d 
A(t) = A ds ie 


= D%168 + DS(t - 168) = 2.688 + 0.0161 . (6.12) 


For time region C where {360 < t < 360+168}, the portion of A°(t) under the 
influence of D®, is still waiting for repair since it only started to accumulate failures 
from 168 hours onwards. On the other hand, depot repair is now possible for other 
portion due to D°, which started from time zero but, like before, there is a residue 


number under repair for a time period of {360+168 - t}. So, 


360 +168 t 
AX%(t) = | Do ds + [25 4 
1 
= Dw1.(360+168 - t) + DS.(t - 168) = 14.208 - 0.016 . (6.13) 


For time region D where {360+168 < t < 720}, only D,° is causing failures of the 
LRU-D. Again only failed units of LRU-D from time {t-360} to time ¢ will still in 


repair at time t. Therefore, 


t 
A%(t) = [ 4 = DS .360 = 5.76. (6.14) 
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4. Calculating the Total Pipeline 
Since there are only one base repair and one depot repair with no 
inventory stocking at the depot, the Total Pipeline is just the sum of the Base Repair 
Pipeline and the Total Depot Pipeline which is 


no] A(t) + A°G) . (6.15) 


The results are discussed later in section E. 


5. Calculating Expected Number of Backorders 
After A(t) were computed for each LRU (i.e., A((t)),the expected number 
of backorders, EB.(t), were calculated based on Equation 3.35. The results are 
discussed later in section E. 
6. Calculating Operational Availability 
The Operational Availability, Ao(t), calculation without cannibalization is 
based on Equation 3.39, and with cannibalization based on Equation 3.47, denoted 
by Ao,(t). Analyses of these results are presented later in section E. 
7. Calculating ENMCS(t) 
The ENMCS(t) calculation without cannibalization is based on Equation 
3.48 and ENMCS,.(t) for cannibalization is based on Equation 3.46. Analyses of these 


results are presented later in section E. 


is 


8. Calculating PEENMCS < NMCS) 
Although not to be used for comparison with the simulation model, the 
alternate MMOE, P(ENMCS < NMCS), has the same calculation for both 
cannibalization policies, based on Equation 3.44. The target NMCS is set at a value 


of 4. However there will be no analysis of these results. 


9. Computations for Analytical Case Two 

The computational steps for Analytical Case Two are straight forward 
Since each LRU has only one pipeline which increases linearly with time. Essentially, 
the pipeline is calculated by multiplying a demand rate and the time period of time 
zero to time t. The peak demand rate D,, computed using Equation 6.1, is used for 
any time t within time region {0 < t s 168} and for time region {168 < t < 720}, the 
lower demand rate D,, based on Equation 6.2, is used. Graphical and numerical 
results of the Total Pipeline and the Backorders for this case are almost the same as 
those obtained under the OPUS-8 Approximation (see Figures 6-6 and 6-13). 


Therefore these results are not repeated here. 
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C. RESULTS FROM SIMULATION RUNS 
Five simulation variations were set up for the example. These variations use 
an input format similar to the one shown in Table 6-2 which list the specific input 


parameters for Simulation Variation Four. 


random seed for base mean repair time 

random seed for interarrival time of failures 
random seed for ratio of NRITS 

random seed for MITR (fault isolation time) 
random seed for depot mean repair time 


5000 = number of replications for each simulation 


20 = number of time intervals (36 hours per interval) 
5 = number of repair stations at the base 


24 = number of aircraft initially deployed for war 


120 transportation time from base to depot 
120 transportation time from depot to base 


= aircraft utilization rates : FROM TO (hours) RATE 
0 168 0.2 
168 720 0.1 


4 = maximum number of aircraft allowed to be cannibalized | 





Table 6-3 provides a summary of the five variations and their differences. Details 
of each variation are given in the following subsections. The results from each 


variation are discussed later in section E. 
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TABLE 6-3: SUMMARY OF DIFFERENCES FOR THE 
FIVE SIMULATION VARIATIONS 


SIMULATION RUNS | MAIN DESCRIPTIONS 


- Exponential interarrival times for failures of 
all LRU. 

- Deterministic values for all repair and trans- 
portation times. 

- No repair capability. 

- No cannibalization allowed. 











Variation One 






















Variation Two - Same as Variation One except that there is 
unlimited repair capability at the base and 


the depot. 






Variation Three - Same as Variation Two except that repair times 
are now made exponential. 


- cannibalization policy is used. 






Variation Four - Same as Variation Three except that the base 
repair resources are limited. 


- FCFS repair priority is used. 





Variation Five - Same as Variation Four except that FCFS is 
replaced by LA/F repair priority. 


1. Simulation Variation One 

This variation is the simplest simulation case of all the five since no repair 
is allowed. The purpose of this variation was to verify the adequacy of using OPUS-8 
as an approximation to the WSM. This variation also helped to verify the algorithms 
of the simulation model. All the input parameters except for the exponential random 
generation of failure arrivals of the LRUs were modeled as deterministic to match 
the characteristics of OPUS-8 input requirements. This also meant that s(t) is kept 
constant for the computation of D(t). In the simulation, a failed LRU caused the 


system to be NMC. Then, if there were spares available, the NMC system was made 
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McC. When no more spares were left, the systems became NMC and when all the 
twenty-four aircraft were made NMC, the replication was considered complete and 


the simulation went to the next replication. 


2. Simulation Variation Two 

The purpose of Variation Two was to ascertain whether the results by 
simulation were close to the results of Analytical Case One when the input 
parameters were kept the same for both. To conform with the assumptions of 
Analytical Case One, the repair resources at both the base and the depot were 
assumed to be unlimited and the policy of no cannibalization was adopted. No 
distribution assumptions were possible for the failure arrivals of the LRUs since the 
value of D(t) was expected to fluctuate widely because ms(t) was allowed to vary 
during the simulation. Indeed, these factors caused differences in results between 
Analytical Case One and Variation Two. Except for the above, the characteristics 


of the other inputs are the same as in Variation One. 


3. Simulation Variation Three 
Variation Three has the same considerations as Variation Two except that 
repair times at the repair base and repair depot were generated from exponential 


distributions. Also, the policy of cannibalization is implemented. 


4. Simulation Variation Four 
Variation Four extends Variation Three to incorporate the realistic 


problem of limited repair resources during an intense conflict period. The simulation 


i 


constrained the repair base to seven repair stations but the repair depot is unlimited. 


FCFS repair priority was used in this variation. 


5. Simulation Variation Five 
The last variation examines the effect of a different repair priority policy 
than Variation Four. Here the LA/F repair priority is implemented; defective LRUs 


with the largest number of backorders have first priority in the base repair queue. 


D. RESULTS FROM OPUS-8 APPROXIMATION 
It was found that the OPUS-8’s sustainability option provided an approximate 
result for the WSM only under the condition that there is no repair capability during 
the anticipated wartime period. However, as explained in Chapter V, a number of 
runs were required to model the anticipated wartime period since each OPUS-8 run 
can only provide the results for one time period and one utilization rate. In the 
example, where 20 time intervals were used to divide the wartime period, the first 
run presented the time period from time zero to the ending time of the first time 
interval, the second also started from zero but ended at the ending time of the 
second interval, and so on up to the twentieth run. In this way, the Total Pipeline 
increased with the run number since no repair was allowed. The approximation 
also required that the utilization rate for each time period be weighted by the 
amount of time spent in the first and second utilization rates. 
The required number of runs and the corresponding time panee and 


utilization rate are shown in Table 6-4. Except for these two parameters, all other 
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input parameters were the same for the twenty OPUS-8 runs. The input format and 
data requirements for the OPUS-8’s sustainability option are provided in Appendix 
E only for run number 10 to illustrate the approximation process. __Results are 


presented later in section E. 


TABLE 64 : OPUS-8 RUNS TO APPROXIMATE WSM. 
LENGTH OF 
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E. ANALYSES OF RESULTS 


This section shows the pertinent results obtained from the three methods 


described in the previous sections of this chapter and specifically gives the analyses 


of the outcomes on the Total Pipeline, Backorders, Ao(t), and ENMCS(t). 


To enhance the analyses, the following legends are used in all the graphs 


presented this section: 
O : 


Al 


S1; 
52 : 
S3 ; 
S4: 


S5: 


OPUS-8 Approximation 

Analytical Case One without Cannibalization 
Analytical Case One with Cannibalization 
Analytical Case Two 

Simulation Variation One 

Simulation Variation Two 

Simulation Variation Three 

Simulation Variation Four 


Simulation Variation Five. 


1. Analyses on Total Pipeline 


For the purpose of selecting a few LRUs for in-depth analysis, the Total 


Pipeline results for all the ten LRUs for the analytical model, the simulation model 


and OPUS-8 Approximation are first examined. Figures 6-4, 6-5 and 6-6 provides the 


these results graphical for Analytical Case One (A1), Simulation Variation Four (S4) 
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and OPUS-8 Approximation (O) respectively. The corresponding numerical values 


for A2 and S4 are tabulated in Tables C-1 and D-1. Numerical results from O are 


identical to A2. 
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FIGURE 6-4: RESULTS OF TOTAL PIPELINE FOR ALL LRUS 
UNDER ANALYTICAL CASE ONE. 
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FIGURE 6-5: RESULTS OF TOTAL PIPELINE FOR ALL LRUS 
FROM SIMULATION VARIATION FOUR. 
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FIGURE 6-6: RESULTS OF TOTAL PIPELINE FOR ALL LRUS 
FROM THE OPUS-8 APPROXIMATION. 


LRU-D, LRU-I and LRU-J are observed to have the largest values. As 
more evident in Figure 6-4, the Total Pipelines of LRU-D and LRU-I peaked at quite 


different times and displayed other contrasting characteristics. Also, these two LRUs 
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displayed other different behavior in their Backorders results, which are shown in 
Figures 6-11 and 6-12 of the next subsection. Therefore from now on, only these two 
LRUs are further analyzed. 

First we consider the input characteristics of these LRUs for some 
explanations of their differing results. As given in Table 6-1 under the NRTS 
column, LRU-D had a high proportion (59%) of its failed units flowing to the depot 
for repair whereas the proportion for LRU-/ is only 6%. Other differences are; 
MTBF of 86.9 hours for LRU-D and 58.8 hours for LRU-D, mean base repair time 
with LRU-D needing 72 hours and LRU-I needing 48 hours, and initial inventory 
levels of 5S and 6, respectively. The rest of the inputs were the same. 

Figure 6-7 shows the results of the Total Pipeline for LRU-D from the 
three versions of the WSM under the condition that there was no repair capability 
during the whole anticipated wartime period. This can be considered the worst case 
scenario in which all the repair resources were severed by the enemy at the onset of 
war. All the results agree extremely well. It is understandable that all the results 
of the Total Pipeline are monotonically increasing since when repair was not allowed, 
more and more failed units of LRU-D became unserviceable as the war progressed. 
The extremely close outcomes suggest that both the analytical model and the 
simulation model are sound and also that OPUS-8 Approximation is a reasonable 
approach. 

Figure 6-8 shows similar trends for LRU-/ although the ordinate values are 


larger. This was mainly due to LRU-J having a lower MTBF and, as a consequence, 
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a higher inherent failure rate than LRU-D. The other differing inputs mentioned 


earlier had no effects since no repair was involved. 
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FIGURE 6-7: RESULTS OF TOTAL PIPELINE FOR LRU-D 
FOR MODEL VARIATIONS HAVING NO REPAIR CAPABILITY. 
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FIGURE 6-8: RESULTS OF TOTAL PIPELINE FOR LRU-I 
FOR MODEL VARIATIONS WITH NO REPAIR CAPABILITY. 


More interesting trends were observed for those model variations having 
a repair capability. Figure 6-9 shows five different results for the Total Pipeline of 


LRU-D. A1 and S2 had more abrupt changes than the others. This is due to the 
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deterministic repair times being used at both the base and the depot. When 
exponential repair times were used the results showed more gradual changes as 
evident from S3, S4 and SS. 

The highest value of Total Pipeline came from A1 which peaked at a value 
of 9.42 at 360 hours. The high proportion of failures flowing to the depot for repair 
is probably the dominant cause of this peak since the total depot turn-around time 
was 360 hours. S2 had a similar profile for the same reason. However, for S4 and 
S5 the limited repair resources at the repair base result in failures having to queue 
for these resources. Their peaks occurred consistently at 252 hours which is an 
indication of a repair "bottleneck" caused by the base repair queue. The number in 
the queue is observed to dominate the number in the repair depot. 

Among the three cases with unlimited base repair capability, the lowest 
range of values between t=0 and t=500 hours was observed when the policy of 
cannibalization was adopted for S3. This indicates that the cannibalization policy is 
an effective means of minimizing the number of demands for spares since good parts 
from NMC systems can be cannibalized as spares. For the two cases with limited 
base repair capability and under the policy of no cannibalization for LRU-D, the 
LAIF priority (S5) has a better result than the FCFS priority (S4) although the 
difference was not as significant as the cannibalization vs. no cannibalization 
comparison. This is a consistent result since the squadron is expected to sustain 


better when its worst LRUs are repaired first. 
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FIGURE 6-9: RESULTS OF TOTAL PIPELINE FOR LRU-D 
FOR MODEL VARIATIONS WITH REPAIR CAPABILITY. 


On the whole for the simulation cases, S3 has the lowest Total Pipeline 
which is evident from the fact that there was unlimited repair capability and parts 


can be cannibalized as spares. 
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All the model variations in Figure 6-9 are observed to reach steady-state 
at about 550 hours. One might say that this is because the lower utilization rate 
became the dominant rate by this time period. This would also cause the value of 
ns(t) be to become stable. Therefore, when computed from Equation 3.1, the failure 
arrival rates of the LRUs become more constant. It has been shown in Equation 3.4 
that when the failure arrival rate of an LRU is constant and follows a Poisson 


process, then the 7Jotal Pipeline is in steady-state. 


Certain results were observed for LRU-/ that differ from those in LRU-D, 
as Shown in Figure 6-10. The most obvious change is that the peak has shifted to 144 
hours for Al and S2. The peak for the other three cases is now at 180 hours. The 
most convincing reason for this result is the fact that a very high proportion (94%) 
of LRU-/ failures went to the base for repair where the demands for spares were 
greatest at about 168 hours. 

Another major change is the value of the peak itself. The highest peak 
value now belongs to S4 closely followed by S5. The results from Al are much 
lower. These observations are opposite to what were observed for LRU-D. These 
can be explained by the same fact that 94% of all LRU-/ failures went to the base 
for repair as compared to LRU-D’s 6%. Also S4 and S5 had limited number of base 
repair resources (5 repair stations). Therefore their Total Pipeline are inflated by the 
number in queue at the base repair. S4 has a larger number than SS since the latter 


used a better repair prioritization policy (i.e., LA/F). 
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FIGURE 6-10: RESULTS OF TOTAL PIPELINE FOR LRU-I 
FOR MODEL VARIATIONS WITH REPAIR CAPABILITY. 
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2. Analyses Of Backorders 

As an illustrations of Backorders results obtained by the three different 
versions Of WSM , results for all the ten LRUs for Al, S4 and O are presented 
graphically in Figures 6-11, 6-12 and 6-13, respectively. The corresponding numerical 
values for Al and S4 are given in Tables C-2 and D-2. Results from O are almost 
identical to that of A2 which is shown later in Figure 6-14. In particular for A2 as 
shown in Figure 6-11, only LRU-D and LRU-J have backorders although these two 
and LRU-I experienced high numbers in the Total Pipeline. Also in the same graph, 
no backorders were observed for the other LRUs although backorders were observed 
for all the LRUs in the case of S4 and O as shown in the other two figures. The 
analysis of these differences are discussed below, using only LRU-D and LRU-/ for 
comparison. 

Without repair capability, Figures 6-14 and 6-15 illustrates results of the 
Backorders from model variations A2, S1, and O for LRU-D and LRU-I, respectively. 
The closeness of these results also reinforce the suggestions made for the Total 


Pipeline results . 
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FIGURE 6-11; RESULTS OF BACKORDERS FOR LRUS 
UNDER ANALYTICAL CASE ONE. 
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FIGURE 6-12: RESULTS OF BACKORDERS FOR ALL LRUS 
FROM SIMULATION VARIATION 4. 
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FIGURE 6-13: RESULTS OF BACKORDERS FOR ALL LRUS 
FROM THE OPUS-8 APPROXIMATION. 
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FIGURE 6-14: RESULTS OF BACKORDERS FOR LRU-D 
FOR MODEL VARIATIONS WITH NO REPAIR CAPABILITY. 


o> 


: 


Wartime Period (Hours) 





FIGURE 6-18: RESULTS OF BACKORDERS FOR LRU-I 
FOR MODEL VARIATIONS WITH NO REPAIR CAPABILITY. 


As betore, more interesting results came from those cases with repair 
capability. LRU-D which has a comparatively large number of units in its Total 
Pipeline as evident from Figure 6-9 of the previous subsection. In fact, Al, S2, $3, 


S4 and S§ have Total Pipeline of at least a value of § when time is more than 100 


Yo 


hours. Although LRU-D has a comparatively high initial inventory level of 5, these 
spares are insufficient to meet a Total Pipeline of more than 5. As a consequence, 
Backorders are expected to occur early in the wartime period and Figure 6-16 


indicates that this is true. 
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FIGURE 6-16: RESULTS OF BACKORDERS FOR LRU-D 
FOR MODEL VARIATIONS WITH REPAIR CAPABILITY. 
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In fact, the Backorders for the five cases are higher than for any of other LRUs. 
Other than this, the Backorders profiles are similar to the corresponding Total 


Pipelines profiles of LRU-D for the same reasons given in the previous subsection. 


In the case of LRU-J, only S4 and SS5 have peak Total Pipelines that exceed 
a value of 6 as indicated in Figure 6-10 of the previous subsection. With an initial 
inventory level of 6, the Backorders results depicted in Figure 6-17 showed the 
expected outcome that only these two cases have significant Backorders while the 
other cases displayed negligible Backorders. Hence, these results suggest that for any 
LRU, the Backorders depends mainly on the Total Pipeline and the initial inventory 


level. 


98 


ad 
o 


5 
a 
2 
© 
i 
S 
S 
o 
ao 
S 
B 
E 
= 
= 
0 
i 


Wartime Period (Hours) 





FIGURE 6-17; RESULTS OF BACKORDERS FOR LRU-I 
FOR MODEL VARIATIONS WITH REPAIR CAPABILITY. 
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3.) Analyses on Operational Availability, Ao(t). 

Under the policies of no repair capability and no cannibalization, the Ao(r) 
results obtained trom A2, Sl and O are shown in Figure 6-18. These Ao(r) results 
are Observed to reach a value of 0.0 at about 450 hours from the start of the war. 
These results retlect the tact that more and more deployed systems became NMC as 
time progressed under the conditions of no repair and limited amount of initial 
inventory supply. Also, the results from the analytical model and the simulation 
model are quite simular up to about 150 hours. The values from S1 change more 
smoothly since they are average values weighted by at teast $000 replications. On 
the other hand, the curve produced by the OPUS-8 Approximation has a more 
uneven protile and an oftset as compared to the other two. The unevenness 1s 
understandable since the results were obtained from an approximation method of 
twenty OPUS-S runs. Nevertheless, they all converged to a value of zero at about 
450 hours. This is a good indication that the three versions of the WSM have similar 


asymptotic characteristics. 


The Ao(t) results for model variations with repair capability are depicted 
in Figure 6-19, An additional result from case Al is obtained by introducing a 
cannibahzation policy and using Equation 3.47, 

The towest 4o(t) value for all six cases occurs at between 300 and 400 
hours from the start of the war. This can be explained by the fact that LRU-D and 


LRU-J, which have the most Backonlers among alt the LRUs (see Figures 6-11 and 
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6-12) peaked at between 200 and 350 hours. The possible causes of this can be 
attributed to effects such as 

a. a long total depot turnaround time of 360 hours for all cases, 

b. limited base repair resources for S4 and S5, and 

c. relatively high demands for spares but insufficient inventory level for LRU-J. 
Therefore, it can be deduced that short repair turnaround times plus adequate repair 
resources, reliable LRUs and sufficient inventory stocks for these LRUs are 
fundamental requirements of a logistics support policy to sustain a military capability 
during the wartime period. 

The results also indicate that Al with cannibalization, or Alc, has a better 
outcome than not allowing cannibalization. 

As observed previously in Jotal Pipeline and Backorders, steady-state 
characteristics for all simulation cases are also noticed near the end of the wartime 
period. Analytical cases stabilized earlier at about 550 hours which is similar to the 
results for Backorders. 

Overall, the Ao(t) results of all the cases are quite close in value. This is 
a good indication that the formulae derived for the analytical model and the 
algorithms used for the simulation model are comparable. However, more complex 
examples should be investigated to draw more definitive conclusions and this is 


recommended as a subject for further development for the WSM. 
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FIGURE 6-18: RESULTS OF OPERATIONAL AVAILABILITY (Ao(t)) 
FOR MODEL VARIATIONS WITH NO REPAIR CAPABILITY. 
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FIGURE 6-19: RESULTS OF OPERATIONAL AVAILABILITY (Ao(t)) 
FOR MODEL VARIATIONS WITH NO REPAIR CAPABILITY. 
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4. Analyses on ENMCS(t) 
It has been shown in the analytical model that ENMCS(t), whether with 
or without cannibalization, is related to Ao(t) as shown in Equations 3.46 and 3.48, 
respectively. The same formulae are also appropriate and were used to compute the 
ENMCS(t) in the simulation model. As a consequence, the results shown in Figures 


6-20 and 6-21, display "mirror images" of the Ao(t) curves. 
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FIGURE 6-20: RESULTS OF EXPECTED NUMBER OF NOT-MISSION- 
CAPABLE SYSTEMS (ENMCS(t)) FOR MODEL VARIATIONS 
WITH NO REPAIR CAPABILITY. 
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FIGURE 6-21: RESULTS OF EXPECTED NUMBER OF NOT-MISSION- 
CAPABLE SYSTEMS (ENMCS(t)) FOR MODEL VARIATIONS 
WITH REPAIR CAPABILITY. 
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VII. SUMMARY, CONCLUSIONS AND RECOMMENDATIONS 


A. SUMMARY 

The development of the Wartime Sustainability Model (WSM) for MINDEF 
began with a discussion of appropriate assumptions, policies and MOEs. These 
provided the framework for the design of both a analytical model and a simulation 
model for the WSM. 

Exact analytical expressions were then derived for the Analytical WSM under 
several scenarios such as a nonhomogeneous Poisson process for the failure arrivals 
for the LRUs and unlimited repair resources at the repair facilities. The analytical 
model was then used on a simple example with two-echelon logistics support 
Organization and one-indenture system structure. However, this analytical model is 
applicable for two-echelon two-indenture problems. 

A simulation version of the WSM was also developed to handle variations on 
repair policies such as limited repair resources and two repair priorities. These 
variations are extremely difficult to model analytically. 

Finally, OPUS-8, a multi-echelon multi-indenture model developed by Systecon 
AB, was used as an approximation to the analytical version of the WSM for the case 
of no repair capability. No repair capability is a limitation of the OPUS-8 


Sustainability option. 
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When no repair capability and no cannibalization are used in an example, the 
results of Total Pipeline, Backorders, Ao(t) and ENMCS(t) are very close for the three 
versions of WSM. Limited repair capability at the base had the most adverse effect 
on these measures. Adverse effects were also caused by LRUs having high demand 
rates but insufficient initial inventory levels. 

The policy of cannibalization performed better than one without and similarly, 
LAIF repair priority was shown to improve the MOEs slightly when compared to a 
FCFS policy. 

The results also exhibited steady-state characteristics which began at about 450 


hours from the start of the war. 


B. CONCLUSIONS 
Based on the work which was carried out in this thesis, the following 
conclusions are made: 


1. The analytical model is preferred to the simulation model if the underlying 
assumptions are acceptable. The exact algebraic expressions permit ease of 
computation and tractability when applied to two-indenture and two-echelon 
problems. The results from the example indicate that the analytical model 
and the simulation model for the WSM have comparable outcomes. 


2. The simulation model allows one to consider alternative complex repair 
policies to sustain a military capability under wartime environment. Very few 
assumptions are needed in its use and therefore it is more Suitable for 
problems for which exact mathematical formulae can not be derived. 
Unfortunately, this approach is time-consuming especially for multi-indenture 
and multi-echelon problems. 


3. The use of OPUS-8 in its present form as an approximation to WSM is only 


valid when analyzing cases which have no repair capability. However, it is 
capable of handling a complex military capability compmising different types 
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of systems, of which each system can have many LRUs. This capability is not 
available in the current WSM. 


4. Finally, the numerical results indicate that the SMOEs and MMOEs are 
significantly affected by the interactions between spares allocation, repair 
resource allocation and repair policies such as cannibalization and repair 
prioritization. 


C. RECOMMENDATIONS 
The following are recommended for further development: 


1. The simple version of the WSM developed in this thesis should be expanded 
into a general multi-echelon, multi-indenture structure. With a more 
complicated version, an analytical model may not be possible. For example, 
the present assumption of independence of all failures should probably be 
relaxed. In addition, the WSM input parameters need to be varied so that 
analyses of the results can be more exhaustive. Since the simulation model 
will undoubtedly be required for such analyses, a more powerful version of 
MODSIM should be used for this purpose. The present WSM program 
written in MODSIM provides a suitable basis for such expansion. 


2. The ultimate WSM will involve optimization . Such a model should be 
analytical if possible because exact algebraic formulae can be more readily 
used in the optimization process. Optimization techniques such as network 
and dynamic programming should be investigated. 


3. OPUS-8 remains the standard spares optimization tool for MINDEF 
personnel until some version of the WSM is fully implemented. It is 
recommended that the developer of OPUS-8 should consider the expansion 
of its current steady-state model into one which can be used to optimize the 
wartime requirements of spares. 
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APPENDIX A. MODSIM-IIT PROGRAM FOR WSM SIMULATION. 


ei ee i 
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{ WARTIME SUSTAINABILITY MODEL (WSM) SIMULATION WRITTEN IN 
MODSIM-II BY LIM HUNG HENG FOR THE COMPLETION M.SC (OA) 
THESIS AT NAVAL POSTGRADUATE SCHOOL } 
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The WSM Simulation implemented by MODSIM has the following modules: 

MWSM.MOD - WSM is the main module. It declares new instances of objects, 
controls the other modules and manages replications by checking 
for Confidence Limit convergence and after each replication, it 
disposes the appropriate objects to clear computer memory for 
the next replication. 


DGLOBALS.MOD AND 

IGLOBALS.MOD- These are the definition and the implementation modules for 
the Globals module respectively (all modules except the main 
have a definition module and a implementation module). The 
Globals module contains the declaration of all the global 
variables used by the program and have three global procedures 
UTILIZATION, Demand and READINPUT. UTILIZATION is 
ensure the use of the correct utilization rate to compute the 
MTBD and Demand is the actual computation of the M7BD for 
any LRU. READINPUT handles all the input data supplied by 
the user. 


DREPAIR.MOD AND 

IREPAIRMOD - The Repair module contains the Station object to manage the 
base repair queuing and servicing. Different policies in repair 
prioritization will be implemented here. 


DEQPT.MOD AND 

IEQPT.MOD - The Egpt module contains the generic Equipment object to 
manage all the arrivals of faulty LRUs and the decision for base 
or depot repair. Different policies in repair prioritization will 
also be implemented here. 


110 


DSTOCKMOD AND 

ISTOCKMOD - _ The Stock module contains the Invitem object to manage the 
recovery of a system either with the use of spare or through 
cannibalization of previous failed systems. 


DREPORT.MOD AND 
ITREPORT.MOD - The Report module collects statistics at fixed interval of time 
during each replication so that statistical convergence can be 
managed. The StatMod module provided by MODSIM is used 
heavily to collect the statistics. It can report the current or the 
final statistics for all the required measures and then dispose of 
all used memories. 
{PROGRAM STARTS 
{ Definition Module for Globals} 
{ Contains data structures which are importable to other modules. 
Objects and variables declared here are essentially the globals 
of the program. } 
DEFINITION MODULE Globals; 


{ Now import objects and variables from built-in routines and WSM 
specific objects } 


FROM IOMod IMPORT StreamObj, File UseType; 
FROM RandMod IMPORT RandomObj; 

FROM StatMod IMPORT SINTEGER; 

FROM Stock IMPORT InventoryQueueObj; 

FROM GrpMod IMPORT QueueObj; 

{ Now declare own objects, procedures and variables } 
ary PE 


TimeArray = ARRAY INTEGER OF REAL; 
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PROCEDURE UTILIZATION(IN Tchange, 
Trate : TimeArray; 
IN NoOfRate : INTEGER) : REAL; 
{ UTILIZATION ensures that the correct utilization is being used 
to compute the demand rate of each LRU } 


PROCEDURE Demand(IN identity: INTEGER; 
IN rate : REAL; 
IN system : INTEGER) : REAL; 
{ Demand uses the utilization rate, the MTBF, the current 
number of available systems and the QPA to compute the 
mean time between arrivals of the demands } 


PROCEDURE READINPUT; 
{ READINPUT opens all appropriate input and output files and read 
in all the essential input parameters } 


VAR 


IdleStationQueue : QueueObj; 
{ Queue object to track number of idle base repair stations } 


ServiceQueue : QueueObj; 
{ Queue object to enqueue and track the EquipmentObj waiting 
for repair } 


InventoryQueue : ARRAY INTEGER OF InventoryQueueObj; 
{ Queue object to enqueue and track the EquipmentObj waiting 
for spares } 


SystemStat : ARRAY INTEGER OF SINTEGER; 
{ statistical variables to collect weighted statistics 
for the number of available systems at specified time intervals } 


BaseDemandStat : ARRAY INTEGER, INTEGER OF SINTEGER; 


{ statistical variables to collect weighted statistics for the 
number of Total Average Pipeline at specified time intervals } 
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InventoryStat : ARRAY INTEGER, INTEGER OF SINTEGER; 
{ statistical variables to collect weighted statistics for the 
number of Expected Backorders at specified time intervals } 


System : INTEGER; 
{ number of available systems at any time } 


Inventory : ARRAY INTEGER OF INTEGER; 
{ number of available spares at any time } 


SystemDown : ARRAY INTEGER, INTEGER OF INTEGER; 
{ number of unavailable systems at any time } 


BasePipeNo : ARRAY INTEGER OF SINTEGER; 
{ number of total average pipeline at any time } 


DepotPipeNo: ARRAY INTEGER OF INTEGER; 
{ number of depot pipeline at any time } 


InvLevel : ARRAY INTEGER OF INTEGER; 
{ number of initial stock for each LRU } 


NoOfltem : ARRAY INTEGER OF INTEGER; 
{ number of LRU in a system } 


ItemName : ARRAY INTEGER OF STRING; 
{ name of each LRU } 


ItemArrmean : ARRAY INTEGER OF REAL; 
{ MTBF of each LRU } 


ItemDptmean : ARRAY INTEGER OF REAL; 
{ mean depot repair time of each LRU } 


ItemBasmean : ARRAY INTEGER OF REAL; 
{ mean base repair time of each LRU } 


ItemRemtime : ARRAY INTEGER OF REAL; 
{ mean removal repair time of each LRU } 


ItemReptime : ARRAY INTEGER OF REAL; 
{ mean replacement repair time of each LRU } 
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ItemNtrs : ARRAY INTEGER OF REAL; 
{ NRTS of each LRU } 


{ input streams } 
InputFile : StreamObj; 
SystemFile : StreamOb); 
SSystemFile : StreamObj; 
SInventFile : StreamOb); 
SBaseFile : StreamObj; 


{ random streams or seeds } 

Base VariateStream, 
DepotVariateStream, 

Arrival VariateStream, 
MTTRVariateStream, 
NTRSVariateStream : RandomOb); 


MaxTime : REAL; 
SendTime > REAL; 
ReturnTime : REAL; 
tolerance : REAL; 


MaxReplications : INFEGER; 
CountReplications : INTEGER; 
TimeSample : INTEGER; 
Tinterval : INTEGER; 
NoOfRate : INTEGER; 
NoOfEqpt : INTEGER; 
NoOfSystem : INTEGER; 
NoOfStation : INTEGER; 
NumberDown — : INTEGER; 
MaxCanned : INTEGER; 


TimeChange, 
Rate : TimeArray; 


END MODULE. 


{ Implementation module for Globals } 
IMPLEMENTATION MODULE Globals; 
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FROM JOMod IMPORT FileUseType(Output,Input); 
{ for file handling } 


FROM StatMod IMPORT SINTEGER; 
{ for statistics collection } 


FROM SimMod IMPORT SimTime; 


FROM RandMod IMPORT FetchSeed; 
{ to obtain random seeds } 


FROM Egqpt IMPORT EquipmentOb), IdentityNo; 


FROM Report IMPORT ReportObj, CountNoOfNew, CountNoOfDispose; 


PROCEDURE UTILIZATION(IN Tchange, 
Trate : TimeArray; 
IN NoOfRate : INTEGER): REAL; 


VAR 
util : REAL; 
count : INTEGER; 
BEGIN 
count := 0; 
LOOP 
INC(count); 
IF SimTime() < Tchange[count] 
util := Trate[count]; 
RETURN util; 
END IF; 
END LOOP; 
END PROCEDURE; 


{seseccocececeereeetetesecceeeeneerccnececerteemmneneceeetetcoaoo} 


PROCEDURE Demand(IN identity : INTEGER; 
IN rate : REAL; 
IN system : INTEGER) : REAL; 


VAR 


demand : REAL; 
BEGIN 
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demand: = ItemArrmean[identity] * (1.0 / ( FLOAT(NoOfltem[identity]) 
* FLOAT(system) * rate)); 

RETURN demand; 

END PROCEDURE; 


PROCEDURE READINPUT, 


VAR 
Report : ReportOb); 
Equipment : EquipmentObj; 
seed], seed2, seed3, seed4, seedS SINLEGER: 
1, Icount, j, count, NSystem, [Level, number : INTEGER; 
DemandTime, DemandMean, Urate : REAL; 
stringdump : STRING; 

(saeccecccseeecceceanneeccnnansseeeseccnaseeeeecccaaaneeseccaaasaess 

BEGIN 


NEW(SystemFile); 
ASK SystemFile TO Open("SYS.out", Output); 


NEW(Ssystemifile ); 
ASK Ssystemfile TO Open("SSYS.out", Output); 


NEW (InputFile); 
ASK InputFile TO Open(“wsm.dat", Input); 


ASK InputFile TO ReadInt(seed1); 
ASK InputFile TO ReadLine(stringdump); 


NEW(Base VariateStream); 
ASK Base VariateStream TO SetSeed(FetchSeed(seed1)); 


ASK InputFile TO ReadInt(seed4); 
ASK InputFile TO ReadLine(stringdump); 


NEW(MTTRVanateStream); 
ASK MTTRVariateStream TO SetSeed(FetchSeed(seed4)); 
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ASK InputFile TO ReadInt(seedS); 
ASK InputFile TO ReadLine(stringdump); 


NEW(DepotVariateStream); 
ASK DepotVariateStream TO SetSeed(FetchSeed(seed5)); 


ASK InputFile TO ReadInt(seed2); 

ASK InputFile TO ReadLine(stringdump); 

NEW(Arrival VariateStream); 

ASK Arrival VariateStream TO SetSeed(FetchSeed(seed2)); 


ASK InputFile TO ReadInt(seed3); 

ASK InputFile TO ReadLine(stringdump); 
NEW(NTRSVariateStream); 

ASK NTRSVariateStream TO SetSeed(FetchSeed(seed3)); 


ASK InputFile TO ReadInt(MaxReplications); 
ASK InputFile TO ReadLine(stringdump); 
ASK InputFile TO ReadInt(TimeSample); 
ASK InputFile TO ReadLine(strngdump); 
ASK InputFile TO ReadInt(NoOfStation); 
ASK InputFile TO ReadLine(stringdump); 
ASK InputFile TO ReadInt(NoOfSystem); 
ASK InputFile TO ReadLine(stringdump); 
ASK InputFile TO ReadReal(SendTime); 
ASK InputFile TO ReadLine(stringdump); 
ASK InputFile TO ReadReal(ReturnTime); 
ASK InputFile TO ReadLine(stringdump); 
ASK InputFile TO ReadInt(NoOfRate); 
ASK InputFile TO ReadLine(stringdump); 
NEW(TimeChange, 1..NoOfRate); 
NEW(Rate, 1..NoOfRate); 


ASK InputFile TO ReadLine(stringdump); 
count := 0; 

WHILE count < NoOfRate 

INC(count); 


ASK InputFile TO ReadReal(TimeChange[count]); 
ASK InputFile TO ReadReal(Rate[count]); 
END WHILE; 


{ the time for the last rate is the maximum simulation time } 
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MaxTime := TimeChange[count]; 


ASK InputFile TO ReadInt(MaxCanned); 
ASK InputFile TO ReadLine(stringdump); 
ASK InputFile TO ReadReal(tolerance); 
ASK InputFile TO ReadLine(stringdump); 
ASK InputFile TO ReadInt(NoOfEqpt); 
ASK InputFile TO ReadLine(stringdump); 


ASK Sinventfile TO Open("SINV.OUT", Output); 
ASK Sbasefile TO Open("SBAS.OUT", Output); 


{ Create arrays to tract inventory levels and # in each pipelines } 


NEW(NoOfltem, 1..NoOfEqpt); 
NEW(ItemName, 1..NoOfEqpt); 
NEW(ItemArrmean, 1..NoOfEqpt); 
NEW(ItemDptmean, 1..NoOfEqpt); 
NEW(ItemBasmean, 1..NoOfEqpt); 
NEW(ItemReptime, 1..NoOfEqpt); 
NEW(ItemRemtime, 1..NoOfEqpt); 
NEW(ItemNtrs, 1..NoOfEqpt); 
NEW(InvLevel, 1..NoOfEqpt); 
NEW(Inventory, 1..NoOfEqpt); 
NEW(BasePipeNo, 1..NoOfEqpt); 
NEW(DepotPipeNo, 1..NoOfEqpt); 


ASK InputFile TO ReadLine(stringdump); 
Icount := 1; 
WHILE Icount < = NoOfEqpt 


ASK InputFile TO ReadString(ItemName[Icount}); 
ASK InputFile TO ReadReal(ItemArrmean{Icount}); 
ASK InputFile TO ReadReal(ItemNtrs[Icount]); 
ASK InputFile TO ReadInt(NoOfltem[Icount)); 

ASK InputFile TO ReadInt(InvLevel[Icount]); 

ASK InputFile TO ReadReal(ItemBasmean[Icount]); 
ASK InputFile TO ReadReal(ItemDptmean{[Icount}); 
ASK InputFile TO ReadReal(ItemRemtime[Icount]); 
ASK InputFile TO ReadReal(ItemReptime[Icount]); 


INC(Icount); 
END WHILE; 
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{ Initialize the Global records for cannibalization } 
ASK InputFile TO Close; 


END PROCEDURE; 


{ Definition Module for Repair } 
DEFINITION MODULE Repair; 
iy PE 

StationOby = OBJECT 


ASK METHOD ObjInit; 
{ method to initialize the StationObj object } 


ASK METHOD FinishRepair; 
{ method to ask a station which just finished repair to engage 
the next waiting equipmentObj for service } 


END OBJECT; 
END MODULE. 
{HSSSSSSSSSSSTSSeSSSSeesSeTSsseseeSeSTsesereessseSeses 
{== SS SSSSSeSSeeSeesSeeSeesessSeseseSeseSeeeereseses 


{ Implementation module for Repair } 
IMPLEMENTATION MODULE Repair; 


FROM Egqpt IMPORT EquipmentObj; 


FROM Globals IMPORT  IdleStationQueue, 
ServiceQueue; 


OBJECT StationOb); 
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ASK METHOD ObjInit; 
BEGIN 
END METHOD; 
ASK METHOD FinishRepair; 
VAR 
Equipment : EquipmentObj; 
number Waiting : INTEGER; 
BEGIN 
number Waiting := ASK ServiceQueue numberIn; 
IF numberWaiting = 0 
ASK IdleStationQueue TO Add(SELF); 
EESE 
Equipment := ASK ServiceQueue TO Remove(); 
TELL Equipment TO StartService(SELF); 
END IE; 
END METHOD; 


END OBJECT; 
END MODULE. 


{ Definition module for Stock } 
DEFINITION MODULE Stock; 


FROM GrpMod IMPORT QueueObj; 


TYRE 
InventoryQueueObj = OBJECT(QueueObj) 


ASK METHOD Objlnit: 
{ initialize InventoryQueueObj } 


ASK METHOD CheckInventory(IN id : INTEGER) ; 
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{ Constantly check for available so that systems can be revived } 
END OBJECT; 


InvitemOb) = OBJECT 
Level : INTEGER; 
id : INTEGER; 
DownNo |: INTEGER; 


ASK METHOD ObjInit; 
{initialize InvIltemObj } 


TELL METHOD Cannibalization; 
{ look for downed systems to cannibalize } 


ASK METHOD TakelInventory; 
{ uSe spare to revive system if available } 


END OBJECT; 
VAR 
Downnumber, identNo : INTEGER; 


END MODULE. 


{ Implementation module for Stock } 
IMPLEMENTATION MODULE Stock; 


FROM Globals IMPORT __ Inventory, 
InventoryQueue, 
System, 
ItemReptime, 
MaxCanned, 
SystemDown; 


{ 
OBJECT InventoryQueueObj; 
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ASK METHOD Objlnit; 
BEGIN 
END METHOD; 


ASK METHOD CheckInventory(IN id : INTEGER) ; 
VAR 
Invitem : InvItemOb); 
numberWaiting : INTEGER; 
BEGIN 
IF Inventory[id] > 0 
IF (ASK SELF numberIn) > 0 


Invitem := ASK SELF TO Remove(); 
ASK InvItem TO TakeInventory; 


END IF; 
END IF; 
END METHOD; 


END OBJECT; 


OBJECT InvItemObj; 


ASK METHOD Objinit: 
BEGIN 


DownNo := Downnumber; 
id := identNo; 


END METHOD; 
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TELL METHOD Cannibalization; 
{ uSe existing inventory to recover system as fast as possible. 


only when there is no spares left, then start to cannibalize. 
However cannibalization takes some time, use constant } 


VAR 
Down, numberzero : INTEGER; 


BEGIN 
Level := Inventory[id]; 


IF Level > 0 


{ since there is spare, can use it to revive system } 


{ time to replace faulty part } 
WAIT DURATION ItemReptime{id] 


ON INTERRUPT 
END WAIT; 


ASK SELF TO TakelInventory; 


ELSE 
{ No spares available Cannibalization Policy START } 


IF DownNo > MaxCanned 


{ Since cannibalization not allowed, must wait for spare to be available } 


ASK InventoryQueue[id] TO Add(SELF); 
TERMINATE; 
ELSE 


{ start cannibalizing on the downed equipment } 


Down := 1; 
IF DownNo > 1 
numberzero := 0; 


WHILE Down <= DownNo - 1 
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IF SystemDown[Down, id] > 0 
SystemDown{Down, id] := SystemDown[Down, id] - 1; 


{ Since able to cannibalize, let the current downed system be up after 
a constant time } 


{ 
WAIT DURATION 3.0 


ON INTERRUPT 
END WAIT; 
{ downed System is now available after cannibalization. } 


INC(System); 
DISPOSE(SELF); 


{ get out of while loop } 
TERMINATE; 


EESE 


{ Count no of downed systems which cannot be cannibalization } 
numberzero := numberzero + 1; 


END IF; 


Down := Down + 1; 
END WHILE; 


IF numberzero = DownNo - 1 
{ Let one more system be unavailable for cannibalization } 
SystemDown[DownNo, id] := 
SystemDown[DownNo, id] - 1; 

DISPOSE(SELF); 
TERMINATE; 

END IF; 

ELSIF DownNo = 1 


IF SystemDown{1, id] > 0 
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SystemDown[1, id] := SystemDown[]1, id] - 1; 
{ Since able to cannibalize, let the current downed system be up after 
a constant time } 
WAIT DURATION 3.0 
ON INTERRUPT 
END WAIT; 


{ downed System is now available after cannibalization } 
INC(System); 


DISPOSE(SELF); 
TERMINATE; 
ELSE 
ASK InventoryQueue[id] TO Add(SELF); 


TERMINATE; 
END IF; 


END IF; 


END IF; 
END IF; 


END METHOD; 


ASK METHOD Takelnventory; 
BEGIN 
{ take inventory now } 
DEC(Inventory[id]); 
{ now can make system available } 
INC(System); 
ASK InventoryQueue[id] TO CheckInventory(id); 
DISPOSE(SELF); 
END METHOD; 


END OBJECT; 
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END MODULE. 


{ Definition Module for Eqpt } 
DEFINITION MODULE Egpt; 


FROM Repair IMPORT StationObj; 


TYPE 
EquipmentObj = OBJECT 
Station : StationOb); 
identity : INTEGER; 
startTime : REAL; 


ASK METHOD Objinit; 
{ method to initialize the EquipmentObj object } 


TELL METHOD EnterSystem; 

{ method to sample for a new equipment to fail and then to 
direct the current failed equipment for appropriate repair 
actions } 


TELL METHOD CheckSystem; 

{ method to direct the failed equipment to check for available 
spare or to cannibalize a good part from downed part. To do 
this a new object InvItem is declared for this purpose. Invltem 
has its own methods for cannibalization. EnterSystem also tell 
the equipment to Chooserepair simultaneous } 


TELL METHOD ChooseRepair; 

{ ChooseRepair decides by uniform random sampling whether an 
equipment goes to base or depot for repair. 
For base repair, it will check for available base repair 
Station and then tell it to StartService } 


TELL METHOD StartService(IN Station: StationOb)); 
{ StartService samples a repair time to complete the repair 
of a failed equipment } 


END OBJECT; 
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VAR 
IdentityNo : INTEGER; 


END MODULE. 


{ Implementation module for Eqpt } 
IMPLEMENTATION MODULE Egpt; 


FROM SimMod IMPORT SimTime; 
{ SimTime shows the current simulation time } 


FROM UtilMod IMPORT ExitToOS; 
{ ExitToOS terminates the program when a condition Is satisfied } 


FROM Repair IMPORT StationOb); 
FROM Report IMPORT CountNoOfNew, CountNoOfDispose; 
FROM Stock IMPORT InvItemObj, identNo, Downnumber; 


FROM Globals IMPORT  InputFile, 
DepotVariateStream, 
Arrival VariateStream, 
NTRS VariateStream, 
MTTRVanriateStream, 
Base VariateStream, 
MaxTime, 

SendTime, 
ReturnTime, 
TimeChange, 
Rate, 
NoOfRate, 
NoOfEapt, 
NoOfltem, 
ItemName, 
ItemArrmean, 
ItemDptmean, 
ItemBasmean, 
ItemRemtime, 
ItemReptime, 
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ItemNtrs, 
NoOfSystem, 
System, 
NoOfStation, 
SystemDown, 
MaxCanned, 
Inventory, 
InvLevel, 
BasePipeNo, 
DepotPipeNo, 
ServiceQueue, 
IdleStationQueue, 
InventoryQueue, 
UTILIZATION, 
Demand; 


ASK METHOD Objlnit; 


BEGIN 
identity := IdentityNo; 
END METHOD; 


TELL METHOD EnterSystem; 
VAR 
Equipment2: EquipmentObj; 
interarrivalTime : REAL; 
DemandMean : REAL; 
Urate : REAL; 
nsys, error : INTEGER; 
BEGIN 


startTime := SimTime(); 


{ Start another object of the same type } 
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Urate := UTILIZATION(TimeChange, Rate, NoOfRate); 
nsys := System; 


{ get out if divide by 0 } 
IF nsys = 0 
ExitToOS(error); 
END IF; 
DemandMean := Demand(identity, Urate, nsys); 
interarrivalTime := ASK Arrival VariateStream Exponential(DemandMean); 
IF ( SimTime() + interarrivalTime ) < MaxTime 
IdentityNo := identity ; 
NEW(Equipmentz2); 
INC(CountNoOfNew); 
{ Schedule the next equipment of the same type to arrive } 
TELL Equipment2 TO EnterSystem IN interarrivalTime; 
END IF; 
TELL SELF TO CheckSystem; 


END METHOD; 


TELL METHOD CheckSystem; 


VAR 
Station : StationObj; 
Invitem : InvitemOb); 
I, numberldle, nsys, ninv, 
numberzero, Down : INTEGER; 
depottime : REAL; 
basetime : REAL; 
ratio : REAL; 
mttremove : REAL; 
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BEGIN 
{ add one to the Basepipeline } 


INC(BasePipeNo[identity]); 


{ wait for diagnosis and removal of failed part in by exponential assumption } 


mttremove := ASK MTTRVariateStream Exponential(ItemRemtime[identity]); 
WAIT DURATION mttremove 


{ downtime experience by a System } 
DEC(System); 

{ Check the current no of downed systems } 
Downnumber := NoOfSystem - System; 


{ cannibalization policy to revive system } 
identNo := identity; 
NEW(InvItem); 
TELL InvItem TO Cannibalization; 


{ send the faulty item for repair now } 
TELL SELF TO ChooseRepair; 


END METHOD; 


TELL METHOD ChooseRepair; 


VAR 
Station : StationObj; 
I, numberldle, nsys, ninv : INTEGER; 
depottime : REAL; 
basetime : REAL; 
ratio : REAL; 
mttremove : REAL; 


BEGIN 
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{ decision point for base/depot repair } 
ratio := ASK NTRSVariateStream UniformReal( 0.0, 1.0 ); 
IF ratio > ItemNtrs[identity] 
{ Base Repair starts } 
numberldle := ASK IdleStationQueue numberIn; 


IF numberldle = 0 
{ Enqueue the customer } 
ASK ServiceQueue TO Add(SELF); 


E ESE 
{ The customer will start service immediately } 
Station := ASK IdleStationQueue TO Remove(); 


TELL SELF TO StartService(Station); 
END IF; 


ELSE 


{ transportation time constant } 
WAIT DURATION SendTime; 
ON INTERRUPT 
END WAIT; 


{ add one to depotpipeline } 
INC(DepotPipeNofidentity]); 

{ repairtime depends on ItemDptmeanfidentity] } 

depottime := ASK DepotVariateStream 
Exponential(ItemDptmean{identity]); 
WAIT DURATION depottime; 
ON INTERRUPT 

END WAIT; 


{ minus one to depot pipeline } 


DEC(DepotPipeNofidentity]); 
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{ transportation time constant } 
WAIT DURATION ReturnTime; 
ON INTERRUPT 
END WAIT; 


{ basepipeline minus one since the part is sent to depot for repair 
completed repair } 


DEC(BasePipeNofidentity]); 
{ downed item is now available for Inventory[identity] } 
INC(Inventory[identity]); 


{ dispose the equipment to clear memory } 
DISPOSE(SELF); 


INC(CountNoOfDispose); 
END IF; 


END METHOD; 


TELL METHOD StartService(IN Station: StationObj); 


VAR 
ServiceTime : REAL; 


BEGIN 
ServiceTime := ASK Base VariateStream 
Exponential(ItemBasmean|identity]); 
WAIT DURATION ServiceTime 
ON INTERRUPT 
END WAIT; 
ASK Station TO FinishRepair; 


{ downed item is now available for Inventory[identity] } 


INC(Inventory[identity]); 
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{ minus one to the Basepipeline } 
DEC(BasePipeNolidentity]); 
DISPOSE(SELF); 
INC(CountNoOfDispose); 
END METHOD; 
{ 
END OBJECT; 


END MODULE. 


{ Definition module for Report } 
DEFINITION MODULE Report; 
TYPE 

ReportObj = OBJECT 


ASK METHOD Objinit; 
{ initialize ReportObj } 


ASK METHOD ReportNow; 
{ asking for statistical collections at specified intervals } 


ASK METHOD FinalReport; 
{ asking for final statistical results with confidence intervals } 


TELL METHOD ReportStatusNow; 
{ asking for statistical collections at specified intervals } 


END OBJECT; 
VAR 


CountNoOfNew : INTEGER; 
CountNoOfDispose : INTEGER; 
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END MODULE. 


{ Implementation module for Report } 
IMPLEMENTATION MODULE Report; 


FROM Globals IMPORT MaxTime, System, Inventory, InvLevel, 
SystemFile, 
Sinventfile, Ssystemfile, Sbasefile, 
Tinterval, TimeSample, 
SystemDown, NoOfSystem, MaxCanned, 
BasePipeNo, DepotPipeNo, 
CountReplications, 
ItemName, 
SystemStat, 
BaseDemandStat, 
InventoryStat, 
NoOfEapt; 


FROM SimMod IMPORT SimTime; 
FROM MathMod IMPORT SORT; 
FROM StatMod IMPORT IStatObj; 


OBJECT ReportObj; 


ASK METHOD Objlnit: 


VAR 
i,j) : INTEGER; 


BEGIN 
{ create new Statistical variables to collect weighted statistics } 
NEW(SystemStat, 1..TimeSample + 1); 


NEW(BaseDemandStat, 1..TimeSample+1, 1..NoOfEqpt); 
NEW(InventoryStat, 1..TimeSample +1, 1..NoOfEqpt); 
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FOR 1:= 1 TO TimeSample + 1 


ASK (GETMONITOR (SystemStat[i], Istatob))) 
TO SetHistogram(0,NoOfSystem, 1); 


FOR j := 1 TO NoOfEaqpt 
ASK (GETMONITOR(BaseDemandStat[1,}],Istatob)) ) 
TO SetHistogram(0,100,1); 
ASK (GETMONITOR (InventoryStat[i,j],[statobj)) 
TO SetHistogram(0, InvLevellj] ,1); 
END FOR; 


END FOR; 
END METHOD; 


ASK METHOD ReportNow; 
VAR 
identity : INTEGER: 
BEGIN 
ASK SystemFile TO WriteReal(SimTime(), 9,2); 
ASK SystemFile TO WriteInt(System, 6); 
ASK SystemFile TO WriteLn; 


SystemStat[Tinterval] := System; 


FOR identity := 1 TO NoOfEaqpt 


IF BasePipeNo[identity] - InvLevelfidentity] > 0 
InventoryStat[Tinterval,identity] := 
BasePipeNofidentity] - InvLevel[identity]; 
BLSE 
InventoryStat[Tinterval,identity] := 0; 
END IF; 
BaseDemandStat[Tinterval,identity] := BasePipeNofidentity] ; 


END FOR; 
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TELL SELF TO ReportStatusNow; 


END METHOD; 


TELL METHOD ReportStatusNow; 
BEGIN 
IF SimTime() < MaxTime 
{ wait a predetermined interval to call report again } 


WAIT DURATION (MaxTime/FLOAT(TimeSample)); 
END WAIT; 
Tinterval := Tinterval + 1; 


ASK SELF TO ReportNow; 
END IF; 


END METHOD; 


ASK METHOD FinalReport; 


VAR 
i, j : INTEGER; 
sysmean, syslowerlimit, sysupperlimit, ssdev, tolerance : REAL; 
basemean, baselowerlimit, baseupperlimit, bsdev, tol : REAL; 
invmean, invlowerlimit, invupperlimit, isdev : REAL; 


BEGIN 
ASK Ssystemfile TO WriteString("FINAL STATISTICS FOR SYSTEM 
AVAILABILITY, 

REPLICATIONS = "); 


ASK Ssystemfile TO WriteInt(CountReplications, 5); 

ASK Ssystemfile TO WriteLn; 

ASK Ssystemfile TO WriteString(" Sim.Time SysMean SysLLimit SysULimit"); 
ASK Ssystemfile TO WriteLn; 
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ASK Sbasefile TO WriteString("FINAL STATISTICS FOR BASE DEMAND 
NUMBERS, 
REPLICATIONS = "); 
ASK Sbasefile TO WriteInt(CountReplications, 5); 
ASK Sbasefile TO WriteLn; 


ASK Sbasefile TO WriteString(" Sim.Time"); 
FOR } := 1 TO NoOfEgqpt 


ASK Sbasefile TO WriteString(" ITEM = a F 

ASK Sbasefile TO WriteString(ItemName|j]); 

ASK Sbasefile TO WriteString(" oe 
END FOR; 


ASK Sbasefile TO WriteLn; 


ASK Sbasefile TO WriteString(" we 
FOR j := 1 TO NoOfEaqpt 
ASK Sbasefile TO WriteString(" Mean LLimit ULimit |"); 
END FOR; 
ASK Sbasefile TO WriteLn; 
ASK Sbasefile TO WriteLn; 


ASK Sinventfile TO WriteString("FINAL STATISTICS FOR INVENTORY 
LEVELS, 
REPLICATIONS = "); 
ASK Sinventfile TO WriteInt(CountReplications, 5); 
ASK Sinventfile TO WriteLn; 


ASK Sinventfile TO WriteString(" Sim.Time”); 
FOR j := 1 TO NoOfEgpt 


ASK Sinventfile TO WriteString(" ITEM = »); 

ASK Sinventfile TO WriteString(ItemName[j]); 

ASK Sinventfile TO WriteString(" "); 
END FOR; 


ASK Sinventfile TO WriteLn; 


ASK Sinventfile TO WriteString(" "); 
FOR j := 1 TO NoOfEqpt 
ASK Sinventfile TO WriteString(" Mean Llimit Ulimit |"); 
END FOR; 
ASK Sinventfile TO WriteLn; 
ASK Sinventfile TO WriteLn; 


FOR 1 := 1 TO TimeSample + 1 
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ASK Ssystemfile TO WriteReal( FLOAT(i-1) 
* (MaxTime/FLOAT(TimeSample)), 9,2); 


ASK Sbasefile TO WriteReal( FLOAT(i-1) 
* (MaxTime/FLOAT(TimeSample)), 9,2); 


ASK Sinventfile TO WriteReal( FLOAT(i-1) 
* (MaxTime/FLOAT(TimeSample)), 9,2); 


sysmean := GETMONITOR(SystemStat[i],Istatobj).Mean(); 
ssdev := GETMONITOR(SystemStat(1],[statobj).StdDev(); 


syslowerlimit := sysmean - (2.0 * 
ssdev/SQRT(FLOAT(CountReplications))); 

sysupperlimit := sysmean + (2.0 * 
ssdev/SQRT(FLOAT(CountReplications))); 


ASK Ssystemfile TO WriteReal(sysmean, 9,3); 
ASK Ssystemfile TO WriteReal(syslowerlimit, 9,3); 
ASK Ssystemfile TO WriteReal(sysupperlimit, 9,3); 
ASK Ssystemfile TO WniteLn; 


FOR j := 1 TO NoOfEapt 


basemean := GETMONITOR(BaseDemandStat{[Ii,j],[statobj).Mean(); 
bsdev := GETMONITOR(BaseDemandStat{i,}],Istatobj).StdDev(); 
baselowerlimit := basemean - (2.0 * 
bsdev/SQRT(FLOAT(CountReplications))); 
baseupperlimit := basemean + (2.0 * 
bsdev/SQRT(FLOAT(CountReplications))); 


ASK Sbasefile TO WriteReal(basemean, 9,3); 
ASK Sbasefile TO WriteReal(baselowerlimit, 9,3); 
ASK Sbasefile TO WriteReal(baseupperlimit, 9,3); 
ASK Sbasefile TO WriteString(" |"); 


invmean := GETMONITOR(InventoryStat[i,j],Istatobj).Mean(); 
isdev := GETMONITOR(InventoryStat[ij],Istatobj).StdDevQ); 
invlowerlimit := invmean - (2.0 * 
isdev/SQRT(FLOAT(CountReplications))); 
invupperlimit := invmean + (2.0 * 
isdev/SQRT(FLOAT(CountReplications))); 
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ASK Sinventfile TO WriteReal(invmean, 9,3); 
ASK Sinventfile TO WriteReal(invlowerlimit, 9,3); 
ASK Sinventfile TO WriteReal(invupperlimit, 9,3); 
ASK Sinventfile TO WriteString(" |"); 

END FOR; 


ASK Sbasefile TO WriteLn; 
ASK Sinventfile TO WriteLn; 


END FOR; 
END METHOD; 


( 
END OBJECT; 
END MODULE. 
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{ Main Module for the WSM Simulation } 
MAIN MODULE WSM; 


FROM SimMod IMPORT StartSimulation, SimTime, ResetSimTime; 
FROM IOMod IMPORT FileUseType(Output,Input); 

FROM StatMod IMPORT Istatobj; 

FROM MathMod IMPORT SQRT; 

FROM Egqpt IMPORT EquipmentObj, IdentityNo; 

FROM Report IMPORT ReportObj, CountNoOfNew, CountNoOfDispose; 
FROM Stock IMPORT InvItemOb); 

FROM Repair IMPORT StationOb); 


FROM Globals IMPORT 
Ssystemfile, 


MTTRVariateStream, 


ah 


Base VariateStream, 
DepotVariateStream, 


Arrival VariateStream, 
NTRS VariateStream, 


MaxReplications, 
CountReplications, 
MaxTime, 
TimeSample, 
Tinterval, 
SendTime, 
ReturnTime, 
TimeChange, 
Rate, 

NoOfRate, 
System, 
SystemDown, 
NoOfEgqpt, 
NoOfltem, 
ItemName, 
ItemArrmean, 
ItemDptmean, 
ItemBasmean, 
ItemReptime, 
ItemRemtime, 
ItemNtrs, 
NoOfSystem, 
MaxCanned, 
tolerance, 
Inventory, 
NoOfStation, 
InvLevel, 
BasePipeNo, 
DepotPipeNo, 
InventoryQueue, 
IdleStationQueue, 
ServiceQueue, 
SystemStat, 
BaseDemandsStat, 
InventoryStat, 
UTILIZATION, 
READINPUT, 
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Demand; 


VAR 
Station : StationOb); 
Invitem : InvItemOb); 
Report : ReportOb); 
Equipment : EquipmentOb); 
i, Icount, j, count, Nsystem, Ilevel, number : INTEGER; 
DemandTime, DemandMean, Urate : REAL; 
sySmean, ssdev > REAL; 
basemean, bsdev, tol : REAL; 
invmean, isdev : REAL; 
countconvergence : INTEGER; 
{ seseececcccccssnsseteeeetteecennnnneneeeeeeeeecnnnnnnnnnananaaaaaaas 
BEGIN 


{ procedure to read in datafile } 
READINPUT; 


NEW(InventoryQueue, 1..NoOfEqpt); 
NEW(SystemDown, 1..MaxCanned, 1..NoOfEqpt); 
NEW(Report); 
{ do replications for statistics collection } 
CountReplications := 1; 


{ loop for replications until statistical convergence is reached } 
LOOP 


ResetSimTime(0.0); 
Tinterval := 1; 


System := (NoOfSystem); 


NEW(ServiceQueue); 
NEW(IdleStationQueue); 
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FOR i := 1 TO NoOfStation 
NEW(Station); 
ASK IdleStationQueue TO Add(Station); 
END FOR; 


FOR Icount : 


1 TO NoOfEgqpt 


IdentityNo := Icount; 
NEW(Equipment); 


INC(CountNoOfNew); 
{ initialize the pipelines } 


BasePipeNo[Icount] := 0; 
DepotPipeNo[Icount] := 0; 


NEW(InventoryQueue[Icount]); 
Inventory[Icount] := (InvLevel[(Icount]); 


Urate := UTILIZATION(TimeChange, Rate, NoOfRate); 
DemandMean := Demand(Icount, Urate, NoOfSystem); 
DemandTime := ASK Arrival VariateStream Exponential(DemandMean); 


TELL Equipment TO EnterSystem IN DemandTime; 
END FOR; 


{ need to initialize the NoOfltem[Icount] for each item } 


FOR 1 := 1 TO MaxCanned 
FOR j} := 1 TO NoOfEapt 
SystemDown{[l, J] := NoOfltemf[j]; 
END FOR; 
END FOR; 


ASK Report TO ReportNow; 


{ ask MODSIM to start simulation } 
StartSimulation; 


{ dispose of all remaining objects to clear memory } 


number := ASK IdleStationQueue numberIn; 
IF number <> 0 
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FOR 1 := 1 TO number 
Station := ASK IdleStationQueue TO Remove(); 
DISPOSE(Station); 
END FOR; 
END IF: 
DISPOSE(IdleStationQueue); 


number := ASK ServiceQueue numberlIn; 
IF number <> 0 
FOR 1:= 1 TO number 
Equipment := ASK ServiceQueue TO Remove(); 
DISPOSE(Equipment); 
INC(CountNoOfDispose); 
END FOR; 
END IF; 
DISPOSE(ServiceQueue); 


{ Dispose InvItem that are not processed } 
FOR i := 1 TO NoOfEapt 
number := ASK InventoryQueue[i] numberIn; 


IF number <> 0 
FOR j := 1 TO number 
Invitem := ASK InventoryQueue[i}] TO Remove(); 
IF InvItem <> NILOBJ 
DISPOSE(InvItem); 
END IF; 
END FOR; 
END IF; 
DISPOSE(InventoryQueue[i]); 
E FOR; 


ASK Ssystemfile TO WriteLn; 

ASK Ssystemfile TO WriteString("# NEW = "); 

ASK Ssystemfile TO WriteInt(CountNoOfNew, 5); 
ASK Ssystemfile TO WriteLn; 

ASK Ssystemfile TO WriteString("# DISPOSED = "); 
ASK Ssystemfile TO WriteInt(CountNoOfDispose, 5); 
ASK Ssystemfile TO WriteLn; 


OUTPUT("REPLICATION NUMBER = ", CountReplications); 


A 
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{ check for statistical convergence } 


IF CountReplications > 5 
countconvergence := 0; 
FOR i := 1 TO TimeSample + 1 


sysmean := GETMONITOR(SystemStat[i],Istatobj).Mean(); 
ssdev := GETMONITOR(SystemStat[i],Istatobj).StdDevQ); 


tol [= (2.0 * ssdev/SQRT(FLOAT(CountReplications))); 
IF (tol <= sysmean * tolerance ) 

INC(countconvergence) 
END IF; 


FOR j := 1 TO NoOfEqpt 


basemean := GETMONITOR(BaseDemandStat{1,j],Istatobj).Mean(); 
bsdev := GETMONITOR(BaseDemandStat{i,j],[statobj).StdDev(); 
tol [= (2.0 * bsdev/SQRT(FLOAT(CountReplications))); 


IF (tol < = basemean * tolerance ) 
INC(countconvergence) 
END IF; 


invmean := GETMONITOR (InventoryStat[i,j],Istatob)).Mean(); 
isdev := GETMONITOR(InventoryStat[i,j],Istatobj).StdDevQ); 
tol [= (2.0 * isdev/SQRT(FLOAT(CountReplications))); 


IF (tol <= invmean * tolerance ) 
INC(countconvergence) 

END IF; 

END FOR; 
END FOR; 
OUTPUT("CONVERGENCE NUMBER = ", countconvergence); 
IF countconvergence >= (TimeSample+1 + 

((TimeSample + 1)*NoOfEqpt*2) ) 
EXIT 


END IF; 
END IF; 


144 


IF CountReplications > MaxReplications 
EXD 
END LE: 
INC(CountReplications); 


{ replications end } 
END LOOP; 


{ Report Final Statistics Now } 

ASK Report TO FinalReport; 
OUTPUT("ENDED NORMALLY"); 

END MODULE. 


{ END OF MODSIM PROGRAM } 
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APPENDIX B. MATHLAB PROGRAM FOR ANALYTICAL WSM. 


tt tees See EES SS Te REE Eee Ke eee ee eee Te eee eee ee ie eee 


TANALYTICAL WAR SUSTAINABILITY MODEL (WSM) PROGRAM 

Jo WRITTEN IN PC-MATLAB BY LIM HUNG HENG FOR THE 

% COMPLETION OF M.SC (OA) THESIS 

Jo AT NAVAL POSTGRADUATE SCHOOL 

ot tt Ses eee e eee eee eee ee eee EE ee ee eee FT ee cia ai olla 
%oThis MATLAB program used the assumptions and analytical solutions based 
Yoon Chapter III - Analytical Model Development and the wartime scenario 
ZYogiven in Chapter V - A Case Study. 

Too ToTo To To Yo Yo To Yo Vo To To Vo To Yo To Ve Po Yo Vo Vo Yo To Yo Yo Yo Vo Yo Yo Vo To Yo Vo Yo Vo Vo Vo 
JDATA INPUT PARAMETERS: 


Opyaccinciiccocncacccesesncakaceigiaren tone eeee eae ene = 
% itemfr = Failure rate (failures / 1,000,000 hours) of each LRU 

Y% itemnrts § = The proportion of repair not repairable at the base 

Jo (nrts = Not-Repairable This Station) for each LRU. 

Jo itemqpa = The quantity of each LRU found in each system. 

Jo (qpa = quantity per application). 

Zo invievel = The initial inventory level stocked at the base for each LRU. 
% brtime = Base repair time for each LRU. 

% dptime = Depot repair time for each LRU. 

Jo totalsys = Total number of available systems at the start of war period. 
Zo nmcs = Maximum number of Not-Mission Capable Systems to be 

%o tolerated by the operational users (an operational target). 

% fortime = Transportation time from base to depot. 

% rettime = Transportation time from depot to base. 

Zo timchg = time at which utilization rate of the system changes. 

% urate = different levels of utilization rates of the system 

%o (to be determined by operational users). 

Yo notimint = number of time intervals to be analyzed. 

Tp os SOT ot vance ca Sea ee 

itemfr = (390 4040 360 11510 4490 6070 570 4480 17000 10920); 

itemnrts = [.07 .06 .21 .59 .04 .04 .24 .09 .06 .28};: 

itemgqpa = {111111111 1); 

invievel = {13152412 62]; 

brtime = [48 48 72 72 48 48 72 48 48 48]; 

dptime = (120 120 120 120 120 120 120 120 120 120}; 
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totalsys = 24; 

nmcs = ae 
fortime = 120; 
rettime = 120; 
timchg = [168 720}; 
urate = {0.2 0.1]; 
notimint = 100; 


JoT0T0 To To To Vo Yo To Yo Yo Vo Yo Yo Yo Yo Yo Yo Yo Yo Vo YoVo Yo Vo Yo Vo Vo To Vo Vo Vo Ve Yo Yo Vo Vo Vo 
YoCalculate equi-distance length between time intervals. 


timeint = 0:(timchg(2)/notimint):timchg(2); 


Y%omake itemnrts into matrix with time intervals as rows and LRUs as columns. 
itemnrts = ones(timeint) * itemnrts; 


%Calculate the demand rate for each LRU at each time intervals. 
ZYobased on the given utilization rates. 


demratl = (1.e-6 * totalsys * urate(1)) * (itemfr); 
demratl = ones(timeint)’ * demrat1; 

demrat2 = (1.e-6 * totalsys * urate(2)) * (itemfr); 
demrat2 = ones(timeint)’ * demrat2; 

SNE 2S. os na dceadeeaaddacnecnaeakwne 


JoF0V%0 0 Yo To Yo To To Yo Vo Vo Yo Yo Yo Vo Yo Vo Vo Yo Yo Vo Vo Vo Vo Yo Yo Yo To Vo Vo Vo Vo Vo Vo Vo Vo Vo 
Yobase repair service time for each LRU at each time intervals. 
basetime = ones(timeint)’ * brtime; 


Yotimedim is just a matrix of time intervals used for this computation. 
timedim = timeint’ * ones(brtime); 

Yocompute the various base service times at which the LRU is still held up. 
JYoservice time before reaching the basetime. 

Ybastim11 depends only on timedim. 


i = (timedim <= basetime); 
bastim11= (i) .* timedim; 


Yoservice time after the basetime but before the change from high 
JYoutilization to low utilization. 
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Yobastim12 has a constant time. 
i = (timedim > basetime) & (timedim < = timchg(1)); 
bastim12 = (1) .* basetime; 


Yoservice time during low utilization but before (timchg(1) + basetime) 
Zosince we know that LRUs that failed under the high-utilization demand 
ZYotime timchg(1) must wait for the whole basetime to be repaired. 


1 = (timedim > timchg(1)) & (timedim <= (timchg(1) + basetime) ); 
bastim13 = (i) .* ((timchg(1) + basetime) - timedim); 
bastim21 = (i) .* (timedim - timchg(1)); 


ZYoservice time after (basetime + timchg(1)) 
Yobastim22 depends only on basetime. 

1 = (timedim > (timchg(1) + basetime)); 
bastim22 = (i) .* basetime; 


Yobastim1 is the matrix addition to be used together with demrat1 the 
Yohigh-utilization demand rate 

Yobastim2 is the matrix addition to be used together with demrat2 the 
Yolow-utilization demand rate 

bastim1 = bastim11 + bastim12 + bastim13; 

bastim2 = bastim21 + bastim22; 

basdemand1 = ((1 - itemnrts) .* demrat1) .* bastim1]; 

basdemand2 = ((1 - itemnrts) .* demrat2) .* bastim2; 


Yobasdemand is the mean demand caused by the depot turnaround time. 
basdemand = basdemand1 + basdemand2; 


Zo To Fo Yo To Vo Yo To To Yo Yo Yo Yo Yo Yo Yo Yo Vo Fo Yo Yo Yo To To To Yo Yo Vo To Vo To To Vo Fo Yo Yo Vo Vo 
Yodepot repair service time for each LRU at each time intervals. 


Yosince dptime is constant, the depot turnaround time from base to 
Yodepot is given below. 

depotime = fortime + rettime + dptime; 

depotime = ones(timeint)’ * depotime; 

Yocompute the possible depot demand for high utilization 


Yocompute the various depot service times at which the LRU is still held up. 
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Joservice time before the change from high utilization to low utilization. 
Yodptim11 depends only on timdim. 

1 = (timedim <= timchg(1)); 

dpttim11= (1) .* timedim; 


JYoservice time during low utilization but before the depot turnaround time. 
Yodptim12 has a constant time and dpttim21 depends on (timedim - timchg(1)). 
i = (timedim > timchg(1)) & (timedim < = depotime ); 

dpttim12 = (i) .* timchg(1); 

dpttim21 = (1) .* (timedim - timchg(1)); 


Yoservice time between the depot turnaround time and (depotime + timchg(1)) 
Zosince we know that LRUs that failed under the high-utilization demand 
Yotime timchg(1) must wait for the whole depottime to be repaired. 

i = (timedim > depotime) & (timedim <= (timchg(1) + depotime) ); 
dpttim13 = (i) .* ((timchg(1)+depotime) - timedim); 

dpttim22 = (i) .* (timedim - timchg(1)); 


Joservice time after (depotime + timchg(1)) 
Yodptim23 depends only on depotime. 

1 = (timedim > (timchg(1) + depotime) ); 
dpttim23 = (1) .* depotime; 


Yodptim1 is the matrix addition to be used together with demrat1 the 
%high-utilization demand rate 

Yodptim2 is the matrix addition to be used together with demrat2 the 
JYolow-utilization demand rate 

dpttim1 = dpttim11 + dpttim12 + dpttim13; 

dpttim2 = dpttim21 + dpttim22 + dpttim23; 

dptdemand1 = (itemnrts .* demrat1) .* dpttim]; 

dptdemand2 = (itemnrts .* demrat2) .* dpttim2; 





Beni ice ve Rs, 2 So a se 7 


Yodptdemand is the mean demand caused by the depot turnaround time. 
dptdemand = dptdemand1 + dptdemand2 ; 


ToV0 Vo To Vo Vo Vo Yo Vo Vo Vo Vo Vo Yo Vo Yo Vo Vo Yo Vo Yo Vo To Vo Yo Vo Yo Vo Po Vo Yo Vo Vo Vo Yo Vo Vo Vo 


Yodemand is the Total Average Pipeline in according to Equation 27 of 
Yoof Chapter III. It is used to compute the SMOEs and MMOEs. 
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demand = basdemand + dptdemand; 


Too Vo To Vo To Vo Vo To Vo Yo Vo Yo Vo Yo Yo Vo Yo Yo Vo Yo To Vo Yo To YoYo Yo Yo Vo Pe To Yo Vo Pe Vo Le Ve 
%Calculate the Expected Number Of Backorders Ebi(t) for each LRU. 


JYomake a matrix out of invlevel based each time intervals. 
invmat = ones(timeint)’ * invlevel; 


% if demand is more than invmat, then there will be backorders 
backord1 = demand - invmat; 
ZYoindex to recognize those elements with backorders. 


b1 = find(backord1 > 0.0); 
invl = invmat(b]l); 


dem1 = demand(b1); 
back1 = backord1(b1); 
Op iwswatenceecankncnasdateete sree soe sine see ee ee 


Yebegin poisson computation by looping from zero to the maximum inventory 
ZYolevel m. 


m = max(invl); 


ZY this is the poisson term when k=0, or when the actual demand k 1s zero. 
ZY poiterm is the number of backorders 

Zosee Chapter 3, Section la for explanation why start from k = 0 

poiterm = (invl) ; 

prob = poiterm; 


for 1 = 1l:m-1 
ll = 1 * ones(dem1); 


Yothe poi(k,x) is a function to compute the poisson probability that 
Yothe actual number is k when the mean number is x. k and x can be 
Yin any form (scalar, vector or matrix). 

probi = poi(ii,dem1); 

YoAt each incremental loop, compute only those which inventory level is 
Yothan or equal to i. 

}) = (invl >= 1); 


poit = j .* ((inv] -1i).* probi); 


Yesum up all the poisson probabilities according to equation 3.35 
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poiterm = poiterm + poit; 
prob = prob + probi; 
end 


ZoToVo To Vo To Vo Vo Po Vo Ve Vo Vo Po Vo Po Vo Lo Po Vo Vo Vo Vo Vo Vo Vo Vo Po Vo To Vo Vo Vo Vo Vo Vo Vo Ve 
% the computation of the expected backorder for each LRU. 

% (see equation 3.35) 

backorder = 0.0 * backordl; 

backorder(b1) = (backl + (poiterm .* exp(-dem1)) ); 


ToToVo Vo Vo To Fe Vo Poo Vo Vo VoVo Vo Vo Vo Vo Vo Vo Vo Vo Vo Yo Vo Vo Vo Vo Vo Vo Vo Vo Vo Vo Vo Vo Vo Ve 
% For Non-cannibalization policy 


Yocalculate Total Expected Backorders (see equation 3.37). 
totaleb = zeros(timeint)’; 

for i = 1:max(size(itemfr)) 

totaleb = totaleb + backorder‘(:,1); 

end 


Too VoTo To To Vo Po Po Vo Yo Vo Vo Yo Vo Yo Vo Vo Yo Vo Yo Lo Yo YoT0 YoYo Vo Vo Vo Vo Vo Vo Vo Vo Vo VoVo 
% For Non-cannibalization policy 


YoCalculate operational availability (see equation 3.39). 
av = 1 - (backorder/totalsys); 
ao = ones(timeint)’; 


for 1 = 1:max(size(itemfr)) 
ao = ao .* av(:,1); 
end 


ooo oe So Ln chee ewwneewecescewes 
%Calculate ENMCS without cannibalization 
J%o(see equation 3.48). 





enmcs = totalsys * (1 - ao); 


JoTo To ToVo Vo Yo Yo Vo Yo Vo Vo Yo To Vo Vo Yo Yo Vo Vo YoYo Vo Yo Vo Yo Yo Yo Vo Vo Vo Vo Yo Vo Vo Yo Vo Vo 
% For Cannibalization Policy 


%Probability of shortages less than or equal to j 
%(see Equations 3.43 and 3.44 of Chapter III). 


Yoprshort is P(j): Probability of having less than or equal to 
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% } NMCS. Form unity matrix for prshort from 
prshort = ones(timeint)’ * ones(1:totalsys + 1); 


for n = O:totalsys 

%if demand is more than invmat, then there will be backorders 
backord3 = demand - invmat+n; 

ZYoindex to recognize those elements with backorders. 


b3 = find(backord3 > 0.0); 
inv3 = invmat(b3); 
dem3 = demand(b3); 


ml = max(inv3 +n); 


Yo start with poisson at k=0 as actual number and demand as mean number. 
poij0 = ones(dem3); 


fori = 1:ml 

ii = 1 * ones(dem3); 

Youse the poisson function (matrix form) 
poic = poi(ii,dem3); 

J) = (inv3+n >= i); 

poic = j.* poic; 

poyO = poi0 + poic; 

end 


probnbo = ones(timeint)’ * ones(itemfr); 
probnbo(b3) = poij0 .* exp(-dem3); 


Yoprshort(:,1) can be interpreted as the probability of all 
Josystems availability 


for j = 1:max(size(itemfr)) 
prshort(:,n+1) = prshort(:,n+1) .* probnbo(:,j); 
end 


Yocalculate Prob(ENMCS <= NMCS) (same result for non-cannibalization 


[sy 


Yoor cannibalization 
pnmesc = prshort(:,nmcs); 


TooToo To Yo To Yo Vo Yo Vo Yo To YoYo Yo Yo To Vo Vo Vo Yo Vo Vo Vo Vo Yo To Vo Vo Ve Yo Vo Vo Vo Vo Vo Ve 
Yocalculate expected number of NMCS (ENMCSc) for cannibalization 


enmesc = zeros(timeint)’; 

for} = 1:totalsys+ 1; 

enmcsc = enmcsc + (1 - prshort(:,j)) ; 
end 


JoTo To To To Vo Vo Vo Vo Lo Yo Vo Yo Yo Yo Yo Vo Vo Yo Vo Yo Yo Vo Yo Vo Vo Vo Vo Yo Vo Yo Vo Vo Yo Lo Vo VoVe 
Yocalculate operational availability with cannibalization 
aoc = 1 - (enmcsc/totalsys); 


JToT0 To Vo To Vo Yo Vo To Yo Yo Yo Vo Yo Vo Vo Yo Yo Yo Yo Yo Yo Yo Yo Yo Yo YoYo Vo Vo Vo Vo Vo Vo Vo Vo Vo Ve 


function b = poi(k,x) 

ZYopoisson function that handles matrix, vector or scalar. 

b = ((x.*k) ./ fact(k)); 

Jo Vo VoVo Vo Too Yo Vo Yo Yo Yo Yo Yo Yo Vo Yo Vo Yo Vo Yo Yo Vo Yo Fo To Yo Vo Yo Vo Vo Vo Vo Vo Vo Vo Vo Ve 


function b = fact(k); 
% factorial function that operates on matrix, vector or scalar. 


m = max(k(:)); 
b = ones(k); 
for1 = l:m 


n = 1 * ones(k); 

j = find(k >= 1); 

c = ones(k); 

cJ) = nQ); 

b=b.*c¢ 

end 

Too To To Yo To YoYo Fo Vo Yo Yo Vo Wo Yo Vo Yo Yo To YoYo Yo Vo Vo Yo Vo Vo Yo Vo Vo Vo Vo Vo Yo Vo Vo Yo Vo 
%END OF MATLAB PROGRAM 
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APPENDIX C. NUMERICAL RESULTS FROM ANALYTICAL VERSION OF THE WSM. 


Tables shown below are numerical results from the Analytical Case One. 


TABLE C-1 : ANALYTICAL RESULTS OF THE TOTAL PIPELINE FOR EACH IBU. 
[=e [ws | we | me 














6.2208000e-02 1.9889280e+00 























0.00 | 0.0000000e+00 | 0.0000000e+00 | 0.0000000e+00 0.0000000e+00 
36.00 | 6.7392000e—-02 6.9811200e-01 

72.00 9.3000960e-02 | 9.5874048e-01 | 1.2441600e-01 3 .9778560e+00 
108.00 9.7718400e—-02 1.0006272e+00 | 1.3747968e-01 9 .1513235e+00 
144.00 1.0243584e-01 | 1.0425139e+00 | 1.5054336e-01 6.3247910e+00 
180.00 9.5921280e-02 | 9.6804864e-01 | 1.5323904e-01 7.1667706e+00 
216.00 6.6942720e—-02 6.6087936e-01 | 1.3519872e-01 | 7.3457741e+00 
252.00 | 6.9301440e-02 | 6.8182272e-01 | 1.2534912e-01 7.6606877e+00 
288.00 | 7.1660160e-02 | 7.0276608e-01 | 1.3188096e-01 8 .2474214e+00 


1.3841280e-01 8.8341552e+00 


324.00 7.4018880e-02 | 7.2370944e-01 























| 360.00 7.6377600e-02 | 7.4465280e-01 | 1.4494464e-01 9.4208890e+00 
| 396.00 7.4018880e-02 | 7.2370944e-01 | 1.3841280e-01 | 8.8341552e+00 
| 432.00 7.1660160e-02 | 7.0276608e-01 | 1.3188096e-01 8.2474214e+00 
| 468.00 6.9301440e-02 | 6.8182272e-01 | 1.2534912e-01 7 .6606877e+00 
| 504.00 6.6942720e-02 | 6.6087936e-01 1.1881728e-01 7.0739539e+00 
| 540.00 6.5370240e-02 | 6.4691712e-01 | 1.1446272e-01 6.6827981e+00 
| 576.00 6.5370240e-02 | 6.4691712e-01 | 1.1446272e-01 6.6827981e+00 

612.00 6.5370240e-02 | 6.4691712e-01 | 1.1446272e-01 | 6.6827981e+00 

648.00 6.5370240e-02 | 6.4691712e-01 | 1.1446272e-01 6.6827981le+00 
- 684.00 6.5370240e-02 | 6.4691712e-01 | 1.1446272e-01 | 6.6827981e+00 | 
| 720.00 


6.5370240e-02 | 6.4691712e-01 | 1.1446272e—-0 6.6827981le+00 


a 
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MIABEe C-1 cont..) 


- 0000000e+00 
- 7587200e-01 
-0551859e+00 
-0862208e+00 
-1172557e+00 
-0189786e+00 
-6207744e-01 
. 7759488e-01 
-9311232e-01 
.0862976e-01 
-2414720e-01 
.0862976e-01 
-9311232e-01 
- 7759488e-01 
-6207744e-01 
.5173248e-01 
-5173248e-01 
-5173248e-01 
-5173248e-01 
-5173248e-01 
-5173248e-01 


0 
a 
1 
1 
1 
a: 
6 
6 
6 
7 
7 
7 
6 
6 
6 
6 
6 
6 
6 
6 
6 


0 .0000000e+00 
7.7414400e-01 
1.0786406e+00 
1.1483136e+00 
1.2179866e+00 
1.1586355e+00 
8.4123648e-01 
8.7607296e-01 
9.1090944e-01 
9.4574592e-01 
9.8058240e-01 
9.4574592e-01 
9.1090944e-01 
8.7607296e-01 
8.4123648e-01 
8.1801216e-01 
8.1801216e-01 
8.1801216e-01 
8.1801216e-01 
8.1801216e-01 
8.1801216e-01 


0. O0O00000e+00 
1.0488960e+00 
1.4264986e+00 
1.4684544e+00 
1.5104102e+00 
1.3775501e+00 
8.9505792e-01 
9.1603584e-01 
9.3701376e-01 
9.5799168e-01 
9.7896960e-01 
9.5799168e-01 
9.3701376e-01 
9.1603584e-01 
8.9505792e-01 
8.8107264e-01 
8.8107264e-01 
8.8107264e-01 
8.8107264e-01 
8.8107264e-01 
8.8107264e-01 


0.0000000e+00 
2.9376000e+00 
4 .0343040e+00 
4.2105600e+00 
4 .3868160e+00 
4.0734720e+00 
2.7809280e+00 
2 .8690560e+00 
2.9571840e+00 


| 3.0453120e+00 


3.1334400e+00 
3 .0453120e+00 
2 .9571840e+00 
2 .8690560e+00 
2. 7809280e+00 
2.7221760e+00 
2. 7221760e+00 
2.7221760e+00 
2.7221760e+00 
2.7221760e+00 
2.7221760e+00 
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0. 0O000000e+00 
9.8496000e—02 
1.9699200e-01 
2.2063104e-01 
2-4427008e-01 
2.5149312e-01 
2.2588416e-01 
2.1275136e-01 
2.2457088e-01 
2.3639040e-01 
2.4820992e-01 
2.3639040e-01 
2.2457088e-01 
2.1275136e-01 
2.0093184e-01 
1.9305216e-01 
1.9305216e-01 
1.9305216e-01 
1.9305216e-01 
1.9305216e-01 
1.9305216e-01 





0.0000000e+00 
1.8869760e+00 
2 .8682035e+00 
3.3965568e+00 
3 .9249101e+00 
4.1387674e+00 
3.7236326e+00 
3 .9878093e+00 
4.2519859e+00 
4.5161626e+00 
4.7803392e+00 
4.5161626e+00 
4.2519859e+00 
3 .9878093e+00 
3.7236326e+00 
3.5475149e+00 
3.5475149e+00 
3.5475149e+00 
3.5475149e+00 
3.5475149e+00 
3 .5475149e+00 


TABLE C-2 : ANALYTICAL RESULTS OF BACKORDERS FOR EACH IRUJ 


Time 


IkU A 


0.0000000e+00 
0.0000000e+00 
0. 0000000e+00 
0. 0000000e+00 
0 .0000000e+00 
0. 0000000e+00 
0.0000000e+00 
0.0000000e+00 
0.0000000e+00 
0. 0000000e+00 
0. 0000000e+00 
0.0000000e+00 
0. 0000000e+00 
0. 0000000e+00 
0. 0000000e+00 
0.0000000e+00 
0. 0000000e+00 
0.0000000e+00 
0.0000000e+00 
0. O0O00000e+00 
0. 0000000e+00 


0.0000000e+00 
0.0000000eE+00 
0. 0000000e+00 
9.6399166e—-01 
1.7570721e+00 
2 -4317595e+00 
2-5837480e+00 
2-8571034e+00 
3 .3835732e+00 
3.9275290e+00 
4.4843017e+00 
3.9275290e+00 
3 .3835732e+00 
2-8571034e+00 
2.3540075e+00 
2.0350727e+00 
2 .0350727e+00 
2.0350727e+00 
2.0350727e+00 
2.0350727e+00 
2-0350727e+00 





0. 0000000e+00 


0.0000000e+00 | 


0.0000000e+00 


| 
| 


- 0.0000000e+00 | 
- 0.0000000e+00 


0. 0000000e+00 


0.0000000e+00 | 


0.0000000e+00 
0. 0000000e+00 


0.0000000e+00 | 


0. 0000000e+00 
0.0000000e+00 
0.0000000e+00 
0.0000000e+00 
0. 0O000000e+00 
0.0000000e+00 
0.0000000e+00 
0.0000000e+00 
0.0000000e+00 
0.0000000e+00 
0.0000000e+00 


0.0000000e+00 
0.0000000e+00 
0.0000000e+00 
0.0000000e+00 
0.0000000eE+00 
0.0000000e+00 
0.0000000e+00 
0.0000000e+00 
0.0000000e+00 
0.0000000e+00 
0.0000000e+00 
0.0000000e+00 
0.0000000e+00 
0.0000000e+00 
0.0000000e+00 
0.0000000e+00 
0.0000000e+00 
0.0000000e+00 
0.0000000e+00 
0.0000000e+00 
0.0000000e+00 
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kU C 


0.0000000e+00 
0.0000000e+00 
0.0000000e+00 
0.0000000e+00 
0.0000000e+00 
0.0000000e+00 
0.0000000e+00 
0.0000000e+00 
0.0000000e+00 
0.0000000e+00 
0.0000000e+00 
0.0000000e+00 
0.0000000e+00 
0.0000000e+00 
0.0000000e+00 
0.0000000e+00 
0.0000000e+00 
0.0000000e+00 
0.0000000e+00 
0.0000000e+00 
0.0000000e+00 


0.0000000e+00 
0.0000000e+00 
0.0000000e+00 
0.0000000e+00 
0. 0000000e+00 
0.0000000e+00 
0.0000000e+00 
0.0000000e+00 
0.0000000e+00 
0.0000000e+00 
0.0000000e+00 
0.0000000e+00 
0.0000000e+00 
0.0000000e+00 
0.0000000e+00 
0.0000000e+00 
0.0000000e+00 
0.0000000e+00 
0.0000000e+00 
0.0000000e+00 
0.0000000e+00 





{ ee 


ieee. C=2 cont. .) 


0. 0000000e+00 
0.0000000e+00 
0. 0000000e+00 
0.0000000e+00 
0 .0000000e+00 
0.0000000e+00 
0.0000000e+00 
0. 0000000e+00 
0.0000000e+00 
0 .0000000e+00 
0.0000000e+00 
0.0000000e+00 
0. Q0O00000e+00 
0. 0000000e+00 
0.0000000e+00 
0.0000000e+00 
0.0000000e+00 
0.0000000e+00 
0. 0000000e+00 
0. 0000000e+00 
0.0000000e+00 


oO OQ O00 00 00 0 O00 00 0000 00 © 


. O0O00000e+00 
. 0O000000e+00 
. 0000000e+00 
. O0O0O0000e+00 
. O0O0O0000e+00 
. O0O00000e+00 
. 0000000e+00 
. O0O0O0000e+00 
. O000000e+00 
. O0O00000e+00 
. 0000000e+00 
-. O0O00000e+00 
. O0O00000eE+00 
- 0000000e+00 
. 0000000e+00 
- O000000eE+00 
- 0000000e+00 
- 0000000e+00 
- 0000000e+00 
- O000000e+00 
- 0000000e+00 


0.0000000e+00 
0.0000000e+00 
0.0000000e+00 
0.0000000eE+00 
0. O0O0O0000e+00 
0. 0000000e+00 
0.0000000e+00 
0. 0O0O00000e+00 
0. 0000000e+00 
0.0000000e+00 
0.0000000e+00 
0.0000000e+00 
0. 0000000e+00 
0.0000000e+00 
0.0000000e+00 
0.0000000e+00 
0.0000000e+00 
0.0000000e+00 
0. 0000000e+00 
0.0000000e+00 
0.0000000e+00 





BESS 


- O0O00000e+00 
- O0O00000e+00 
» 1447217e+00 
-59772787e+00 
-0418910e+00 
-2366346e+00 
-8618360e+00 
-0988250e+00 
-3409888e+00 
-5873900e+00 
-8372476e+00 
-5873900e+00 
. 3409888e+00 
-0988250e+00 
-8618360e+00 
. 7072617e+00 
1.7072617e+00 
1.7072617e+00 
1.7072617e+00 
1.7072617e+00 
1.7072617e+00 
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APPENDIX D. NUMERICAL RESULTS FROM SIMULATION VERSION OF THE Woh. 
Tables below show the numerical results for the case of Simulation 


Variation Four. 





TABLE D-1 : SIMULATION RESULTS OF THE TOTAL PIPELINE 
FOR EACH LRU WITH STATISTICAL CONFIDENTIAL LIMITS 


IRU A at 95% C.I. IRV B at 95% C.I. 
Sim.Time Mean Lhimit ULimit Mean LLimit ULimit 


- 000 0.000 0.000 





PRPRPRP RPP PrPOOO 





OO 0000 0 
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(TABLE D-1 cont..) 


324.00 
360.00 


720.00 





. 000 
sO>3 
093 
- 126 
- 158 
7 
226 
0292 
-2/0 
~274 
2/3 


SS Sere © OOo) OOO © O° S70.08 


. 000 
. 046 
084 
-116 
147 
- 184 
213 
238 
296 
- 260 
259 
ofl 
236 
219 
203 
186 
0.181 
0.167 
0.157 
0.155 
0.143 


© ©O-@ @ @ 2 O71 @ @ Oo ' © 24 e7® 


IRJ C at 95% C.I. 


0.000 
0.059 
o- 102 
0.136 
0.169 
0.209 
0.239 
weZoo 
0.284 
0.287 
0.287 
0.279 
0.262 
0.245 
0.228 
Oi9211 
0.204 
0.190 
0.179 
0.177 
0.165 


IRU D at 95% C.I. 


0.000 
> LO 
3.448 
4.967 
6.344 
72424 
7.952 
8.216 
8.194 
7.960 
7.503 
6.991 
6.449 
6.020 
5.687 
5eoc9 
5.234 
5.141 
5.093 
5.072 
5.074 


0.000 
02 
3.396 
4.904 
6.274 
7.350 
7.879 
8.143 
SB. 122 
7.889 
72434 
6.922 
6.381 
2953 
5.622 
5.324 
5.169 
5.076 
5.029 
5.007 
5.009 





0.000 
1.848 
3.500 
2.029 
6.413 
7.498 
8.026 
8.288 
8.266 
8.031 
Yeol2 
7.060 
Geol 
6.086 
5.153 
5.454 
2.0299 
5.205 
3.158 
Delo 7 
55.139 


—— IRUJ E at 95% C.I. IRD F at 95% C.I. 


0.00 
36.00 
72.00 

108.00 
144.00 
; 180.00 
216.00 
252.00 
288.00 
324.00 
360.00 
396.00 
432.00 
468 .00 
504.00 
540.00 
576.00 
612.00 
648.00 
684.00 
720.00 


| SS eee ee 


0.000 
0.555 
0.924 
1.246 
1.525 
1.740 
1.770 
1.715 
1.562 
12593 
1.209 
1.027 
0.897 
0.776 
0.694 
0.625 
0.600 
O59 1 


| 0.587 


0.583 
0.560 


0.000 
0.534 
0.896 
1.214 
1.489 
1.703 
1.733 
1.680 
1.529 
1.360 
ee 
0.997 
0.869 
0.749 
0.668 
0.601 
0.577 
0.568 
0.565 
0.561 
0.538 


0.000 
0.576 
0.952 
1.278 
1.560 
1.778 
1.806 
1.751 
1.596 
1.426 
1.240 
1.057 
0.925 
0.803 
0.719 
0.649 
0.624 
0.615 
0.610 
0.605 
0.581 


0.000 
0.786 
1.288 
1.693 
2.046 
2.315 
2.334 
2.198 
1.983 
e729 
1.512 
1.301 
1.146 
1.011 
0.927 
0.860 
0.811 
0.798 
0.785 
0.760 
0.742 


0.000 
0.760 
1.254 
1.655 
2.004 
2e2tk 
2.290 
2.156 
1.942 
1.690 
1.476 
1.266 
1.113 
0.980 
0.897 
0.831 
0.784 
0.771 
0.759 
0.734 
0.716 


0.000 
0.811 
Pec) 
e732 
2.089 
2.360 
2.377 
2.240 
2.023 
1.767 
1.549 
1.336 
1.180 
1.043 
0.957 
0.889 
0.839 
0.825 
0.811 
0.786 
0.768 





(TABLE D-1 cont..) 





- 000 
085 
~ 154 
-216 
- 263 
314 
358 
-400 
~420 
-421 
-409 
soo9 
3/71 
349 
320 
298 
282 
- 268 
- 260 


oooooo0oo0o0o 00°00 0 902 0 00 0000 


0.000 
2.204 
3.579 
4.766 
5.759 
6.358 
6.046 
5.509 
4.991 
4.497 
4.006 
3.561 
3.239 


lie@oo71 


2.760 
2.585 
2.470 
2.428 
2.355 
2.343 
2.347 


OOO O00 00 00 00 000 00 00 0 
r= 
© 
ee 


0.000 
Zeloz 
3.521 
4.697 
5.683 
6.279 
5.968 
5.434 
4.918 
4.426 
3.937 
3.496 
3.177 
2.911 
2. /03 
Pa EO 
2.417 
e577 
2.306 
2.295 
2-279 


0.000 
0.093 
0.166 
0.230 
0.277 
0.329 
0.374 
0.417 
0.437 
0.439 
0.426 
0.416 
0.387 
0.364 
02535 
0.313 
0.296 
0.282 
0.274 


0.000 
2.246 
3.636 
4.835 
5.834 
6.436 
6.124 
5.585 
5.065 
4.568 
4.075 
3.626 
3.301 
3.031 
2.817 
2.639 
2.522 
2.479 
2-404 
2.391 
2.395 
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0.000 
1.447 
2-621 
3.630 
4.523 
5.199 
5.342 
5.252 
5.049 
4.710 
4.344 
3.976 
3.625 
3.340 
3.131 
2.964 
2.849 
2.778 
2.740 
2.J2/ 


| alae IRU G at 95% C.I. IRU H at 95% C.I. 


0.000 
1.516 
2./14 
3.740 
4.644 
5.324 
5.465 
5 
5 
g 






IRU A at 95% C.I. 
ULimit 
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TABLE D-2 : SIMULATION RESULTS OF THE BACKORDERS FOR EACH 
LAO WITH STATISTICAL CONFIDENCE LIMITS 


IRV B at 95% C.I. 


Mean 


0.000 
0.001 
0.013 
0.041 
0.074 
0.113 
0.115 
0.103 
0.083 
0.061 
0.042 
0.029 
0.021 
0.015 
0.012 
0.009 
0.006 
0.005 
0.004 
0.006 
0.004 


LLimit 


0.000 
0.000 
0.009 
0.034 
0.065 
0.100 
0.103 
0.091 
0.073 
0.053 
0.035 
0.023 
0.016 
0.011 
0.009 
0.006 
0.003 
0.003 
0.002 
0.003 


0.002 


TES r= 


ULimit 








(TABLE "D=2 ‘Gont a. ) 


-———_. 
| 


| 0.00 
36.00 
72-00 
108.00 
144.00 
180.00 
216.00 
252.00 
288.00 
324.00 
360.00 
396.00 
432.00 
468.00 
504.00 
540.00 
576.00 
612.00 
648.00 
684.00 
720.00 


| 0.00 
36.00 
72.00 

108.00 
144.00 
180.00 
216.00 
252.00 
288.00 
324.00 
360.00 
396.00 
_ 432.00 
| 468.00 
504.00 
540.00 
576.00 
612.00 
648.00 
684.00 
720.00 





IRU C at 95% C.I. 


- 000 
- 002 


0.000 
0.000 
0.003 
0.005 
0.005 
0.010 
0.012 
0.014 
0.017 
0.017 
0.017 
0.017 
0.015 
0.011 
0.010 
0.009 
0.007 
0.005 
0.004 
0.005 
0.005 


0.000 
0.003 
0.006 
0.010 
O.011 
0.017 
0.020 
0.022 
0.026 
0.025 
0.025 
0.026 
0.023 
0.018 
0.016 
0.015 
0.013 
0.009 
0.009 
0.010 
0.010 


IRU E at 95% C.I. 


0.000 


0.027 


0.000 
0.018 
0.081 
Ou171 
0.286 
0.380 
0.372 
0-325 
0.262 
0.213 
0.155 
0.116 
0.088 
0.065 
0.048 
0.038 
0.033 
0.032 
0.030 
0.024 
0.021 


0.000 
0.028 
0.102 
0.202 
0.326 
0.426 
0.416 
0.377 
0.299 
0.246 
0.184 
0.142 
0.110 
0.085 
0.066 
0.053 
0.047 
0.046 
0.043 
0.036 
0.033 


IRV D at 95% C.I. 
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0.000 
0.015 
0.233 
0.852 
Lee: 
2.645 
3.084 
3.316 
3.289 
3.074 
2.662 
2-253 
1.814 
1.503 
1.283 
1.108 
1.026 
0.963 
0.938 
0.940 
0.941 


0.000 
0.011 
0.213 
0.813 
1.716 
2.580 
3.016 
3.248 
3.222 
3.008 
2.600 
2.174 
1.759 
1.452 
122356 
1.064 
0.984 
0.921 
0.896 
0.899 
0.899 


0.000 
0.020 
0.253 
0.892 
1.826 
2.709 
3.5m 
3.384 
3.356 
Sel oe 
Pa PAS, 
2ec92 
1.869 
1.554 
1.330 
15153 
1.069 
1.005 
0; SE 
0.982 
0.982 


IRJ F at 95% C.I. 


0.000 
0.001 
0.017 
0.053 
0.103 
0.156 
0.142 
0.112 
0.088 
0.056 
0.038 
0.025 
0.020 
0.013 
0.014 
0.010 
0.007 
0.005 
0.003 
0.003 
0.003 


0.000 
0.000 
0.013 
0.044 
0.090 
0.140 
0.128 
0.099 
0.077 
0.048 
0.031 
0.019 
0.015 
0.009 
0.010 
0.006 
0.004 
0.002 
0.001 
0.002 
0.001 


0.000 
0.000 
0.021 


0.061 


0.115 


O.771 
0.157 
O7125 
0.099 
0.065 
0.044 
0.030 
0.025 
0.017 
0.019 
0.014 
0.010 
0.007 
0.005 
0.005 
0.005 


| 
| 
| 
| 


(TABLE D-2 cont..) 


0.00 
36.00 
72.00 
108.00 
144.00 
180.00 
216.00 
z52.00 
288.00 
324.00 
360.00 
396.00 
432.00 
468.00 
504.00 
540.00 
576.00 
612.00 
648.00 
684.00 
720.00 


i] altOwO0 
36.00 
72.00 

108.00 
144.00 
180.00 
216.00 


288 .00 
324.00 
360.00 
396.00 
432.00 
468 .00 
504.00 
540.00 
576.00 
612.00 
648.00 
684.00 
720.00 


oooooo0ooooeoo0ooo 00 00 00 


IRJ G at 95% C.I. 


000 
003 
~014 
~024 
~031 
044 
049 
057 
-061 
-061 
~055 
048 
-040 
032 


IRU I at 95% C.I. 


0.000 
0.011 
0.161 
0.502 
0.954 


| 1.287 


1.124 
0.848 
0.632 
0.467 
0.339 
0.235 
0.175 
0.138 
0.105 
0.084 
0.065 
0.061 
0.046 
0.038 
0.037 


0.000 
0.001 
0.011 
0.019 
0.026 
0.037 
0.042 
0.050 
0.053 
0.053 
0.047 
0.041 
0.034 
0.026 
0.021 
0.019 
0.016 
0.018 
0.018 
0.019 
0.018 


0.000 
0.007 
0.144 
0.470 
0.909 
1.233 
1.074 
0.804 
0.594 
0.434 
0.312 
0.212 
0.155 
0.121 
0.090 
0.070 
0.053 
0.050 
0.036 
0.029 
0.028 





0.000 
0.005 
0.018 
0.029 
0.036 
0.050 
0.056 
0.064 
0.069 
0.068 
0.062 
0.055 
0.046 
0.037 
0.030 
0.028 
0.025 
0.027 
0.027 
0.028 
0.027 


2° ee 


0.000 
0.015 
0.178 
0.535 
1.000 
1.340 
1.174 
0.891 
0.669 
0.499 
0.366 
0.257 
0.194 
0.156 
0.120 
0.097 
0.077 
0.073 
0.056 
0.046 
0.046 


LRU H at 95% C.I. 


0.000 
0.023 
0.110 
0.212 
0.332 
0.450 
0.474 
0.430 
0.376 
0.303 
0.236 
0.180 
0.127 
0.102 
0.081 
0.065 
0.060 
0.051 
0.043 
0.045 
0.040 


0.000 
0.019 
0.098 
0.195 
0.311 
0.426 
0.449 
0.406 
0.354 
0.283 
0.219 
0.166 
0.114 
0.091 
0.071 
0.056 
0.051 
0.044 
0.036 
0.038 
0.033 


Tine = of 
0.000 
0.028 
©. 12a 
0.228 
0.353 
0.475 
0.499 
0.454 
0.398 
0.323 
0.252 
0.195 
0.139 
0.113 
0.091 
0.074 
0.069 
0.059 
0.049 
0.052 
0.047 


IRU J at 95% C.I. 


0.000 


0.277 
0.999 
1.840 
2.655 
3.300 
3.431 
3.341 
3.138 
2.810 
2-464 
2.128 
1.813 
1.584 
1.411 
1.283 
1.194 
1.140 
1.104 
1.101 
1.111 


= 


0.000 
0.258 
0.962 
1.791 
2.998 
3.239 
a, 571 
3.282 
3.080 
2.754 
2.410 
2.076 
1.765 
1.537 
1.367 
1.241 
1.153 
1.099 
1.065 
1.062 





0.000 
0.296 
1.035 
1.890 
2.712 
3.360 
3.491 
3.400 
3.196 
2.866 
Zooks. 
2-179 
1.862 
1.632 
1.456 
1.325 
1.235 
1.180 
1.144 
1.140 
1.151 





APPENDIX E. OPUS-8 INPUT REQUIREMENTS AND RESULTS 
A. OPUS-8 INPUT FORMAT AND REQUIREMENTS 


*PROBLEM DEFINITION rev.2 
a 
'WRITE PROBLEM HEADER BELOW: | 
a ee 
! ANALYTICAL WSM USING OPUS-8 APPROXIMATION WITHOUT REPAIR \ 
Se 
i i er 

! PROBLEM TYPE (MARK YOUR CHOICE WITH AN "X"): ! 
a enema 
'X! ANALYSIS OF A GIVEN ALLOCATION OF SPARES | 
en tf 

! ! INITIAL OR REPLENISHMENT PROCUREMENT OF SPARES | 
rt ef 

! ! REALLOCATION OF A GIVEN ASSORIMENT OF SPARES | 
a re 

! ! REALLOCATION FOLLOWED BY REPLENISHMENT PROCUREMENT ! 

RR RR RE RSE SSRN St emi 
of ee ci Sn rr er mre een enn meemematane 

! INITIAL STOCK (MARK YOUR CHOICE WITH AN "X"): ! 
Pf 

! ! ZERO STOCK ! 
Ce nm ee a net Se nN EDS A nc emf 
'X! A GIVEN ALLOCATION OF SPARES : 
ee ee ee a a ee eee 

! ! AVERAGE NUMBER OF DEMANDS DURING RESUPPLY TIME (OR TIME T) ! 
I a 
a 
'MEASURE OF EFFECTIVENESS (MOE) (MARK YOUR CHOICE WITH AN "X")! 
a RR NR SN eh 

! STEADY STATE PROBLEMS: \ 
ee i ee i mene er meh 

! ! WAITING TIME (WTI) ! 
a RS RE 

! ! PROBABILITY OF MISSION SUCCESS (PMS) | 
CRS Reem Sr NE em meee re ean 

! ENDURANCE PROBLEMS: ! 
Ph tee 

! ! PROBABILITY OF NO BACKORDER DURING TIME PERIOD T (PNBO) ! 
RSE RM SES A SS NS need = 


! ! SUM OF BACKORDER TIME DURING TIME PERIOD T (SBT) ! 
SO a aan eS 
'X! NUMBER OF BACKORDERS DURING TIME PERIOD T (NBO) ! 


$4} 
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*CONTROL PARAMETERS rev .3 
eae cee arene f 


! LIMITS FOR INVESTMENT AND/OR MOE: <DEFAULT> ! 
: { MINIMAL INVESTMENT << ec Oa > ! 
! ! MAXIMAL INVESTMENT < 1E+18 > ! 
: ! MINIMAL MOE <MOE DEP.> ! 
: { MAXIMAL MOE <MDE DEP.> ! 


! ! NUMBER OF POINTS IN FINAL CURVE < 30 > ! 
! FORMAT FOR ALL OUTPUT: <DEFAULT> ! 
: ! NUMBER OF LINES PER PAGE < 64 > 
: ! NUMBER OF CHARACTERS PER LINE < 80 oe 
—-_—$— 
! USE OF THE COMPOUND POISSON DISTRIBUTION IN THE CALCULATIONS: ! 
! (MARK SELECTED MODE WITH "X") (VMR = VARIANCE / MEAN) 
{ ! AUTOMATIC COMPOUND POISSON TO HANDLE WAITING TIME VARIANCE “4 
!X! STANDARD POISSON DISTRIBUTION (VMR = 1 FOR DEMAND AT ALL STATIONS) ! 
! ! COMPOUND POISSON DEMANDS (VMR VALUE ENTERED EXPLICITLY BELOW) 1 
! ! VMR FOR DEMAND AT STATIONS DIRECTLY SERVING SYSTEMS ! 
en af 
! ! EXCLUDE INITIAL STOCK IN RESULTS ("X" IF YES) : 
SH, 
'MULTIPLE REMOVAL HANDLING (MARK YOUR CHOICE WITH AN "X"); ! 
! (ONLY RELEVANT IF REMOVAL RATE FACTOR EXCEEDS 1 FOR SOME ITEM)! 
! ! POSSIBLE MULTIPLE REMOVALS INCLUDED IN GIVEN REPAIR TIMES ! 
! ! POSSIBLE MULTIPLE REMOVALS WILL INCREASE GIVEN REPAIR TIMES! 
! ! CREATE INPUT FILE TO PROGRAM "OPUS-8W" Cee YES) © 
! ! FILE NAME” <"infile".W>! 
: ! OPUS 8W GROUP IDENTIFICATION ! 
pa a ree 
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*STATION DATA 
rT 


eee (2) 1! (3) !(4) 1 (5) 3 
! V ! QUAN- !SUB !GROUP ! 
'ID !! DENOMINATION ! TITY !SET! OF ! 
ITAG !! ! PER ! OF ! STA- ! 
1 !SUPPOR-!STA-! TIONS ! 
1 ! TING ! TION! ! 
! i 'sTa- ! ! 
1 'TION(S)!(*) ! ! 
at 
! DEFAUL?: Po avtete- |! 
fp pt 
!IOC1!!DEMAND GEN ST ! 1 |! ! 


pt 


*SYSTEM DATA 
en nr enc a ene fe 
! (2) 1 (2) 1! @) 1 3) ! 
! DEPLOYED!! QUAN- !UTILIZA-! MEAN! 
! SYSTEM ! aT !! TITY !TION PER! TIME |! 
1 | 
| ’ 


ID TAG’ !FOLLOWING!! ' CALENDAR! TO 
!STATION (S) ! ! TIME ! REPAIR 
oe 


: DEFAULT : LP = SO : 0.0 : 


{SYSTEM : LOC 24 ! 0.1467! 2.0 i 


*ITEM DATA 

rr rrr fener ememenef 
» (1) it (2) ! (3) : (4) ! (5) 
1 ! ! FAIIURE ! ITE |! 
! ID TAG !! DENOMINATION ! UNIT ! RATE  ! (REPAIR ! 
1! ! PRICE ! 10 TO-6 !AND STOCK)! 
! 1! ! ! ! CATEGORY ! 
ta nnn nn ne nnn nnn nnn rn nnn t 
! DEFAULT: !! ! - ! -  !<ITEM ID ! 
Hp a pt 
!LRUO1 1! oA ! 1000. ! 390. ! 
! LRUO2 1! OB ! 1000. ! 4040. ! 
!TRUO3 oe ! 1000. ! 360. ! 
!LRUO4 1! 6D ! 1000. !11510. ! ! 
! LRUOS lho £ ! 1000. ! 4490. ! ! 
! LRUO6 1! OF ! 1000. ! 6070. ! 
! LRUO7 eee ! 1000. ! 570. ! 
! LRUO8 1!) 6H ! 1000. ! 4480. ! 
!TRUO9 Ll eee ! 1000. ‘17000. "2 
!LRU10 3 ! 1000. ! 10920. ! 
a et 
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*TTEM STRUCTURE rev.2 


i ()) (2) tc (4 eee! 
| ! MOTHER !! !REMOVAL !ENVIRON-! 
! ITEM ! ITEM !!QUAN-! RATE ! MENT ! 
!IDTAG ! oR  !!TTTY ! FACTOR ! FACTOR ! 
! ! SYSTEM !! ! ! ! 
! (ID TAG) !! jee)  ! 
po mf fe np at 
! DEFAULT: ipl eed Oar et! 
nr nn tat 
!LRUO1 IsySTEM ss}! «12 fo (21.o ! ! 
!TRUO2 ISYSTEM ist}! «1! (21.0) ! ! 
!TRUO3 isyYSTEM ss}! 1} 2100) ! ! 
! TRUO4 ISYSTEM st}! «2 6! «21.60 ! ! 
!LRUOS ‘SYSTEM is}! 2 of 21.0) ! ! 
! LRUO6 ISYSTEM ss}! «2! «2100 ! 
! LRUO7 ISYSTEM st! 1 to (21.Oo ! 
! TRUO8 IsySTM os !! 12 6! «21.0 ! ! 
!LRUOS IsysSTEM st}! «2 fo (21.o ! ! 
!LRU1O ISYSTEM ss}! «12! «(2100) ! 


*STOCK POLICY AND TIME PERIOD 
cp 


! TO BE USED IN ENDURANCE PROBLEMS 


! (1) (Qe (ayeet 64) 
! !STOCKED AT!! ! ! 
! ITEM !FOLLOWING !! INITIAL! TIME ! 
! CATRGORY ! (GROUP OF)!! STOCK !PERIOD T ! 
! ! STATIONS !! LEVEL ! ! 
a el 
! DEFAULT: re 
qe fp en ee | eet 
!TRUO1 Tere] et! i! 360m! 
!LRUO2 ! focl !! 3 !360 ! 
!TRUO3 ! foc: !! 1 #!360. ! 
!LRUO4 ! jock !! 5 !360 ! 
!LRUOS ! focl !! 2 !360. ! 
!LRUO6 Peeper t! deee!’ 3608 2 
!LRUO7 ! jock !! 1 !360. ! 
!LRUO8 tr ioe fh 2 360. |! 
!LRUOS ! toc !! 6 !360. ! 
!LRU10 ioc tees! 360. sa 
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*ALLOCATION OF SPARES ee 


STOCK LEVELS PER ITEM AND STATION 
a 
! = (1) ! (2) ! (3) TOTAL NUMBER OF EACH STATION ' 
! TEM ! OLA ————_— rn a rr ee 
| ! NUM-! (4) STATION ID TAG ! 
! ID TAG ! BER! ! 
wp ap a a a a pe ep ft 


s/f s ffs tut YE 1 1 ft 1 2 
ey / 


LT TA 
//////// tor! 


ee ne ae ef et 


> LRUO1 
: LRUO2 
> LRUO3 
- LRUO4 


- LRUOG6 
> LRUO7 
- LRUO8 


t 

t 

{ 

I 

! TRU0S)—! 
j 

t 

( 

! LRUOS | 
t 


| t 
0 | 
( | 
1 5! 
aoe 
1 
| | 
| ! 
| ! 
t 


NANF LPN UH WW! 


LRU10 


*OUTPUT SELECTION rev.2 
a 


MARK WITH "'/" TO SUPPRESS OUTPUT (BLANK MEANS YES) : 


A cn 
!/! ID: INPUT DATA ! 
ma ee 


! ! ID.1 ‘INPUT FORMS ‘ 
! ! ID.2.1 INITIAL STOCK LEVELS ! 
! ! ID.2.2 DEPLOYMENT OF SYSTEMS : 
! ! ID.2.3 UTILIZATION OF SYSTEMS ! 
! ! ID.3 SIGNIFICANCE LEVEL DATA ! 


em me ea EE ER SE te 


! IR: INTERMEDIATE RESULTS : 


mm an NN RES SASS a mete = 


! ! IR.1 SYSTEM SUMMARY DATA ! 
a REESE A NN 
! ! IR.2 ITEM AND STATION RELATED VARIABLES : 
Sn a pe 
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! ! IR.2.1 DEMAND RATE 

IR.2.2 RESUPPLY TIME GIVEN NO SHORTAGE (OR TIME PERIOD T) 
IR.2.3 THE COMPONENTS OF RESUPPLY TIME 

! ! IR.2.4 AVERAGE NUMBER OF DEMANDS 


1 

1 

‘ 

1 

! 

! 

\ 

1 

! ! IR.2.5 INITIAL ALLOCATION OF SPARES 

| 

!/! IR.2.6 SIGNIFICANCE LEVELS ! 

ppp 

!7! IR.3 SYSTEM AND DGS RELATED VARIABLES ! 

ee 

! ! IR.3.1 MDT GIVEN NO SHORTAGE { 

St 

!/!TR.4 SIGNIFICANCE LEVEL DATA ! 

ee eR SE pc aE a 

! ! IR.4.1 NUMBER OF POINTS ! 

cf 

! ! IR.4.2 PRINTING CALCUIATED POINTS OF THE C/E-CURVE(S) ! 

ppm pana 

{ ! IR.4.2.1 ALL CALCULATED POINTS 

+—+ 

{ ! IR.4.2.2 ALL OPTIMAL POINTS ! 

+—+ 

! ! IR.4.2.3 ALL SELECTED POINTS ! 

a 

{ ! IR.4.3 PLOTTING CALCULATED POINTS OF THE C/E-CURVE(S) ! 

mes 
IR.4.3.1 ALL CALCULATED POINTS 


a 


I 


: 
! ! IR.4.3.2 ALL OPTIMAL POINTS : 


{ ! IR.4.3.3 ALL SELECTED POINTS : 

{ ! FR: FINAL RESULTS : 
! ! FR.1 RESULTS PER POINT (IN FOLLOWING INTERVAL IF NOT ALL POINTS)! 
: ! INVEST. LOWER LIMIT : 

: ! INVEST. UPPER LIMIT : 

'/! FR.1.1 TOTAL STOCK LISTING 

! ! FR.1.2 ITEM AND STATION RELATED VARIABLES : 

! ! FR.1.2.1 ALLOCATION OF SPARES : 

tat : 
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!/! FR.1.2.2 VALUE OF OPTIMIZATION MOE ! 
+—+ 
!/! FR.1.2.3 EXPECTED SUM OF BACKORDERTIME DURING TIME PERIOD T ! 
+—+ 
!/! FR.1.2.4 PROBABILITY OF BACKORDER DORING TIME PERIOD T ! 
+—+ 
! ! FR.1.2.5 EXPECTED NUMBER OF BACKORDERS (*) ! 
+—+ 
'/! FR.1.2.6 RISK OF SHORTAGE (*) ! 
+—+ 
!/! FR.1.2.7 TOTAL RESUPPLY TIME INCL. WAITING TIME (*) : 
+—+ 
!/! FR.1.2.8 INVESTMENT ! 


Pe et 
! ! FR.1.3 SYSTEM AND DGS RELATED VARIABLES 


FR.1.3.1 VALUE OF OPTIMIZATION MOE 


= 


, aia 


FR.1.3.2 OPERATIONAL AVAILABILITY (*) 


os 


FR.1.3.3 EXPECTED NUMBER OF NOR (*) 


FR.1.3.4 EXPECTED NUMBER OF BACKORDERS (*) 


Peo 


mn 


RISK OF SHORTAGE (*) 


IS] 


OV 


FR.1.3.6 TOTAL MEAN DOWN TIME INCL. WAITING TIME (*) 
pp ft 
! 


<< 


d 
d 
: 


! ! FR.1.5 MEAN VALUES PER STATION FOR EACH RELEVANT MEASURE 
! ! FR.2 PRINTING FINAL C/E-CURVE FOR INV. VS ALL RELEVANT MEASURES ! 
! ! FR.3 TRACEBACK OF FINAL POINTS TO INTERMEDIATE CURVES : 
! ! FR.4 MARGINAL COST EFFECTIVENESS : 


!/! FR.4.1 FOR INVESTMENT VERSUS WAITING TIME(*) ! 
+—+ 
!! FR.4.2 =" " " EXPRCTED NUMBERS OF BACKORDERS(*) ! 
+—+ 
' !'§FR.4.3 (8 " " NOR (*) $ 
ae 
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B. RESULTS USING OPUS-8'S ENDURANCE OPTION TO APPROXIMATE WSM 


OPUS 8 (rev 3.2 ) MINDEF Singapore 
ANALYTICAL WSM USING OPUS-8 APPROXIMATION WITHOUT REPAIR 


JAN 08 1991 15-45-29 


a 
! OPUS-8 version 3.2 1990-03-21 
! Copyright (C) Systecon AB 1990 ! 
! This copy of OPUS-8 is licensed for use by ! 
! MINDEF Singapore (Major Lim) ! 
! according to the terms and conditions of ! 
! "Right-of-use to Systecons proprietary software" ! 


nn 


meen Te RMEDIATE RESULTS 
| 


Number of different items: 10 
Number of different stations: 1 
Number of stock positions: 10 


Number of individual systems: 24 
Operating time per system: 360.0 
Total operating time: 8640. 


NOTE! The figures above do not consider utilization factors. 
All "Sum Of Backorder Time" values calculated by OPUS-8 must 
be compared against operating times as presented above! 


Investment of initial stock: 27000.0 


*SYSTEM SUMMARY DATA 


(1) Vey ol 3) Sh (4) OP (5), -! ey 4 
! '! TOTAL ! TOTAL ! ! ! ASYMP- | 
! SYSTEM !! FATIURE ! ITEM ! MTBF ! DEMAND ! TOTIC ! 
' IDTAG !! RATE ! PRICE ! ! FLOW ! AVAL ! 
] !1(10 TO-6)! ! 1(10 TO-6)! ABILITY ! 
op cf nc 


! SYSTEM !! 59830.0 ! 10000.0 ! 16.71 ! 8777.1 ! : 


fj pt 


*ALLOCATION OF SPARES 
pede - et 


STOCK LEVELS PER ITEM AND STATION 
a ee et 
! (1) ! (2) ! (3) TOTAL NUMBER OF EACH STATION ! 
TEM TOTAL +--+ ee! 
] ! NUM-! (4) STATION ID TAG ! 
| 1 DiiAe es BERS! ! 


SSS SSR ! ! ! ! 
ey, 


en af 
YJ SSSI S TOC | ee ee eee 
eet 
Pees cee ee ee i, 
ec cy ee) 2 SOR PU hg, 
1tego3 tt ea 0 ae 
1 TROOA $5 8 5 i st eee 
\tmos5  «! at 2 ee ee ee 
(TRUO6G 80) 4 ge ee 
et. 6cy A eel ee 
{ TRUOS 9! 2 2 es ee 
| TH09 = $6} OG I! se 
eae lL ee LE eC 
Oe a ep at 
*DEMAND RATE 

cm ee SSS etme Sm SSS ue ecm 
! TOTAL DEMAND RATE PER ITEM AND STATION (10 TO-6) ! 
cr ec er ef 
! (1) | (2) ‘TOTAL NUMBER OF EACH STATION ! 
1 TTEM +--+ 
! ID TAG ! (3) STATION ID TAG ! 
rr nna t 
\// 77 ! ! ] ! ! 
ee fp nf a 
Vi five | ! ! ! ! ! 


ee pe ee ef rt 
! LRUO1 WIS73.. 11s! 
! LRUO02 '14224.0 ! 
! LRUO3 !1267.49 ! 
! LRUO4 140524.4 ! 
! LRUOS5 '15808.4 ! 
! LRUO6 12157123! 
! LRUO7 '2006.86 ! 
! LRUO8 WNS773.2 
! LRUO09 159853.6 ! 
! LRU10O '38447.1 ! : 
tT 


eam Om Om Om == Cop 6 CoP Cop 6 == 


2 


*AVERAGE NUMBER OF DEMANDS 

ft 
! AVERAGE NUMBER OF DEMANDS DURING RESUPPLY TIME / TIME PERIOD T! 
4 ee + et 


ae Ge ! (2) TOTAL NUMBER OF EACH STATION ! 
{TEx 
! ID TAG ! (3) STATION ID TAG ! 
a ff nn ee 
Tay // | ‘ge! ! ! ] ! 
tf fff f+-———— 4 
way // tod ! ! ! ! ] ! 
tp rt ne np 
! TRUO] 1 0.49 ! ! ! ! ! ! 
| TRUO2 seis! ! 
! TRUO3 ! 0.46! ! ! ! 
! TRUO4 1! 14.59 ! ! ! ! ! 
! TRUOS I 556e! ! ! ! ! 
! TRUO6 1 7.69 ! ! ! ! ! 
! LRUO7 ony 2 a! ! ! ! ! ! 
! TRUO8 1 5.68 ! ! ! 
! LRUOS ! 21.55 ! ! ! ! 
! ! ! ! ! ! ! 


: LRU1O 1c 4e 
So a 8 ne 


i Foe Ae Racer 5 


| | 
WEES IMI] «occ ccc cee cc ccccees 2/7000.00 ! 
| | 
! SUM OF BACKORDERTIME......0e- 6805.2935 ! 
! AVERAGE NORS DURING T........ 18.9036 ! 
) ENDOORANCE FACTOR... cc cccccee ! 
| | 
! NUMBER OF BACKORDERS........2. 50.8480 | 
\ | 
! PROB. OF NO BACKORDER........ 0.000000 ! 
' PROB. AT LEAST ONE BACKORDER. 1.000000 ! 
| 1 
$a 


iS 


*ALLOCATION OF SPARES 
ae 


STOCK LEVELS PER ITEM AND STATION | 
nt 
! (1) ! (2) ! (3) TOTAL NUMBER OF EACH STATION ! 
| TYEM | 'TOTAL+--~--~-----------_---_------_--_---___------- + 
! ! NUM-! (4) STATION ID TAG 
! IDTAG ! BER! ! 


Bn pe ee ee ee ee pet 
PISS SLL ! ! ! ! ! ! ! ! : 
! a af ef af 
SSL LST SA te : ! ! ! ! ! ! ! ! 


Be a i er ee et 


! IRU01 ! Peles ! ! : ! ! : 
! ILRU02 ! 30 oe: ! : : ! ! ! : 
! LRUO3 ! lie ! ! ! ! ! ! 
! LRU04 : 5 ie ! ! ! ! ! ! 
! LRUOS5 ! 2 acer ! ! : ! ! ! ! ! 
! LRUO6 ! 4! 4! ! ! ! ! ! ! ! 
! LRUO7 : Do eee ! ! ! ! ! 
! LRUO08 ! 2) 25 ! : : ! ! ! ! 
! LRUOS ! Gan oe: ! ! ! ! ! ! 
! I1RU10 ! 2 lees ! ! ! ! : ! ! ! ! 
ap ee eee Le ee] tt 


*EXPECTED NUMBER OF BACKORDERS /ITEM & STATION 
fn 


! EXPECTED NUMBER OF BACKORDERS PER ITEM AND STATION ! 
of 
! (1) ! (2) TOTAL NUMBER OF EACH STATION ! 
{0 [TEM +-=----- st 
! ID TAG ! (3) STATION ID TAG ! 
a 
V/s a ee! ! ! ! ! ! 
a 
OC) re ! ! ! ! 
fp ee pe et 
! IRUOl  ! 0.1043 ! ! ! ! ! ! 
! TRUO2~=s !:- 2.2780 ! ! ! ! ! ! 
! TRUO3 = !-:«0.0899 ! ! ! ! ! ! 
! TRU04 = !:-:« 9.5903 ! ! ! ! ! ! 
! IRUODS) = !:« 3.7170! ! ! ! ! 
! IRUO6 ! 3.7676 ! ! ! ! ! 
! IRUO7 ! 0.2080 ! ! ! ! ! 
! IRUD8 ~=—s !:« 3.7046 ! ! ! ! ! ! 
! TRU09) =—s«: 115.5473 |! ! ! ! ! ! 
! IRUIO =—_:!11.8410 ! ! ! ! ! ! 
wo pe a pee fe ee pet 
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* ENDURANCE MEASURES /ITEM 
pt st 


(1) fh (2) te (3) (4) ES) 

!!EXPECTED ! PROBA- !EXPECTED ! INVEST- ! 

ITE !! SUMOF ! BILITY !NUMBER OF! MENT $! 

ID TAG !!BACKORDER! OF NO ! BACK- ! ! 
! 


!! TIME ' BACKORDER! ORDERS  ! ! 
|} 


! 

! 

! 

! 

+” 

! TRUOl = !! 13.0 !0.911514 ! 0.1043 ! 1000.0! 
! IRUO2  !! 249.1 !0.248497 ! 2.2780 ! 3000.0 ! 
! IRUO3 ss $!— «11.2 !0.922747 ! 0.0899 ! 1000.0 ! 
! IRU04 = !!-_—«:1196.1 !0.003708 ! 9.5903 ! 5000.0 ! 
! IRUOS !! 492.3 !0.077263 ! 3.7170 ! 2000.0! 
! TRUO6 =!! 408.1 !0.118566 ! 3.7676 ! 4000.0 ! 
! IRUO7 !! 26.4 !0.836348 ! 0.2080 ! 1000.0! 
! IRUOS = }!— 490.4 !0.077958 ! 3.7046 ! 2000.0 ! 
! TRUOQ =s- !!——- 2069.4 !0.000083 ! 15.5473 ! 6000.0 ! 
! IRU1IO !! 1849.4 !0.000108 ! 11.8410 ! 2000.0 ! 


i er i 


*ENDURANCE MEASURES /STATION 
ee 


Meme (2)  § (3) !. (4) 3! (S) ! 
'STA-!!EXPECTED !| PROBA- !EXPECTED ! INVEST ! 
'TION!! SUM OF ! BILITY !NUMBER OF! MENT |! 
! ID !!BACKORDER! OF NO ! BACK- ! | 

| 


'TAG !! TIME ' BACKORDER! ORDERS ! J 
et ne Fe ee 
'Toc1!! 6805.3 !0.000000 ! 50.8480 ! 27000.0 ! 
(pe a ee Le a aac eee 


*MARGINAL COST EFFECTIVENESS (AVERAGE NORS) 
4 
!ME/C IS THE CHANGE IN AVERAGE NORS (AVERAGE OVER TIME |! 
!PERIOD T) ACHIEVED BY A UNIT INCREASE IN INVESTMENT. ! 
!MC/E IS THE CHANGE IN INVESIMENT ACHIEVED BY ALLOWING ! 
!NORS TO INCREASE ONE UNIT. ! 
1 


: ME/C = 1/(MC/E) : 
nr cr ce reenact naemeteneenmmeeminchs 
' (ay) t! (2) $03) E44) OS) ES) 


1 
! V 'EXPECTED ! AVERAGE !MARGINAL !MARGINAL ! 
'POINT!! INVEST- ! SUM OF ! NORS !COST PER! EFFEC- ! 
! NO. !! MENT ~ ! BACKORDER! (DURING ! EFFEC- !TIVENESS ! 
! V ! TIME ! TIME $!TIVENESS !PER COST ! 
! V ! !PERIOD T)! (MC/E) ! (ME/C) |! 
pe 
! 1 !! 27000.0 ! 6805.3 ! 18.9036 ! ! 
Mp 
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*MARGINAL COST EFFECTIVENESS (BACKORDERS) 
et 
! ME/C IS THE CHANGE IN NO. OF BACKORDERS 

! ACHIEVED BY A UNIT INCREASE IN INVESTMENT. 

! MC/E IS THE CHANGE IN INVESTMENT ACHIEVED BY 
! ALLOWING THE NO. OF BO TO INCREASE ONE UNIT. 
: ME/C = 1/(MC/E) 

a 


1) '! (2) 2 eye ea ee) 
! 1 'EXPECTED !MARGINAL !MARGINAL ! 
!POINT!! INVEST- ! NUMBER !COST PER ! EFFEC- |! 
! NO. !! MENT !OF BACK- ! EFFEC- !TIVENESS ! 
i ! ORDERS !TIVENESS !PER COST ! 
VW ! (MC/E) ! (ME/C) ! 
ae a reef 
! 1 !! 27000.0 ! 50.8480 ! 


SS 


*FINAL C/E-CURVE (SBT) 
pe ee See fe pe et 


1 (yt! (2) -! (3) !$ () =! «(S) & (Cel Ge 
! 1! !EXPECTED ! ! PROBABILITY  !EXPECTED ! 
!POINT!! INVEST- ! SUM OF !ENDURANCE! ! OF AT !NUMBER OF! 
! NO. !! MENT !BACKORDER! FACTOR! OF NO !LEAST ONE! BACK- ! 
! i ! TIME ! ! BACKORDER! BACKORDER! ORDERS ! 
po ee | | ee ef a 
! 1 !! 27000.0 ! 6805.3 ! !10.000000 !1.000000 ! 50.8480 ! 


nt rr rt at 


END OF OUTPUT OPUS 8 
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